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Berechnung des Niederhalterdrucks beim Tiefziehen. 
Von H. Jung, Stuttgart. 
Mit 6 Textabbildungen. 


Zusammenfassung. Mit den in Teil I gewonnenen Ergebnissen wird gezeigt, daB die bisher 
vernachlassigte Aufstauchung der Platine einen wesentlichen EinfluB auf den Niederhalterdruck 
beim Tiefziehen hat. Im Anschlu8 daran wird die Bedingung fiir die Faltenbildung beim Tief- 
ziehen mit vereinfachenden Annahmen aufgestellt. Der zur Berechnung des Niederhalterdruckes 
benétigte Rechenaufwand wird skizziert. Die ermittelten Zahlenwerte fiir den Niederhalter- 
druck stimmen mit den Versuchswerten gut tiberein. Die bei der Berechnung des Niederhalter- 
druckes gemachten Annahmen wurden durch Versuche jeweils iiberpriift. 


Summary. The results derived from the contents of Part I are such as to reveal that the 
hitherto neglected upsetting of the flat bar exercises a considerable influence on the retaining 
pressure in deep drawing. In connexion thereto the conditions for the formation of folds in deep 
drawing are set forth on basis of simplified assumptions. The method of calculating the retaining 
pressure is shortly discussed. The values ascertained in this manner for the retaining pressure 
are in line with those obtained by tests. The assumptions on which the retaining pressure is 
calculated are in each and every case checked by adequate tests. 


Résumé. Les résultats obtenus dans la Premiére Partie démontrent que le refoulement, 
jusqu’aé présent négligé, de la plaque exerce une influence considérable sur la pression de retention 
dans les travaux d’emboutissage. Profitant des nouvelles notions auteur formule les conditions 
concernant la formation de plis au cours de ’emboutissage a la base de données simplifiées, en 
indiquant en méme temps la méthode de calcul qu’il recommande pour se rendre compte de la 
pression de retention. Les valeurs de la pression de retention établies de cette maniére coincident 
assez bien avec les résultats des essais. Les suppositions a la base desquelles la pression de retention 
est calculée sont contrdélées par des essais en chaque cas particulier. 


I. Problemstellung. 


Gegeben sei eine runde Blechscheibe vom Durchmesser 27, und der Dicke h. 
Aus dieser Blechscheibe soll ein Napfchen vom Durchmesser 27, und der Wand- 
stirke h hergestellt werden. Der Defor- 
mationsvorgang wird mit dem durch 
Abb. 1 in vereinfachter Form gegebenen 
Werkzeug durchgefiihrt. 

Bewegt sich der Stempel S nach 
unten, so wird das Blech (Platine) in die 
Bohrung eingezogen. Dieses Hinziehen 
fiihrt in dem auf der Matrize aufliegen- 
den Teil zu Radial- und Tangential- 
spannungen. Uberschreiten diese Span- 
nungen einen kritischen Wert, so tritt 
Faltenbildung auf. Beim weiteren Hin- 
ziehen der Platine miissen diese Falten 
wieder ,,ausgebiigelt‘‘ werden. Dieses Abb. 1. Werkzeug. 
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Ausbiigeln vergroBert die Stempelkraft und ergibt eine Verschlechterung der Ober- 
flachenbeschaffenheit des Werkstiickes. Setzt man auf die Platine einen Niederhalter 
mit der Kraft Py auf, so 148t sich, sofern Py einen Minimalwert tberschreitet, die 
Faltenbildung, die von der Biegesteifigkeit der Platine und dem Ziehverhiltnis r)/r, 
abhingt, verhindern. 

Es soll nun versucht werden, mit der in Teil I aufgestellten Theorie einen Hin- 
blick in die Vorginge beim Tiefziehen zu gewinnen. In Ziffer If wird angenommen, 
daB die Platine wihrend des Ziehvorganges ihre Dicke h beibehalt und die an der 
Ziehkante auftretenden Biegespannungen auf den Niederhalterdruck keinen EinfluB 
haben. Es lassen sich dann die Spannungen und Dehnungen naherungsweise be- 
stimmen. Aus den Spannungen kann dann der Niederhalterdruck berechnet werden. 
Hierbei zeigt es sich, daB die Voraussetzung der konstanten Dicke h zu auBerordent- 
lich hohen Niederhalterdriicken fiihrt, die technisch nicht realisiert werden kénnen. 
Ziffer III gibt eine Naherungstheorie zur Bestimmung der Beulgrenze, wobei eine 
Aufstauchung der Platine zugelassen wird. Der auferordentlich umfangreiche 
Rechnungsgang zur Bestimmung des Niederhalterdruckes wird dann in Ziffer IV 
skizziert. Fiir das in Ziffer V angegebene Zahlenbeispiel wurde die Rechnung durch- 
gefiihrt und die Endergebnisse angegeben. Es zeigt sich hierbei, daB eine gute Uber- 
einstimmung mit den mir bekannten Versuchsergebnissen erreicht wird. 


II. Ziehvorgang bei konstanter Dicke h. 


Die Platine wird entlang des Kreises r= 7, aufgeschnitten (Abb. 2) und dort 
eine Verschiebung c, vorgeschrieben. Es wird nun angenommen, dai diese Ver- 
schiebung c, klein von erster Ordnung gegen- 
tiber den sonstigen Abmessungen der Platine 
Cy eee! ba ist. Am AuBenrand der Platine ist die Radial- 
spannung a ee 


Es mu in den Punkten 


| r=, : | r=T, 
| 


die Druckspannung 
Abb. 2. Werkstiick. 


OPA, 


angenommen werden, wie aus den Gleichgewichtsbedingungen hervorgeht. Die 
Randbedingungen des vorliegenden Problems lauten in gemilderter Form 


w= @) fiir bond 0) und a h im Bereich Yo = i = Tas (1) 
b 
\o,dz=0 fir r=r, (2) 
0 
b a 
jude=—ce fir r=7%, (3) 
0 
o,= 0) firz=h, rer und r=r, (4) 


Laut Voraussetzung wurde cy klein von erster Ordnung angenommen, so da die 
in Ziffer VI, Teil I, aufgestellte Theorie zur Anwendung kommen kann. Die Losung 


der Differentialgleichung (37); wird aus zwei Teillésungen aufgebaut, wobei die 


Grundlésung ®,=c,Inr (5) 


die Hauptverschiebung und 
®, = sin Az [aS 9 (tAr) + 4,No(tAr)] + cos dz [b, Jy (t Ar) +b, N, (iAr)] (6) 


: 


a ee ee Se eS ee Te ee ee em 
: 
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die drtlichen Stérungen beriicksichtigt. Mit (5) und (6) erhailt man nach Ziffer VI, 
Teil I die Verschiebungen und Dehnungen 


u=—tAsin Az [a,J,(tAr) + a,N, (tar)]— } 
— bcos Az [b) Jy (i Ar) + by Ny (GAr)] +, 


w= Acosiz[a, J kana’ o(t@Ar)]— Asin Az [b, J) (@Ar) +6,N, (tar), 
e, = Msin dela, (J (6.47) — ge Ix (Ar) + ay(Q (iar) — Ge, Ar))] + 
+ eos 2 z|b, (Far) — 4, (iar) + 


+b, (No (iar) — > N, Ar))]—%, 


fe sin Az|5 - J, sed + 


+ 72 cos A2|o J, (tar) + 


oy pe? 


N, (i Ar)| + 
&,= — MsindAz[a,J,(tAr)+a,Ny(tAr)]—A2 cos Az [b,J (tar) +b,N,(tAr)]. 


Setzt man den in (6) auftretenden Parameter der Zylinderfunktionen 


so fiihrt die Randbedingung (1) auf 
G, =a, = 0. 
Die Randbedingung (3) ergibt, da hierbei die Storfunktion ®, bei der Integration 
verschwindet, Cy 
i= | 
Ks ist nun zweckmiBig, an Stelle der Bessel-Funktionen J, und N, die Hankel- 
Funktionen H," und H,) einzufiithren, um einfache Niherungsausdriicke fiir die 
Spannungen ableiten zu kénnen. Mit den Hankel-Funktionen und (7) erhalt man 
die Spannungen in der Form 


cos Az by (Ho (iar) — ap HL Ui ar)) fs 


=) 


S deanery ee 


op= {St + Boos Az |e H, (ar) Hy (Ar) 9 (2K) +H. : 
o, ={— # cos Az [b, Hy (Ar) + b,H™ (i Ar)}}g (|/2 /K ) + xo. 
Geht man mit der ersten Gl. (8) in (2) ein, so wird 
9/2 VK ), = %- (9) 


In (8) und (9) ist %) eine Konstante, die nach Ziffer VI, Teil I, den hydrostatischen 
Spannungszustand beriicksichtigt. Die in (9) auftretende Funktion K ist 


K= he Jy = 5 (& +6 + &)® — & Ey — &p & — Ey Ey. (10) 


Fiir die Durchfiihrung der Rechnung ist es nun zweckméfBig, an Stelle der in (7) 
gegebenen Dehnung nur die aus (5) folgenden Dehnungen 


a — Cc ea 
g= — 4, Ege  & = 0 (11) 
in (10) einzusetzen. 
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Geht man mit (8) in die Randbedingung (4) ein, dann erhalt man 


Sie ) 
UNS ts H,®) (i Ar) —Hy@) (é Ar») 
b, g (V2 VK) r=Te 
a es 2 Hi Gre. ) Hy(®) (6 Ary) Ho) (Ar) Ho (GA Tg)’ a2) 
gly? Vilar H, (é Ar) —Hy® (¢ 479) 
jp a ee g (V2 VK)r=n 
2 Fay 2 ae (i274) Hy® (i Ar) —H,® Ar») Ho) (EArg) J 


Beachtet man, daB Ar 1 ist, so kénnen an Stelle der Hankel-Funktionen ihre 
asymptotischen Anniherungen 


(Gy ae —Ar CU rice (2) wg r 
Hi, wil|/re » Ay yas , HO wil sae 
in (12) eingesetzt werden. Damit geht (12) tiber in 
g We \K)r=rq V/ 2 aha if 2 its ) 
b Cy g (V2 VK)r =1o TAs | TAT) 
1 Ue Sa ? 
ae re ae JEEP ie? (tq—T)__ 9 ("q— 10) 
A Tete 
BAG TE, (13) 
g (V2 VK)rarq 2 eta /_ 2 ge? % 
E Cy g (V2 VK )r=1% TAT y AT, 
Sim Orne Tat 2 I 2 
eer (a=) s_ ef ("a —70)) 
wah Tare J 
Die Ausdriicke (13) lassen sich noch weiter vereinfachen, wenn 
e-4*"g—%) | gegen = e* (7g) 
vernachlassigt wird. Man erhalt dann . 
b ey A Cy g (V2 VK) T= Tey 1 Vr 
iy tr, 2 a 
ils VK)r=m V2nd (14) 
—Z —Ar 
by & oe s Wet ea 


Geht man mit (14) in die drei Gl. (8) ein und ersetzt dort ebenfalls die Hankel- 
Funktionen durch ihre asymptotischen Annaherungen, so wird 


2VK)r= 
— 72 a a ee A(rmy—r) 
CG, =i cos 2 | 7 (V2 VEros, ~e nr) 4 
+x) 2 Hea) 9 (V2 ‘K) +x. (15) 


Die Gl. (15) ergibt nun einige wesentliche Aufschliisse fiir die Berechnung des Nieder- 
halterdruckes. 


Die durch die Stéorfunktion bedingte Spannung 


MELE On ern >) 
SWiyEuan Ve" of Berea] 9 VIVE) ate (18) 


hat keinen KinfluB auf den Niederhalterdruck, da sie auBerordentlich rasch abklingt. 
Der Niederhalterdruck hangt im wesentlichen nur von dem konstanten hydrostatischen 
Spannungszustand x) ab und ist 


g(/2/K),_ —r) [kg]. (17) 


Ge #t00812| 


— toe 


Py= 
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Aus (17) ergibt sich der sogenannte spezifische Niederhalterdruck zu 
Pyr— ad (/2|/K),.,, Use/em*). (18) 


Rechnet man mit (17) und (18) das in Ziffer V gegebene Zahlenbeispiel durch, so 
zeigt es sich, dafs diese Niederhalterdriicke mit den im Versuch ermittelten nicht 
tibereinstimmen. 

Wird eine Platine nur etwa ein Drittel in die Matrize eingezogen und dann aus 
dem Werkzeug entnommen, so lassen sich die Aufstauchungen im entlasteten Zustand 
messen. Es zeigt sich hierbei, daB die 
Dehnungen «, an der Ziehkante ein Maxi- 
mum haben, das nach aufen abklingt. 

Die bei der Rechnung gemachten An- 
nahmen stimmen mit der Wirklichkeit 
nicht tiberein, zeigen aber, da8 es nur 
darauf ankommt, die ,,Hauptverschie- 3] ae ey — 
bung’ zu erfassen, da nur diese einen “* ~ 
wesentlichen Einflu8 auf den Nieder- 
halterdruck hat. 


Z 


Up 


LL EEE 


lil. Aufstellung der Beulbedingung. 


Aus den Versuchen ergibt sich, daB Abb. 3. Ziehvorgang. 
der Niederhalter nur in dem Bereich 
ry Sr Sr; auf der Platine aufliegt (Abb. 3). Der AuBenrand der Platine 7; Sr Sr, 
ist frei und bildet Falten, sofern die Spannungen o, und o, einen kritischen Wert 
tiberschreiten. Nach dem Gesetz von RoS-Hichinger gilt nun 

VE (o; — oo)? =f /Z(ei — 0)? (6 =1, 2, 3). (19) 
Setzt man den Niederhalter auf die Platine auf, so tritt auBer der Radial- und 
Tangentialspannung noch eine Druckspannung auf. Durch diese Druckspannung 
wird o, und o, im Bereich r, =r <1; herabgesetzt, so daB auf diese Weise die Falten- 
bildung verhindert werden kann. 

Um den Bereich, auf dem der Niederhalter aufliegt, zu bestimmen, mu8 die ,, Beul- 
bedingung“ fiir eine auf dem Innenrand r = r; eingespannte Kreisringplatte (r; S 
<r<r,) aufgestellt werden. Bezeichnet man die Radialverschiebung auf dem 
Kreis r =r, mit c, so ist [Teil I (66) und (67)] 

ee ae? 2 | 


ya Bay L 
ha EN ie maa Ty 
1 (2% +1—/2xe+ 1)r+e 
Bre irae ae re =H 
hi 2 V2%+1r—e 


fhe 


1 / ——— ¢ 
po=5(1— /2xt1+5], 


: 
a 


~3(5—*), (21) 
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Vernachlassigt man noch +e gegen c, so erhalt man aus (21) 


/2/K= ee + [5 Ko) 7? +6 (1+ Ko) (22) 


r (r2 —¢) - 


Mit der ra des Wurzelausdruckes aus (22) wird in erster An- 
naherung 


2/K w faite 2 
2 VK y2 r? (72 —¢) ° (23) 
Die Verfestigungsfunktion nimmt damit die Form 
rae 
a 2 
t(y2 VK) = 4,4 B; fe 72 (r2 —c) or (24) 


an und es wird 


a2 VK )= AF OR ag 


bf 
2 


Um das Ausbeulen der Kreisplatte zu behandeln, macht man fiir die Dehnungen 
und Schiebungen den Niaherungsansatz, wenn mit £ die Ausbiegung der Platine 
bezeichnet wird: 


1 xrte arg ) 
Sr tee att @ Gp2? 
ie 1 res a (= ars 1 S | 
PI Ie 2 By? ex er |? 
ye pede (25) 
eee IG: 
1 (4 1 | 
Qe ROA ig or dy ean (enna. 


Die Annahme, daB die Verfestigungsfunktion nur von den Grundverzerrungen é,9, 
€y) und é,, abhangt, fiihrt auf 


; lxr+e aE \ / x; y2 r? (72 — ¢) 
a =(—5 72 —¢ + z | Ic| a Cc cere | 
rt + —— 
2 
; 1 c—xr 1 oe 1 a&\) / a, V2 72 ( ptt 
oy = |p tel Se +5 wll we ote rg a 
v2 > 
(26) © 
Diane ape { 
es @ Ve ann: 
72 “+ oy 
1 arg 1 a&\[ V2a,; 12 (r?—c) 
re *e Gap oe a wae ge celle) 
2 +- pest | 
Aus (26) ergibt sich der hydrostatische Spannungszustand zu 
‘ SN 1 Pt, 6 
Oj =e ier iy r sia dr + A,. (27) 


Verwendet man bei der Reihenentwicklung der Verfestigungsfunktion nur die beiden 
ersten Glieder und wertet das Integral (27) aus, so wird 
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er er es x; G Oy ea ) 
Fo. (1 + x) =In|r — *_In it 
e| 2/2 2 aes tac z 
c B; ate a 
Pah re ae eet a a et aes 
1 c—x#r lxrte V2% 12 (2+ Je 
Fo =(5 r VAN hee |c| L D B+ rm 
teen eae 
2 
1 xr?te\|(V/2 a, 7? (r? + |e 
a0 = (Ko si 212 S76 |c| ‘2 {oh ! + Bi| + ore, 
72 — . 
SE Ee Se Se Ee ed 
7a =| Foe r2 Ay? r Or ss r all lc| Apa e? t Bi ae (28) 
_d 2 
+ A, 
LOE 1 0 aE\ [V2 % 2 (r2 + |e|) 
oe og? ard or = || , C t Bi] + Om, 
by 


In (28) wurden mit oj) (i = r, y, z) die durch den Grundverschiebungszustand hervor- 
gerufenen Spannungen bezeichnet, wahrend o;, die zusitzlichen Spannungen sind, 
die in der Platte Biegemomente hervorrufen. 


Aus (28) ergibt sich, wenn h die Dicke der Platte ist, das Radialmoment 


Ha oe tale +6, (29) 
72 — 5 


i ls HS le 
rn) Or r or og? 


das Tangentialmoment 


M, a Cc 
pe ey 


2 


me pl. ee Woe eae Qo, 9? (2 + 
( a AE Ba pt A +a (30) 


das Drillmoment 
are é 2 2 (72 
eam eakale afin ia th (ai) 


12\r Gép r2 Op Pa ee 


2 


Leis. a [ig Ie BE 1 ae (array hee 
A= x \ \ [ae + Ui, (= dy? r r aS : 2H, (= “Orép or ae)|r ar dy. (32) 
Geht man mit (29) bis (31) in (32) ein und setzt 
E=sinmg:¢(r), (33) 
so erhalt man 
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Ta 


4n— 3 (pe (S+1E— Moraes SPEEA) 


r ar r Ts 


r dr ip Paes ie 


a2l 1 dt m? m dl m ..\2) | V2 0, (2+ [el) L 2. 34 
— + .. )+a (eae) [Me + B,\rdr. (34) 


Wie aus Versuchen bekannt ist, wird die Faltenzahl 
m > 10. 


Damit kénnen in (34) die Glieder mit m”, wobei n < 2 ist, vernachlassigt werden. 
Gl. (34) geht mit dieser Vereinfachung tiber in 


VICE wm one \ 
Ae = ea Fs - BS 2 [ 7 IONE ee | | r dr, (35) 
ae 
wobei zur Abkiirzung gesetzt wurde: 
1 dé 2 V2 x, 12 (r2 + |e}) dr 
us= ||, dr re | Me n—|s B r 1 6. (36) 


Die in der Platte aufgespeicherte potentielle Energie ist 


Apa “El (aera Gs (37 


Ts 


4 


Beriicksichtigt man in (37) nur die Glieder mit m?, dann erhalt man 


Ta 


5} \ Gyo ar. (38) 


Die Minimalbedingung fiir die Formanderungsarbeit fiihrt mit (35) und (38) auf 


Ta 


272 C - 2 ‘V2 a0, 7 Bie 5 ce ‘ 
oe 53 za (Tet Cale a ae | % dr = Min. (39) 


r 


Aus (39) laBt sich ni&iherungsweise mit Hilfe 
des Ritzschen Verfahrens der Eigenwert be- 
rechnen. Hierbei ist zu beachten, daB die 
Integrationsgrenzen 7; und 7, ebenfalls noch 
von der kritischen Verschiebung c abhingen. 

Wird die Verfestigungsfunktion entsprechend 
Abb. 4 angenommen und durch einen Polygon- 
zug angenihert, so ergibt sich 7; aus der 


xleichung 
/9 r2(r2+ |e 
= (o; — O41) sete eb + Bi — Biss = 9. (40) 
els 
Fir 7, erhailt man 
Vg — Ta 18k (41) 
[eea-d : 


wobei 7, der urspriingliche AuBendurchmesser 
Abb. 4. Verfestigungsfunktion. der Platine ist. 
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Um die Minimalbedingung (39) zu befriedigen, wird der eingliedrige Ansatz 
C= ay (r — 14) 
verwendet, der die geometrischen Randbedingungen 
de 
f=; ee = 0 
erfiillt. Geht man mit (42) in (36) ein, so laBt sich das Integral auswerten. Es er- 


gibt sich = eas 
M / 3) / Te (+V2)e-V 3) 
es ? V2 «, ; - a Cc a 
ga ard tng) gn (2) (n+ |.) 


: Tr “a 


Y 
2 pe a 
file ie 


r2 : In ic In e +. 2B; l27 


Die in (36) auftretende Konstante wird so ermittelt, daB die Randbedingung 
Pat. Me 0 


a 


erfiillt wird. 
Mit den dimensionslosen Verinderlichen 
r |c| r; jit 
aan = X, 7.2 pan! Fn Xs ra = Xy 
a a a a 


erhalt man 


— i) parlecty 
i ere = 4 («+//$) (= lee. 


='2 L;| ¢©— 2x, In = —— |— 
a q co Dy, (2— 12 (= a 2 
é : ®, 5 
3 rs 2 2 
; 2 v2 ean G2 — & Gon 
i id 2 2 In ea q In 2 ++ 2B, (2a2 ee Gq? a as x ), (44) 
t_° — > 
xt — 2 ‘ 
Opo= = et In A 2 In— ( tant 4 
\ 2 tq — > a as 
CR alla CUE Sl a Bee a) 
r 2 |2 is (%4" + q) x? r 2 a a” : x? (a + @Q) s f) 
a) 3K, 22—2K 
; — Z x C2 — —— 
Sg = | In eget raet med) 
V2 vq? — a a aia 
Be {a (a? + q) «,? Ca x errs 
1 2 LS % In (ae g) e- : Peitihd pine Aus e+ g (46) 
Die Minimalbedingung fiihrt mit (44) bis (46) aut 
va d 
Ti | Fg (@ — a) — 
m?* h? gargs : Bett ate : ie: (47) 
Linas va F 
ns (w—wa,)? (c— @,)* | V 2a; x? (x? + Q) +B, ee 
ay ao? oe q wes ; 
“ 2 
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Um aus (47) den kleinsten Eigenwert zu ermitteln, wahlt man ein bestimmtes q und 
wertet die Integrale bei festem «; und f; aus. Es ergibt sich dann ein Wert 


m? h? 
fees 
Wird nun q variiert, so erhalt man 
1 =F (9). (48) 
Der Eigenwert ergibt sich dann aus der Bedingung 
a — Min. (49) 
q 


Bei der Berechnung des Niederhalterdruckes ist zu beachten, daB auch die hoheren 
Kigenwerte vermieden werden miissen. 


IV. Berechnung des Niederhalterdruckes. 
Mit Abb. 3 wird der Niederhalterdruck 


x. 
27,7 : 

PN = Grn \ 0,90 dx. (50) 
Xo 


Die untere Integrationsgrenze x, ist fest, waihrend x, von dem jeweiligen Deformations- 
zustand abhangt. Sie muB so ermittelt werden, daB der freie AuBenrand der Platine 
(x; <u <u,) keine Falten bildet. 


Wahlt man die Verfestigungsfunktion f (|/2\//K ) entsprechend Abb. 4, so wird 
eine erste dimensionslose VerschiebungsgréBe g so angenommen, dafi der AuBenrand 
der Platine gerade noch im rein elastischen Bereich der Verfestigungsfunktion ver- 
bleibt (x) = 0). Ebenso soll bei diesem ersten Deformationsvorgang die Aufstauchung 
im elastischen Teil der Platine vernachlissigt werden. Damit wird 


Mix IRC 2 Pics 0s “on | A 
dy = 2Bp |2 a, gta?  “§ ee |? (51) 
Bod [e2®—2? sees 59 
Co 4 | ata  ~ a®(a® + q)|’ ye 
: d 
ee at} x 
me Fo (% — fe 
m2 h2 uv, ~ 
Sr ietep aera se eee (53) 
[fat oat 0) ph gy 
on x x J 


Da es sich hierbei noch um sehr kleine Verschiebungen handelt, kann «, = 1 gesetzt 
werden. 

Die Rechnung wird nun wesentlich vereinfacht, wenn eine Grenzfaltenzahl mg 
festgelegt wird, die gerade noch vermieden werden soll. Man wiahlt nun drei Werte q,, 
dq. und q; und erhalt fiir diese drei Werte aus (40) 2;,, 2;, und x;,. Wertet man (53) 
numerisch aus, so erhéilt man m,, m, und ms. 

Der zu mg gehorige Wert gg kann nun auferordentlich einfach mit Hilfe der 
projektiven Unterteilung gefunden werden (Abb. 5). 

Man tragt auf einer Geraden G, die Werte m,, m., ms linear ab und wihlt einen 
Pol, zieht die Polstrahlen und paBt in das Strahlenbiischel eine Gerade G, so ein, 
daB gq, auf den Polstrahl nach m,, q. auf den Polstrahl nach m, und q, auf den Pol- 
strahl nach m, zu liegen kommt. Die GrdBen q,, g, und q, werden auf der Geraden G, 
ebenfalls linear aufgetragen. Der zu mq gehérige Polstrahl schneidet dann auf der 
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Geraden G, den gesuchten Wert qa, aus. Mit diesem Wert gg, bestimmt man das 
dazugehorige Wis Die in (46) auftretende Aufstauchung K, wird nun so bestimmt, 


da die Bedingung Oo=0 fir xwz=2 
20 m4 


erfiillt wird. Dies fiihrt auf die Gleichung 


By {1 ui (2,7 + gy) @,? Lq? — 2,2 ie 3 Koy ee 2% Koy— UT G, 
2 fa Ee ile a qq) we qi-9 v2 @,2 222+ 9; 
be ay 1— x, | ago . 2] x; th 3 Ko; «2 — 2 q, Ko; — UY 
Liye F n = 40 —— ty 54 
V2 2 eet a % (54) 
x, —— x,;2 — — 
7 BY a 9 


kh 
07 Q015 402 O05 PF 
se th Si at aes ee 4 
006 006i Gi? Qh Qi ale QQ Gith 


Abb. 6. Niederhalterdruck. 


Abb. 5. Ermittlung von q,. 


Mit dem aus (54) bestimmten Ay, geht man in (46) ein und kann dann das Integral (50) 
auswerten. 

Man nimmt nun die Verschiebungsgr6Be g beim zweiten Schritt so an, daB der 
AuBenrand der Platine im Bereich «,, 6, der Verfestigungsfunktion bleibt und be- 
stimmt auf dieselbe Weise wie beim ersten Schritt die kritische Verschiebung q,. 
Die in (45) auftretende Aufstauchung A, wird in erster Anniherung Ko,,= Ko, ge- 
setzt. Hierbei ist nun zu beachten, dai beim zweiten Schritt x, noch von g abhiangt. 

Es ist nun moglich, daB bei fest angenommenem mg und h 


t Bd = ae il. pps Fes = 


Ax: 4 
A xe x q ress 
2 


le (z—a,)* — (e—a,)4 | V2.0, 2 (x* + 9) re | ode 


ist. In diesem Falle bleibt die Platine noch stabil und konnte ohne Niederhalter 
eingezogen werden. Diese Erscheinung zeigt auch der Versuch bei groBen Platinen- 
dicken. Die Rechnung wurde an einem Zahlenbeispiel durchgefiihrt. 
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V. Zahlenbeispiel. 


Als Zahlenbeispiel wurde eine Platine mit den folgenden Abmessungen gewahlt: 
ry = 20cm, r, = 4cem, 0:02 <A S0°24cem. Fir M, (63) wurde die Verfestigungs- 
funktion bestimmt und durch einen Polygonzug entsprechend Abb. 4 angenahert, 
wobei 

Xt = 0, pa =~ 10-104 [kg/em-], 


wo, = 40-10, B,= 2°0- 104 [kg/em?], 
X, = 1731-102, B, = 33°3- 10? [kg/em?] 
angenommen wurde. Die Niederhalterdriicke sind dann fiir eine Grenzfaltenzahl 
Ng = o2 


ermittelt worden. Sie sind in Tab. 1 zusammengestellt und in Abb. 6 aufgetragen. 


Tabelle 1. Niederhalterdruck (ry) = 2, r, = 4 [cm)). 


| | | | | 


| | | | 
PARSIAN| 5 S60 oF 0°02, 0°04 | 0°06 | 0708} .0°10} 0°12) O-14) 0:16) O18) 0:20 


Py(ke/em®].| 75°6 | 621 | 53:3 | 876 | 260 | 21-4 | 155 |. 16-7 | 13:2°/9 os 
Py(kg]..... 2848°6 |2340 2008 |1417 | 979°7 80635 584-0 |629°2 497-4 358-0 


Es ist nun auSerordentlich interessant, daB von h = 0°15 cm ab die Verkleinerung 
des AuBenradius der Platine einen Einflu8 auf den Niederhalterdruck gewinnt. 

Die hier gewonnenen Ergebnisse stimmen mit den Versuchsergebnissen von Herrn 
Dr. BeiBwenger, soweit sie mir bekanntgeworden sind, gut tiberein. Dies diirfte 
darauf zuriickzufiihren sein, daB der in Ziffer III verwendete Ritz-Ansatz zu groke 
Werte liefert, die durch die Wahl einer Grenzfaltenzahl herabgemindert wurden. 
Um genaue Werte zu gewinnen, mtiBte die Rechnung mit den wirklich auftretenden 
Faltenzahlen, die nach Ziffer III ermittelt werden, nochmals durchgefiihrt werden. 
Da aber die an der Ziehkante auftretenden Biegespannungen nicht mitberiicksichtigt 
wurden, diirfte eine Wiederholung keine wesentlichen Verbesserungen mehr bringen. 


(Eingegangen am 4. November 1952.) 


Die Biegefestigkeit von Balken 
auf zwei Stiitzen aus bildsamen, spréden und Verbundwerkstoffen. 


Von K. Sehaden, Wien. 
Mit 6 Textabbildungen. 


Zusammenfassung. Fir Werkstoffe mit beliebiger Spannungs-Dehnungslinie o = f (e) und 
fir Verbundwerkstoffe (z. B. Stahlbeton) wird unter Beniitzung allgemeiner Beziehungen nach 
A. Nadai das Biegebruchmoment eines Balkens auf zwei Stiitzen mit J, 7 und Rechteckquer- 

; oe : dM 
schnitt berechnet. Dabei ist der Bruch bzw. FlieBbeginn durch ae. 0 definiert, wobei # den 
Neigungswinkel zweier Querschnitte bezeichnet. Der Bruchvorgang wird an einem unbewehrten 
und stark bewehrten Betonbalken eingehend erliiutert, wobei aus den gemessenen Dehnungen 


am Biegezug- bzw. Biegedruckrand und dem Biegemoment die Spannungs-Dehnungsbeziehungen 
berechnet werden. 


Summary. For building materials of an optional tension-elongation curve o = f(e) and for 
compound-type building materials, such as steel-reinforced concrete, the rupturing bending 
moment of a beam on two supports of 7, Z, and rectangular cross-sections has been calculated, 
utilizing the general relations, as established by A. Nadai. The point of rupture, or the yield 


point are defined by the equation Sean 0, whereby ® signifies the angle of inclination of two 
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cross-sections. The phenomenon of rupture is explained in detail with reference to a concrete 


beam, not reinforced, and a heavily reinforced concrete beam, the values of the elongations 
measured on the bending-tension edge and the bending-compression edge, and that of the 
bending moment being used for calculating the relations between tensions and elongations. 


Résumé, Pour les matiéres de construction d’une courbe quelconque tension-élongation 

o =f (e) et pour les matiéres composées, p.e. le béton armé, d’acier, le moment de flexion A 
> . . . . . 

rupture d’une poutre sur deux appuis & section en forme de I et de 7T', ainsi qu’& section 

rectangulaire, est calculé, en utilisant les relations générales selon A. Nadai. Le point de rupture 

ou bien la limite d’étirage est défini par la fonction Ss 0, & désignant l’angle d’inclinaison 


de deux sections transversales. Le phénoméne de la rupture est expliqué en détail par rapport 
& une poutre en béton non-armé et une poutre en béton fortement armé et les relations entre 
tensions et 6longations sont calculées en se servant des élongations mésurées au bord de la flexion- 
tension et au bord de la flexion-compression, ainsi que du moment de flexion. 


I. Einleitung. 


Wir betrachten den Fall der reinen Biegung (ohne Querkraft) eines Balkens 
aut zwei Stiitzen aus sprédem oder bildsamem Werkstoff, dessen Spannungs-Dehnungs- 
linie nicht linear oder linear (Hook- 
sche Annahme) verliuft und die Be- 
anspruchung im iiberelastischen oder 
plastischen Bereich liegt. Der Stab habe 
einen tiber die ganze Stabliinge gleich- 
bleibenden rechteckigen, T- oder I-fér- 
migen Querschnitt und sei durch ein 
Biegemoment in der Haupttriagheits- 
ebene senkrecht zur Querschnittsbreite 
b (Abb. 1) belastet. AuBerdem gelten 
folgende Voraussetzungen : 


Abb. 1. Doppelt bewehrter Stahlbetonquerschnitt. 


1. Die Bernoullische Hypothese, das heiBt die Langen der Faserstiicke 
zwischen zwei benachbarten Querschnitten verhalten sich wie deren Abstand vom 
Kriimmungsmittelpunkt (Ebenbleiben der Querschnitte, was im folgenden nach- 
gewiesen wird). Dadurch wird das Vorhandensein von Schubkraften im Querschnitt 
ausgeschlossen, was im Falle der reinen Biegung zutrifft. Nach K. Jager kann 
die strenge Giiltigkeit der Bernoullischen Hypothese bei reiner Biegung fiir ein 
beliebiges Formanderungsgesetz aus einem raumlichen Spannungszustand ab- 
geleitet werden, wobei nur kleine Formanderungen vorausgesetzt sind!. Durch sorg- 
faltige Versuche von C. Bach? und A. Féppl (nach A. Probst’) mit Stoffen, deren 
Forminderungslinie nicht dem Hookschen Gesetz entspricht und durch Versuche 
von E. Meyer? iiber die Durchbiegung von Staben aus FluBeisen nach Uberschreitung 
der FlieBgrenze wurden diese Annahmen bis weit in den FlieBbereich bestatigt. 
Aus den an Stahlbetonquerschnitten durchgefiihrten Messungen® ist zu ersehen, 
daB im rissefreien Balkenbereich sich die Dehnungen gentigend genau proportional 
zur Entfernung von der Nullinie verhalten. Die Dehnungen am Zugrand nach dem 
Auftreten von Rissen wurden an einer Mefstrecke gemessen, innerhalb welcher 


“pba erred) — Zur Bernoullischen Hypothese bei nicht linearem Formanderungsgesetz. Osterr. 
Bauzeitschrift 9, 135 (1947). . 

2 ©. Bach und R. Baumann: Elastizitét und Festigkeit, 9. Aufl., S. 259. Berlin: Julius 
Springer. 1924. j 

3 3 ‘A Probst: Vorlesungen tiber Hisenbeton, Bd. I, 8. 249. Berlin: Julius Springer. 1917, 

Siehe auch B. Kirsch und R. Saliger: Tetmajers Versuche, Armierter Beton 2 bis 4 (1911). — 
K. Hajnal-Konyi: Tests on concrete beams. Magazin of concrete Research 9, 113 (1952). 

4 7, Ver. dtsch. Ing. 8. 197 (1908). é 
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Risse vorhanden waren, wodurch das Ergebnis unsicher erscheint. Jedenfalls kann, 
wie die verglichenen Versuche in Abschn. E zeigen, zunachst das Ebenbleiben der 
Querschnitte auch fiir den gerissenen Stahlbetonbalken vorausgesetzt werden. 

2. Der Einflu8 der Querkontraktion kann vernachlassigt werden und die 
Querschnittsform bleibt unverandert. Aus den Ergebnissen der strengen St. Venant- 
schen und der daran anschlieBenden Biegetheorien ist zu ersehen, da die oben ge- 
troffenen Annahmen der technischen Biegelehre fiir alle praktischen Falle ausreichend 
sind’ (S. 177). Ausgenommen ist der Fall des Biegungstragers mit hohem Steg oder 
breiten Flanschen (s. Abschn. D). Diese Ergebnisse sollen auch fiir den tiberelastischen 
Bereich als giiltig vorausgesetzt sein. 

3. Die Spannungs-Formanderungs-Beziehungen o =f (e) in den Fasern 
des Querschnittes sind bei Biegung die gleichen, wie die bei einachsiger Zug- und 
Druckbeanspruchung. Die Zulassigkeit dieser Annahme kann nur durch Ver- 
gleich von Versuchen mit den theoretischen Biegewerten tiberpriift werden, da eine 
ausreichend genaue Methode der direkten Spannungsmessung nicht bekannt ist. 
Gemessen werden die einem Biegemoment entsprechenden Dehnungen am oberen 
und unteren Balkenrand. Die zugehérige Spannung wird aus den Formeln der Biege- 
theorie berechnet und die so erhaltene Forméanderungslinie mit der bei einachsiger 
Beanspruchung erhaltenen, verglichen. Nach Abschn. C stimmen beide Formanderungs- 
linien gut tiberein, so daB die oben getroffenen Annahmen als berechtigt erscheinen®. 

Mit den Voraussetzungen 1 und 2 ergeben sich die folgenden Abhangigkeiten 
von Formanderung und Querschnittsgr6Ben: Die spezifische Dehnung in der Ent- 
fernung 7 von der Nullinie ist mit den Bezeichnungen von A. Nadai’ bekanntlich 

q 
ate. (1) 
wenn o den Kriimmungshalbmesser der Nullinie des verbogenen Balkens bezeichnet. 
Der Winkel 0, um den sich zwei Querschnitte im Abstand / gegeneinander neigen, 
ist bekanntlich’ 
l 
Bezeichnet d, die Hohe des Querschnittes und 7,,, den Abstand der Nullinie vom 
unteren bzw. oberen Rand, so gilt 


m1 + Nz = do, (3) 
wobei nach Gl. (1). 
M1 = 0, e= @ &. (4) 
Aus Gl. (3) und Gl. (4) wird 
1 & +e 
Soe eh ; (5) 


5 J. W. Geckeler: Elastostatik. Handbuch der Physik, Bd. VI. Berlin: Julius Springer. 
1928. 

6 Fir den Bruchzustand liegt nach Untersuchungen an Glasstiben von A. Smekal die 
Biegefestigkeit eines Zylinderstabes mit spr6dem Stoffverhalten nicht unbetrachtlich héher als 
die Zugfestigkeit, da beim Biegestab nur etwa ein Viertel der Gesamtoberfliche eine bruchaus- 
lésende Inhomogenititsstelle lefert, wahrend beim gleichférmig gedehnten Zugstab diese Stelle 
von der Gesamtoberflache geliefert wird. — Wieweit diese qualitative Feststellung mit Riicksicht 
auf die unvermeidlichen Unterschiede in der Stoffgiite von Zug- und Biegeprobe quantitativ 
von KinfluB8 ist, kénnte mit Hilfe der hier abgeleiteten Biegeformeln geklirt werden. Die 
quantitativen Vergleichswerte von A. Smekal beruhen offensichtlich auf elastizititstheoretischer 
Grundlage. — Siehe A. Smekal: Zum Bruchvorgang bei sprédem Stoffverhalten unter ein- 
und mehrachsigen Beanspruchungen. Osterr. Ingenieur-Arch. 7, 49 (1953). 

7 A, Nadai: Der bildsame Zustand der Werkstoffe. Berlin: Julius Springer. 1927. Theory 
of Flow and Fracture of Solids. New York: Mc Graw-Hill Book Company. 1950. 
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und mit Gl. (2) 


U(e, + &3) 
6 Suen (6) 
Die Lage der Nullinie ist dann nach Gl. (4) und (5) 
apd, d 
Ro ana Ne Resear (7) 


Das Gleichgewicht zwischen den Schnittkraften M,Q — 0, N = 0 und den Spannungen 
im Querschnitt o, t= 0 verlangt 


| odq=0, (8) 
|ondg=M, (9) 

wobei in der Faser 7 des Querschnittselementes von der Breite 6 = ® (n) 
dq = ® (n) dy (10) 
die Normalspannung o = / (e) vorhanden ist. Mit Gl. (1) und (10) wird Gl. (8) und (9) 
[f (2) P(e e)de=0, (8a) 
M =o | fle) ®(e €) ede. ; (9a) 


Die kiassische Biegetheorie berechnet auf Grund der Hookschen Annahme iiber 
die Proportionalitat zwischen Spannung und Dehnung die Randspannungen und 
bezeichnet fiir den Bruchzustand diese Spannungen als Biegefestigkeit, welche sich 
als querschnittsabhangig und verschieden von den im einachsigen Zug- und Druck- 
versuch ermittelten Werten ergibt. Die Plastizitatstheorien nehmen eine zur ein- 
achsigen Spannungs-Dehnungslinie affine voll-* oder teilplastische® Spannungsver- 
teilung tiber den Querschnitt an und berechnen aus dem Moment der inneren Krifte 
das Bruch- bzw. FlieBmoment. Der Hebelsarm der Innenkrafte wird dabei nach 
Versuchen festgelegt oder auf Grund einer angenaherten Spannungs-Dehnungslinie 
rechnerisch ermittelt. Die Gro8e der beim Bruch vorhandenen Randspannung und 
Randstauchung wird von den einzelnen Verfassern unterschiedlich angenommen. 
Im Gegensatz zu diesen Naherungsansaitzen wird hier der Bruch bzw. FlieBen 
dadurch definiert, da eine Zunahme des Neigungswinkels zweier 


Querschnitte # bei gleichbleibendem Moment erfolgt, das heifbt = 0, 


wobei in den Randfasern die Spannungs-Dehnungsbeziehungen des ein- 
achsigen Versuches gelten. Eine Behinderung der Formanderung wird als ver- 
nachlassigbar nachgewiesen. Einen guten Hinblick in den Bruchvorgang ergibt die 
Untersuchung eines unbewehrten Betonbalkens (Abschn. C). Nach dem Erreichen 
der Bruchdehnung (gemessen im einachsigen Zugversuch bzw. berechnet aus dem 
Verdrehversuch) und Ri®bildung am Balkenzugrand, steigt das Biegemoment noch 
um 17% bis zum Bruch. Beim ausreichend bewehrten Stahlbetonbalken erfolgt 
nach dem Erreichen der Prismenfestigkeit im Druckrand ebenfalls noch eine geringe 
Zunahme des Biegemomentes bis zum Bruch. Wird die Stauchung am Druckrand 
nach dem Erreichen der Prismenfestigkeit bzw. der Bruchstauchung bis zum Hintritt 


8 Zum Beispiel W. Prager: Die FlieBgrenze bei behinderter Formaénderung. Forsch. Gebiete 
Ingenieurwes. 88, 95 (1933). ee 

9 KF, Fritsche: Tragfahigkeit von Balken aus Stahl mit Beriicksichtigung des plastischen 
Formanderungsvermégens. Bau-Ing. 11, 851 (1930). — Zum Beispiel K. Haberstock: Die n-freien 
Berechnungsweisen. Deutscher Ausschu8 fur Stahlbeton, H. 103. Berlin: Wilhelm Ernst und 
Sohn. 1951. 
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des Bruches weitergemessen, so ergeben sich die vielfach gemessenen hohen Werte 
der ,,Biegedruckstauchung“, die jedoch nur fiktive Bedeutung haben. Kine ein- 
gehende Untersuchung des BiegeflieBvorganges an Versuchen mit Balken aus Stahl 
mit ausgeprigter FlieBgrenze ist hier nicht méglich, da die bekanntgewordenen Ver- 
suche groBe Streuungen im Werkstoff des Biegebalkens und der Zug- und Druck- 
versuche aufweisen bzw. ausreichende Versuchswerte nicht angegeben sind. Jedenfalls 
tritt auch beim Stahlbalken nach dem Erreichen der FlieBgrenze an den Randern 
noch eine Zunahme des Biegemomentes bis zum Flie8moment ein, wodurch das in 
der Literatur erwihnte Auftreten von FlieBfiguren vor dem Erreichen des Flieb- 
momentes geklart erscheint. 

Bei Verwendung von Gleichungen nach A. Nadai wird aus dem Biegeversuch 
eines stark bewehrten Stahlbetonbalkens die Druck-Stauchungslinie des Betons 
berechnet, welche mit dem Druckversuch an einem Prisma tbereinstimmt. Ebenso 
zeigt die aus dem Biegemoment und den Dehnungen am Zugrand eines unbewehrten 
Betonbalkens berechnete Spannungs-Dehnungslinie gute Ubereinstimmung mit den 
im Zugversuch ermittelten Werten. 

Mit den allgemeinen Beziehungen von A. Nadai wird im Abschn. B die Bruch- 


bedingung = = 0 fiir den Rechteckquerschnitt abgeleitet. Das Bruchmoment des 


doppelt bewehrten I-Stahlbetonquerschnittes ergibt sich nach Abschn. D. Durch 
Nullsetzung der einzelnen QuerschnittsgroBen folgt das Bruchmoment fiir die daraus 
abgeleiteten bewehrten und unbewehrten Querschnittsformen aus beliebigem 
Werkstoff des Balkens und der Bewehrung, sofern die oben getroffenen Voraus- 
setzungen zutreffen. Fur den Stahlbetonquerschnitt kann durch Beriicksichti- 
gung der Betonzugfestigkeit o, = k,oz oder deren Nullsetzung o, = 0 die Mit- 
wirkung der Betonzugzone beriticksichtigt oder ausgeschaltet werden. Der Faktor k 
beriicksichtigt die Eigenspannungen aus Schwinden und Kriechen, die auch die 
Festigkeit am Druckrand o, = k, op und die Spannung in der Bewehrung beeinflussen. 
In Abschn. E sind die theoretischen Werte und die Ergebnisse von Versuchen mit 
schwach und stark bewehrten Stahlbetonbalken verglichen, wobei der Unterschied 
zwischen dem berechneten und gemessenen Bruchmoment héchstens 9°% betragt. 


II. Rechteekquerschnitt. 
Ks soll zunaichst die Bruchbedingung des Biegebalkens = 0 am Rechteck- 


querschnitt niher erlautert werden. Hier ist b = ® (j) = konst. Wird nach A. Nadai’? 
die Integralbedingung Gl. (8a) und (9a) nach # differenziert, so wird 


Gye, — 6, Oe, (11) 
und 
1 ee Gi 
b dg? mae = 5 Gp (MO). (12) 
Mit 
0 8 
Fi=|f (ede,  Fy=\f(e)ds (13) 
le 0 
lautet die Bedingung Gl. (8a) 
f,=Ff,, (14) 


Wird fiir ein bestimmtes Forminderungsgesetz o = f (e) die Summenlinie der Flachen- 
inhalte aufgetragen (Abb. 2), so entspricht nach Gl. (14) und Abb. 2 jeder Spannung 
am Zugrand — o, eine Spannung o, am Druckrand. Nach Gl. (8) muB auch bei einer 
Belastungszunahme des Stabes die Zunahme der Flicheninhalte dF, = 6, de, und 
dF, = o, de, entgegengesetzt gleich grof sein. Gl. (12) gibt die Abhangigkeit der 
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Querschnittsabmessungen, Randspannungen und Neigungswinkel # zweier Quer- 
schnitte vom Biegemoment an. 

Die Erschépfung der Tragfihigkeit eines Balkens soll dadurch gekenn- 
zeichnet sein, da eine Zunahme der Forman derung und damit des Winkels # 


bei gleichbleibender Belastung erfolgt, das heiBt a“ 0. Diese Bedingung 
fiihrt bei spréden Stoffen zum Bruch, wihrend 
sich Balken aus bildsamen Stoffen dabei weit- 
gehend verformen und so als Tragelement un- 
brauchbar werden. Die Voraussetzungen nach 
Abschn. A treffen fiir die bis zum Bruch der 
spréden Stoffe geringen Formainderungen noch 
zu; Bildsame Stoffe werden je nach der GréBke 3 
ihrer Formanderungsfihigkeit, Abweichungen 
von den theoretischen Werten zeigen. Mit 
= 0 wird Gl. (12) zur Bruch- bzw. FlieB- 
bedingung bat nam Ga er cia a 
oe a iar ae are ue) 
wobei o, und o, nach Gl. (14) bzw. Abb. 2 einander entsprechen miissen. 


Bei fiir Zug und Druck gleichartiger Formanderungslinie wird ¢, = ¢€,, 6, = 0, = 0 
und in Abb. | F,=F,, so daB Gl. (15) 


tO 


o-flel 


u=**, (16) 


die bekannte Formel fiir plastische Biegung wird, die sich dort aus der Annahme 
einer vollplastischen Spannungsverteilung tiber den Querschnitt ergibt. 


ILI. Ermittlung der Formanderungslinie in den Randfasern eines Rechteckquerschnittes. 


Nach A. Nadai’ kann die Formanderungslinie aus einem Biegeversuch berechnet 
werden. Aus den gemessenen Werten M, ¢,, «, ist nach Gl. (6) der Winkel # und 


damit aus der Linie ¢« = /(#) auch S bekannt. Aus Gl. (11) ergibt sich das Ver- 


haltnis a und damit aus Gl. (12) die den Randdehnungen entsprechenden Rand- 


spannungen. Ahnliche Beziehungen hat H. Herbert abgeleitet und aus dem Biege- 
versuch mit guBeisernen Stiaiben die Formanderungslinien des Biegezug- und Biege- 
druckrandes berechnet!°, deren Ubereinstimmung mit den im einachsigen Zug- und 
Druckversuch ermittelten als befriedigend bezeichnet wird. Fir Stahl mit ausge- 
pragter Flie8grenze wurden ahnliche Untersuchungen von E. Chwalla™ durchgefihrt. 
Hier soll noch ein Rechteckbalken aus Beton untersucht werden. In den Versuchen 
von Bach-Graf 190712 wurde der unbewehrte Balken Nr. 66 durch zwei Einzel- 
lasten auf reine Biegung in der bekannten Art stufenweise beansprucht. In den 
einzelnen Laststufen wurden die Formainderungen am Zug- und Druckrand auf einer 
MeBstrecke von 70 cm gemessen. In Tab. 1 sind aus Abb. 3 nach der oben erléuterten 
Methode einzelne Stufen berechnet und die Formanderungslinien fiir den Zug- und 
Druckrand in Abb. 4a aufgetragen. 


10 H. Herbert: Uber den, Zusammenhang der Biegungselastizitat des GuBeisens mit seiner 
Zug- und Druckelastizitat. Mitt. Forschg. Ing.-Wes., H. 89. Berlin: Julius Springer. 1910. 

11 BE. Chwalla: Uber die Erhohung der FlieBgrenze in prismatischen Balken aus Baustahl. 
Stahlbau 6, 149 (1933) und Zuschrift an die Schriftleitung. Stahlbau 1, 8 (1936). 

12 GC. Bach: Versuche mit Eisenbetonbalken, 2. Teil. Mitt. Forschg. Ing.-Wes., H. 45—47. 
Berlin: Julius Springer. 1907. 
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Yooh M10 “kg cm 
“ i) My, = 52.94.10 7%kg ei 
pale ke is it 


=— 


O00 


Q05- 


G3 @. 103 


Abb. 3. Randdehnungen und Biegemoment eines Betonbalkens in Abhangigkeit vom Neigungs- 
winkel ® zweier Querschnitte. 
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Abb. 4. Spannungs-Dehnungslinien an den Randern eines Biegebalkens aus Beton. 


Tabelle 1. Betonrechteckquerschnitt: 
d = 31°00 cm, b = 15°07 cm, 1 = 69°9 em, M, = 5340 kgem, bd? = 14480 cm’. 


M "loo "low 8103 | dey déy | oor aM oy Oy 


kg cm e138 6,38 dd dd oy dd ¢ kg/cm? kg/em2 


| 
Higengew. | 0°007 | 0007 | — | — ete re ae 
24040 | 0°034 | 0°035 | 0°077 | 0°471 | 0430 , 0912 | 18°60 8°78 9°62 
42800 | 0:077 | 0°067 | 07162 | 0°528 | 0°358 | 0-680 | 28:00 | 13:18 | 19:40 
493404 | 0:098 | 0°083 | 0°204 | 0°549 | 0°338 | 0°618 | 30°00 | 14:35 | 23°30 
55340 | 0°127 | 0:099 | 0°254 | 0°610 | 0294 | 0-483 | 23°50 | 14:10 | 28°60 


8 Die Dehnung infolge des Eigengewichtes wurde auf elastizititstheoretischer Grundlage 
berechnet. 


14 Wasserflecke an der Balkenunterseite. 
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Die dort eingetragene Spannungs-Dehnungslinie des Zugversuches stimmt mit 
der aus dem Biegeversuch berechneten nahezu iiberein. Die Spannungs-Stauchungs- 
linie des Druckversuches weicht bis zu 8°6°, von der aus dem Biegeversuch berechneten 
ab. Wird auf Grund der im Zug- und Druckversuch ermittelten Formanderungslinie?® 
nach Abb. 4b und Gl. (15) das Bruchmoment berechnet, so ist” 

14480 = 13°819°8 

My marry 13°8 + 19°8 
gegeniiber dem gemessenen Bruchmoment von 57:940 kgcm (Fehler 1°5%). Dazu 
ist folgendes zu bemerken: Nach den Beobachtungen von Bach-Graf, Turneaus® 
und anderen Verfassern treten nach dem Erreichen der Bruchdehnung bzw. Zug- 
festigkeit des Betons am Zugrand (Ri®- 
bildung) an der Balkenunterseite Was- Vie — 
serflecke auf. Die hier gemessene Bruch- sh 


= 58°800kg em 


dehnung nach dem Auftreten der ersten lad 
Wasserflecke betragt 0°09°% gegen oe 
0°08°%, beim Zugversuch. Die zugehérige 
Betonstauchung an der Balkenoberseite — yyl_yy 


wurde mit 0°08°% gemessen, gegen 
0°073°% fiir die nach Abb. 3b aus dem 
Druckversuch entnommene Spannung 
oq = 19°8 kg/em?. Diese Vergleiche 
zeigen, daB das Bruchmoment 
auf Grund der im. Zug- und 
Druckversuch ermittelten Form- 
anderungslinien sehr zutreffend 
berechnet wurde. Dazusind keinerlei 
Messungen des Biegeversuches not- 
wendig. Dagegen zeigen die mit Hilfe 


100 200 Ve 
der gemessenen Dehnungen am Zug- 3 
und Druckrand berechneten Formande- pb. 5. Spannungs-Stauchungslinie des Druck- 
rungslinien Abweichungen von den- randes eines Stahlbetonbalkens. 


jenigen des Zug- und Druckversuches. 

Eine Erklirung dafiir ware, neben der bekannten Unsicherheit von Dehnungs- 
messungen, die Anordnung der Mefstrecken, die bei einer Balkenhéhe von 30 cm 
nur 15cm von der Einzellast entfernt beginnen, wodurch insbesondere die Druck- 
meBstrecke kaum gleichmaBig beansprucht ist. , 

Der Bruchvorgang ist folgender: Nach dem ersten RiB des Betons an der Balken- 
unterseite schreitet die RiBbildung ins Innere des Balkens bis zum Bruch fort. Dabei 
steigt das Biegemoment vom Rimoment auf das Bruchmoment um 17°4% an. Dies 
steht in Ubereinstimmung mit der in der Literatur erwihnten Feststellung von FlieB- 
linien vor dem Erreichen des Flie8momentes in auf Biegung beanspruchten Stahl- 
tragern. Dazu sei ausdriicklich vermerkt, daB hier keinerlei Annahme tber 
die Spannungsverteilung im Bruchzustand vorausgesetzt ist. 

Es soli noch die Spannungs-Stauchungslinie des Druckrandes eines 
stark bewehrten Stahlbetonbalkens aus einem Biegeversuch ermittelt werden. 


15 Im Zugversuch wurde die Messung der Dehnungen etwa bis zur halben Bruchdehnung 
durchgefiihrt. Die Abhangigkeit der Bruchdehnung von der Zugfestigkeit des Betons wurde 
vom Verfasser aus Drehversuchen berechnet. Daraus ist der hier eingesetzte Wert ¢, = 0°085°/o, 
in Ubereinstimmung mit der Messung Bachs (s. Anm. 12, 8. 99) fur o, = 13°8 kg/cm? ent- 
nommen. Siehe K. Schaden: Die Ri®- und Bruchlast des auf reine Verdrehung beanspruchten 
Stahl- und Spannbetons. Osterr. Bauzeitschrift, Jg. 1953. 

20* 
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Von den Versuchen Bach-Graf 19106 sind aus der Reihe 2 fiir den Balken Nr. 328 die 
gemessenen Randdehnungen ¢,, €, in Tab. 2 eingetragen. Dort sind, wie oben nach 
Abb. 5a, die Randspannungen oj, 0, berechnet. Die Zugrandspannung ist jetzt infolge 
der Eiseneinlagen und RiSbildung im Beton eine fiktive Spannung, die hier zunachst 
nicht weiter untersucht wird. In der Abb. 5b ist die Beziehung o, = f (é) aut- 
getragen und mit der im Druckversuch ermittelten Forménderungslinie verglichen. 
Die Ubereinstimmung ist ausreichend!’, wobei zu bemerken ist, da8 wie beim un- 
bewehrten Betonbalken eine ungleichméBig beanspruchte Druckstrecke gemessen 
wurde (MeBstrecke 70 cm, Entfernung der Einzellasten 100 cm). Der Bruchvorgang 
ist hier folgender: Nach der RiB®bildung am Zugrand wird das Biegemoment von der 
Bewehrung und der Betondruckzone aufgenommen. Nach Erreichen der Prismen- 
festigkeit am Biegedruckrand (angenaihert Punkt 4 der Abb. 5) schreitet der Bruch 
des Betons ins Balkeninnere bis zur Erschédpfung der Tragfihigkeit des Balkens 
(Punkt 5) fort. 
Tabelle 2. Stahlbeton-Rechteckquerschnitt; 


d = 25'1 cm, b = 18°2cm,/ = 70cm, op = i. M. 105 kg/cm’, M, = 14400 kgem, bd? = 11460 cm’. 


M 10-3 | "loo "loo 108 dey des | cel | aM, 7 Gy 
kg/em &y Cy ado do | OG, dd | kg/cm? kg/cm? 
14418 | 0:029 | 0029 = a = Dales hl —s = 
114°4 0°132 0°213 0°482 0°273 OC: StOR Ee G2 ea OO an 74°6. 46°2 
214°4 0°278 0°478 1°058 0°246 0°458 | 1°86 150. | +1440 | 78°8 
314°4 0°465 0°885 1°885 | 0°206 0°498 | 2°42 170 242°0 | 98°0 
414°4 0°739 1°749 37490 || O73 7a nO75 725 |) 4220 132 419°0 | 104°0 
439°4 0°861 2°349 4°490 | 0°083 0°615 7°40 0 688°0 | 86°5 


Wie auch aus Abschn. E zu ersehen, ist der Ansatz der Prismenfestigkeit am 
Druckrand zutreffend. Eine ,,Stiitzwirkung“ ist im Rahmen der erreichbaren 
Versuchsgenauigkeit nicht von Bedeutung, da beim Bruch bereits in einem breiteren 
Randstreifen die Prismenfestigkeit erreicht ist und somit o, = op gilt. Das gleiche 
dirfte fiir Metalle gelten, wo auch die Festigkeit des Werkstoffes in der Biegeprobe 
verschieden von der in der Zug- bzw. Druckprobe ist (20°, Unterschied bei den 
Versuchen von E. Chwalla mit Baustahl''). Die Bruchstauchung am Druckrand 
steht in Ubereinstimmung mit den Prismenversuchen!® und betragt beim minder- 
festen Beton etwa 1°8°/),, beim hdher- und hochwertigen Beton etwa 2°f,). Ge- 
messene héhere Werte am Biegedruckrand schlieBen bereits die Verformungen des 
Druckrandes nach der Erschopfung seiner Tragfiihigkeit (Prismenfestigkeit) ein und 
haben daher nur scheinbare Bedeutung. 


IV. Der I-Querschnitt aus Verbundwerkstoff. 


Die Voraussetzung des Abschn. A, 2 ist nur dann berechtigt, wenn die Defor- 
mationen des Tragers gering sind, was bei spréden Stoffen, insbesondere Beton 
im allgemeinen zutrifft und auch fiir die im Bauwesen zugelassenen Durchbiegungen 

6 C. Bach und O. Graf: Versuche mit Eisenbetonbalken. Mitt. Forschg. Ing.-Wes., H. 90—92. 
Berlin: Julius Springer. 1910. 

17 Die Formanderungslinie o, = f (e,) wurde aus Biegeversuchen unter der in der Theorie 
des Stahlbetons tblichen Annahme einer gerissenen Betonzugzone und einer zur einachsigen 
Druckstauchungslinie affinen Druckspannungsverteilung im Querschnitt berechnet. Damit 
ergeben sich gegen die Bruchgrenze stark abfallende Forminderungslinien, die beim ebenen 
Spannungszustand durechaus unwahrscheinlich sind. Siehe z. B. H. Hamann: Berechnung der 
Druckspannungskurve im Stahlbetonbiegequerschnitt: Schweiz. Bau-Ztg. 44, 629 (1952). 

8 Eigengewicht, ¢ berechnet auf elastizititstheoretischer Grundlage mit Ez = 260: 10%. 

” K. Schaden: Die Verdrehungsfestigkeit des Betons. Osterr. Ingenieur-Arch. 6, 12 (1951). 
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angenommen werden kann. Bei hohen Biegungstrigern mit breiten Flanschen aus 
bildsamen Stoffen mit gréBerer Formanderung entsteht im Steg eine Radialkraft 
und zusiitzliche Formanderungen in den Flanschen2°, wodurch die nachfolgenden 
Formeln nur angenihert gelten. 

So wie im Abschn. B fiir den Rechteckquerschnitt, wird nachfolgend fiir den 
I-Querschnitt, der durch Einlagen im oberen und unteren Stegrand verstirkt ist 
(z. B. ein doppelt bewehrter Stahlbetonquerschnitt), nach Abb. 1 das Bruchmoment 
berechnet. Hier ist die Querschnittsbreite ) — © (y) nur streckenweise konstant. 
Sie springt nach Abb. 5 bei Ha Von by auf b, baw. b,. Die Integration von Gl. (8a) 
bzw. (9a) kann daher in drei Bereichen durchgefiihrt werden. Damit lautet die 
Bedingung Gl. (8a) 


ara te & 
\ br f(e) de + \ by f (e) de + \ by f (e)de ! Ba Te 4 Pal (eed 


55: 5 
1 dy “ds 


= 0. (17) 


Mit 


Say 0 ®de 


[jth (e)de = — Fy, \f(e)de = —F,”, \f(e)de= Fy, (/(eae = Fy’,| 


at — ea, 0 ed, +(18) 
Potts) _ _ Facalert+ es) _ pe Fefley) _ Facaler+*s) _ py | 
‘ be by do ; bo @ by do , 
yo a=—f(-a),  94=—f(—e4)  %,=—f(—&) nr 
Oo, =f (€), oa, =f (£a,)s oe, = f (€,) 
bezeichnen, erscheint die Bedingung XH = 0 nach Gl. (8) in der Form 
FPP +P +P," —F4+ Fe’ = 0. (20) 


Vor Behandlung der Bedingung Y M=0 Gl. (9) sollen die dazu erforder- 
lichen Formanderungsbeziehungen erledigt werden. Aus der Proportionalitat zwischen 
Dehnungen und Nullinienabstand ist nach Abb. 1 


fy ees &2— €q, a + q, bins &, + &, nt &g— &, 
d = d Sade h se h 2 
woraus . ; a tt a a 
| dy (€, + €) hy 
Eq, — €2 fe] 7e,) ad,” By a E9, 
0 0 (21) 
_— (&1 + &2)(dgo—4y) aoa nGae fa) he, | 
En, Fe dy Eg, Ee, wt ee do is ate 
Nach Gl. (21) und Gl. (6) wird 
20 20 20 
Eq, 1 2 = (dy — d,) ay a, Cie0 ar oh, hy —26,,, Cds Hoe dy i 
(22 
20 20 26 26 
Eq, — €2 = — 7 a Ey + &2 = 7p; Eg — &, = 7 ha, Ea + &, = 7 hy. 
Durch Differentiation nach di? wird aus Gl. (4), (21) und (22) 
dea =(dy—)>—de, de,=de,—Thy  dey,=de,— 7 dy, | 
2 2d 
de, => h,+de,, dey= 7 — dé, | (23) 


¢ 4 2 
ef — de, = dé, — To +2de,, de,= ve h, — dé. 


20 I. Kober: Stegbeanspruchung hoher Biegungstrager. Abhandlg. aus dem Aerodynamischen 
Institut a. d. T. H. Aachen, H. 7, 8. 34. Berlin: Julius’Springer. 1927. 
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Wird noch Gl. (21) nach @ differenziert und durch by dividiert, so ist mit Gl. (2) nach 
Umformung und Beriicksichtigung von Gl. (19) 


b 
ze 0, d&, = — $ (o4,déq, — 01 d€,) + o4,¢€a, — Fa, deg, + Bae Og, 4&4, + 
2F, 2F,, 
ms byt 4, Net O Gide, aeaurpre= by l * (6 Fe, tdog). (24) 
Die Summe der Momente um die neutrale Achse Gl. (9a) wird 
= OE “ds & , 
| f(e) eb, ede + | f(e) aby erde + | f(e) eb, @de + Fy, %, &,0 + Fe, 5, &,0 = M, (25) 
—& —€q, eq, 
oder mit Gl. (2) 
‘ —eé ‘L ( Ey ys oF, fT : 
M # = Fb \ f(e)ede + bo \ P(e) ede + b, | F(6) ede + —7~ 90, £5, TF Fe, Eeaf- (26) 
— & Sai, ed, 


Durch Differentiation nach dd in Gl. (26) wird mit Gl. (19) und Gl. (24) 


d ie 
ag (MB) = | {by 04, dea, (6a, + &2) — B10, 4ea, (60, + &2) + B10 dey (€, + &2) + 
2 


+ by 04, Léa, (€4,— Ey) — by oq, deg, (&4,— &3)'-+ ey ket 2 Oe, (& + é,) & (es 6,58 
2F, 2F, 2F, 2F, ) 
Pr Oleg mr On lbs | Sakae oe, eas ieee Fe,0 dee, (27) 
Unter alien ee der Gl. (22) und (23) wird 
# (M8) = vi 5. | dod 7 — de4| coy o4,|(do tag? > — de, 
2d, 20 2 206 
(dy a pa =| +8101 rt — dey) = dy — bo 04, (dex — 7 dy) =] ds + 
2 20 2) Pela 2F,. 292 2F,, 
HER See (de, = dali : ty 04, yy + 7 de, ~— yt "a, Ode, — 
8 28 - ae i 1 2F,, 29 
i etsy do,, h,+—7~ &,8 de, T %, 7 hel. (28) 


Nach tbe der Glieder mit # wird 
d 
ay (M8) = 7 {de * (oa, (do —d 1) (by — bo) — 6, 4, dy + 04,4 (by — by) —F,, 6, + Fe, 4¢,] + 
4 
zs ae o [o4, (dy — dy)? (by — 6,) +b, 0, di? + 4,49" (by — bg) + 2F,, Oe, hy — 2F',. onal}: (29) 


Mit Gl. (24) wird aus GI. (29), wenn 
A= Oa, (dy — dy) (by — 4), 


Er ==67 6, 0. 
C = a4, dy (by — 5), (30) 
D=F,,6,,, 
H=F, o,, 
bezeichnen und fiir den Bruch wie beim Rechteckquerschnitt oA ws = gilt, das 


Bruchmoment 
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Mp 
_ Lfj—(4?+ Bi + C? + D*— Et)—2[A(B+ 0 + D—E) + B(@' + D—E) + C(D—E)— DE] 
: eee A gem af 
ie dy dy— dy dy 


+A (d>—d,)+Bd,+Cd,4 2 (Ff, o.,h1 — F,, o,, hs)». (31) 


Fiir den Rechteckquerschnitt wird mit d, = d, = d,=d, 6b, =b,= 6b) = 8, 
Ff, =F; =0 in Gl. (30) 
A=0, .B=b60,d,.0=0,..D=0,.F=0, 


bd? a, a, 
2 o, +0, (32) 


wie in Gl. (15), wobei die Bedingung YH = 0 nach GI. (14) gilt. Fiir den T-Quer- 
schnitt wird mit d, = 0, b, = bo, F, =F,’ = 0 in GL. (30) 


Ot opSeby od MOS ody byes y, inp sig, 


Mess 


1 /— B?—C?—2BC 
M, == (——=— e+ Ba + Cay, (33) 
ial Sac 
0 2 


Die fiir die Berechnung von A bis C erforderlichen Werte von o werden durch Iteration 
aus Gl. (19), (20) und (21) mit der bekannten Bruchranddehnung und angenommenen 
Werten der zugehérigen Bruchrandstauchung oder umgekehrt erhalten, je nachdem 
der Bruch durch das Versagen des Zug- oder Druckrandes eintritt. Fiir einen Werk- 


stoff mit gleich groBer Zug- und Druckfestigkeit und c = 7a = 1 wurde “mit 


zx = ), wodurch Gl. (20) tiberbestimmt wire. 
0 


Daraus ist zu ersehen, daB auch bei einem solchen Werkstoff die gleichzeitige Erreichung 
der Festigkeit am oberen und unteren Rande nur bei zur Schwerlinie symmetrischen 
Querschnitten méglich ist, was mit dem Ergebnis der elastizitatstheoretischen 
Berechnung iibereinstimmt. 

Fir den I-Querschnitt gilt Gl. (31) ohne die Gleder F,, und F,. Die Be- 
dingung Gl. (20) wird so wie oben durch Iteration aus Gl. (21) fiir gegebene Werte 
von ¢, bzw. e, erfiillt. Fiir den in Hohe und Breite symmetrischen I-Querschnitt 
mit d, =d,=—d, b, = 6,=5 und c=1, das heiBt ¢, = «, = e, wird aus Gl. (21) 


Ep nares & =a (1— de ). (34) 
Die Gl. (18) ergibt — F,'’ = F,', — F,” =F", womit Gl. (20) erfiillt ist. In Gl. (30) 
wird mit o4,=04,= %a 

A =o,(d,—4@) (6, — 6), B=bo,d, C=0,d (by — 9), (35) 
womit sich das Bruchmoment nach Gl. (31) ergibt. Fiir lineare Spannungs-Dehnungs- 
beziehungen bis zum FlieBen bzw. Bruch und c = 1 wird mit ¢ = e- # in Gl. (35) 


( 2d 2d 
A=e,E(1-,) (d,—d)(b,— 6), B=«e,Hbd,, O=E(1-= | (b)>—). (35a) 
Fur d=0 oderd=d, b=), & H=o wird in Gl. (35a) 

AO. bes St Oe, a) 


und Gl. (31) ayes bd? 
Pome. 


é) = & = & aus GI. (21) a, = «& (1 _ 


2 


Eq = 


oO 


wie fiir den Rechteckquerschnitt nach Gl. (16). 
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Im zugbewehrten Rechteckquerschnitt wird nach Gl. (30) mit d=d,=d 
oder. d.=,0, 0) =) = ye 


A=0; JB=0iyd JC= 0". Delo eae, 


YR 
M,= 5 |b @—272_ + F,o,\2h -—__——] . (36) 


0, + oO 0, + G2 id 
Wird die Zugbeanspruchung im Grundwerkstoff nicht in Rechnung gestellt (Aus- 
schlu8 der Betonzugzone im Stahlbetonbalken), so wird Gl. (36) mit o, = 0 


M,=5F,0,(2h — 522"). (36a) 
Aus Gl. (20) folgt mit 

0 &y ) 

\f(e\de=—F,, \f(e)\de= —F,,-. a) | 
—& 0) (37) 

FF, ¢ (Ey + €o) - 

ad 2 — Fe, b) } 

Fe 2 UP ei ee (38) 


Die rechte Seite dieser Gleichung ist aus den bekannten Bruchdehnungen ¢, und e, 
gegeben. Aus Gl. (37b) ergibt sich daraus die beim Bruch in der Bewehrung vor- 
handene Spannung o,, wobei zu beachten ist, daB die nach Gl. (21) fiir einen an- 
genommenen Wert von ¢e, berechnete Dehnung der Bewehrung e, bereits aus der 
Beziehung o, = f (€,) die Spannung in der Bewehrung festlegt. Da F, > F,, kann 
F,—F, =F, gesetzt werden, womit 
ue bd : 

Bogue CaP FP? wobei o,=f(€,). (39) 
Durch Iteration ergibt sich daraus die beim Bruch des Stahlbetonbalkens in der 
Bewehrung vorhandene Spannung o,. 


Im zug- und druckbewehrten Rechteckquerschnitt mit dj =d oder 


d= 0, b) = 6b ist mit Gl. (30) d=0, B=bo,d, C=0, D=F,6,, E=F,,6,, 
LN Ghee) ‘ Fo 
$00; = id o, Fas 
forte ES ety e y | Foe (34 ai y } 
2 Ot Oy he eee Ga Uy de te ee 0 + Op 
(40) 
Bei Ausschlu8 der Betonzugspannungen wird mit o, = 0 aus Gl. (40) 
1 Ff, o, i Go, 
M, =F Fe, %, (22, — = pee.) see Pa (2 h,— apt (40a) 
Mit den Bezeichnungen nach Gl. (37) bzw. (18) wird Gl. (20) 
F,—F,=F* —F’, (41) 


woraus wie vorhin die Spannungen o,, o,, und o, fiir den bekannten Wert von o, 
gegeben sind. Wegen fF, > F, wird wie oben 


F, = arts (ie, Oe, — Be Oe,)s | 
I 
Cn ee (& + &) + seb, \ (42) 


Mit den bekannten Beziehungen o,, = f (¢,,), Oo, =f [&, (€,)] stellt Gl. (42) eine 
durch Iteration eindeutig lésbare Bestimmungsgleichung fiir o,, und o,, dar. 
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' In der gleichen Art bestimmt die Ql. (20) die Spannungen fiir Gl. (31), woraus 
sich noch die Bruchmomente fiir den T- und I-Stahlbetenquerschnitt ergeben. Fiir 
den T-Stahlbetonquerschnitt wird bei Ausschlu8 der Betonzugzone mit 
d,=0, b= b, o4,=0, h, =h in G1. (30) A=0, B=0, C= G4, da (bo — 52), 
D=F,o,, E=0 und in Gl. (31) 
—0(2D+0) 


C 
by Og + = 


M,= +Cd,+2Dh, (43) 


1 
2 


wobei o4, und o, nach Gl. (19), (20) und (21):durch Iteration fiir den bekannten 
Wert «, und den angenommenen Wert von ¢, zu bestimmen ist. 


V. Versuche mit Stahlbetonbalken. 


Fiir einen unbewehrten ‘Betonbalken mit rechteckigem Querschnitt wurde das 
Bruchmoment bereits in Abschn. © berechnet. Weitere Versuche mit unbewehrten 
Betonbalken und ausreichender Angabe der Betonfestigkeiten im einachsigen Versuch 
sind nicht allgemein bekannt. Von den zahl- 
reichen Balkenversuchen mit einfach bewehrtem 
Rechteck- und Rippenquerschnitt sowie doppelt 


Op-ta 


100 


700 


4 2 2 Eco G2 OY G6 4 kgfem? 


Abb. 6. Flacheninhalt der Spannungs-Stauchungslinien von Beton. 


bewehrtem Rechteckquerschnitt wurden diejenigen ausgewertet, welche ausreichende 
Versuchsprotokolle ausweisen und wo die Versuchsdurchfiihrung keine Mangel aufzeigte, 
soweit die Literatur dem Verfasser erreichbar war. Versuche mit Stahlbetonbalken von 
doppelt bewehrtem Rippenquerschnitt sowie I-Querschnitt sind bisher nicht allgemein 
bekanntgeworden. Die Druck-Stauchungslinien des Betons sowie die Werte von F, 
sind fiir mehrere GréBen der Prismenfestigkeit in Abb. 6 aufgetragen. Da eine Angabe 
der Zugfestigkeit des Betons fast ausnahmslos fehlt, wurde diese mit ein Zehntel der 
Prismenfestigkeit angenommen. Die Versuche wurden durchwegs in der bekannten 
Art ausgefiihrt, daB stufenweise belastet und wieder entlastet wurde. Zum Teil er- 
folgte kurz vor dem Bruch wieder eine Entlastung (Bach-Graf, Gehler), zum Teil 
wurde ab einer weit niedrigeren Laststufe bis zum Bruch ohne Entlastung belastet 
(Empa-Versuche). Dadurch ist bei der gemessenen Bruchlast die Mitwirkung der 
Betonzugzone mehr oder weniger ausgeschaltet. Dies wurde bei der Berechnung 
der Bruchlast dadurch beriicksichtigt, da aus dem Unterschied der theoretischen 
Werte des Bruchmomentes mit und ohne Betonzugzone der Anteil der Betonzugzone 
am Bruchmoment berechnet wurde. Dieser Anteil ist angenaihert proportional dem 
Verhaltnis zwischen der letzten Laststufe vor der Bruchlast und Bruchlast. Durch 
Addition des tatsichlichen Anteiles der Betonzugzone zum Bruchmoment ohne 
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Betonzugzone ergibt sich der theoretische Wert des Bruchmomentes eines Versuchs- 
balkens. Damit konnte auch das Bruchmoment schwach bewehrter Balken zutreffend 
berechnet werden. Nach R. Saliger®! entstehen durch das Schwinden des Betons 


E, ; 
im einseitig bewehrten Rechteckbalken bei ¢, = 10: 10° n = ae 15d, = Ole 


je nach dem Bewehrungsgrad ps = a von 4 = 0°005 bis 0°040, am Betonrand Zug- 
1 


spannungen von .2°9 bis 10°4 kg/em?, Druckspannungen von 1:2 bis 4°3 kg/cm? und 
in der Bewehrung Druckspannungen von 170 bis 25kg/em?. Der Einflu8 dieser 
Spannungen auf das Bruchmoment wurde zunachst in der Rechnung nicht bertick- 
sichtigt, da dieser wegen der nur angeniherten Angaben iiber die Baustoffgtiten 
(Prismenfestigkeit nicht aus der Balkendruckzone bestimmt, Zugfestigkeit geschatzt, 
teilweise Ausschaltung der Betonzugzone usw.) innerhalb des hier erreichbaren 
Streuungsgebietes liegt. 

In den Versuchen von Bach-Graf 1907” ist fiir den einfach bewehrten Recht- 
eckbalken Nr. 52 mit den Abmessungen d = 30°54cm, h = 29:04cm, 6 = 15°09 cm 
und den Werkstoffkennwerten op= 146 kg/em?, oz = 13°6 kg/cm?, F, = 2°36 cm’, 
o, = 3200 kg/cm? aus Gl. (39) o, = 6, = 3200 kg/cm?, womit bei Mitwirkung der 
Betonzugzone nach Gl. (36) das Bruchmoment M,’ = 274000 kgem wird. Bei Aus- 
schluB der Betonzugzone wird das Bruchmoment nach Gl. (36a) M,’’ = 205000 kgem. 


Dabei ist das Verhaltnis der letzten Laststufe vor dem Bruch zur Bruchlast 
8000 


3500 — 0°94. Der Anteil der Betonzugzone am Bruchmoment ist 274000 — 205000 = 
= 69000 kgem, wovon nur noch (100 — 94)% = 6% wirksam sind. Damit ist das 
Bruchmoment M, = 205000 + 0:06 - 69000 = 209150 kgem gegentiber dem_ ge- 
messenen Bruchmoment von 218800 kgcm. 

In den Versuchen von M. Ro8* ist fiir den mit Torstahl 40 bewehrten Platten- 
balken TS 25014 mit den Abmessungen 6, = 17cm, b,= 43cm, d, = 32cm, 
d,= 12cm, h,= 29cm und den Werkstoffkennwerten op= 216kg/em?, o, = 
= 21:6 kg/em?, F&F, = 4:'74cm? und nach Gl. (20) o, = 4700 kg/em?, womit nach 
Gl. (30) A = 0, B = 11800, C = — 6750, D = 22300 und bei Mitwirkung der Beton- 
zugzone nach Gl. (31) das Bruchmoment M,’ = 760000 kgem und bei Ausschaltung 


der Betonzugzone das Bruchmoment M,’’ = 590500 kgem wird. Das Verhiltnis 


; 3000 : A 
der Laststufen ist “350 = 0°34, womit wie oben das Bruchmoment M, = 590500 + 


+ 112000 = 702500 kgem wird, gegentiber dem gemessenen von 708000 kgem. 

In den Versuchen von W. Gehler, H. Amos und E. Friedrich” ist fiir den 
doppelt bewehrten Rechteckbalken Nr. 103 mit den Abmessungen b = 30 cm, 
d = 205 em, h, = 17cm, h, = 3°5cm und den Werkstoffkennwerten op = 123 kg/em?, 
oz = 123kg/em’?, F, =F, = 628cm*, o, = 2425kg/em? nach Gl. (41) o, = 
= 2815 kg/cm’, o,, = 590 kg/em?, womit das Bruchmoment bei Mitwirkung der 
Betonzugzone nach Gl. (40) M,’ = 296500 kgem und bei Ausschlu8 der Betonzug- 
zone nach Gl. (40a) das Bruchmoment M,’’ = 244850 kgem wird. Mit dem Last- 
stufenverhaltnis ae = 0°92 wird wie oben das Bruchmoment M, = 244850 + 
+ 4100 = 248950 kgem gegeniiber dem gemessenen von 245200 kgem. Von diesen 
Versuchen wurde noch der Balken Nr. 114 ausgewertet, der den theoretischen Wert 


“1K. Saliger: Der Stahlbetonbau, 7. Aufl. Wien: Franz Deuticke. 1949. 

~ M. Ros: Festigkeit und Verformung von auf Biegung beanspruchten Hisenbetonbalken. 
EMPA-Bericht Nr. 141. Zurich. 1942. 

*8 W. Gehler, H. Amos und EK. Friedrich: Versuche an Stahlbetonbalken zur Bestimmung 


der Bewehrungsgrenze. Deutscher AusschuB fiir Stahlbeton, H. 100. Berlin: W. Ernst u. Sohn. 
1943. 
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fiir das Bruchmoment bei Ausschlu8 der Betonzugzone von M, = 176300 kgcm 
ergibt. Der gemessene Wert von 162430 kgem liegt wesentlich niedriger, was mit 
der tatsaichlich erfolgten Zerstorung des Balkens an den Auflagern und nicht durch 
reine Biegung in Ubereinstimmung steht. In der gleichen Art wurden die anderen 
auswertbaren und kennzeichnenden Balken der soeben genannten Versuche und 
weitere der Versuche von Bach-QGraf 191016, A. Brandtzaeg™ und der Schmidt- 
stahlwerke mit Torstahl 6025 ausgewertet. Die theoretischen Werte sind in Tab. 3 
mit den gemessenen verglichen und zeigen eine Abweichung von héchstens 9°%. 


Tabelle 3. 
. MRect gp. |" ¢@ 
Vv h | : echnung P Berechnet 14, 
ersuch von Probe Nr. Querschnitt op RE kelon® | ke hae Le moe ar Stahl 
1 | | 
Bach-Graf 66 |Rechteck | 1015 | 146 | Heh U8 
1907 52 » | 0960 | 146 | 3200 | 0:0027 | 36, 36a | — 
Bach-Graf 328 | Rechteck | 0°953 104 3400 | 0°075 36, 36a 
1910 395 Rippen | 1:090 | 120 | 3400 | 0°018 31, 43 — 
se) se ts Re | es oa i 2 S| ees Prete cat me EWES 
EMPA TS 300 Rechteck 0°990 | 244 | 4140 | 0°0045 36, 36a | Torstahl 40 
1942 | O°972 | 204 | 3840 | 0°0089 
TS 250 - | O978 | 244 | 3970 | 0°0183 
— 0913 157 | 4140 | 0°0045 
| 0°930 | 160 | 3840 | 0°0089 
S | 07965 | 209 | 3970 | 0°0183 
ST 52/250 | | 0°968 150 | 4520 | 0°0087 | St 52 
St N/300 — ac 1080 | 279 | 3250 | 0°0142 St N 
TS 25014 | Rippen | 0°993 216 | 3955 | 0:0038 31, 43 | Torstahl 40 
Gehler 103 | Rechteck 1:020 | 123 2425 | 0°0123 = yp’ | 40, 40a 
1943 75 | se | 0°923 120 | 2925 | 0:0017 36, 36a -_- 
88 6 07933 | 118 | 2425 | 0°0123 
99 os _ 07998 117 | 2425 | 0°0308 
= | <— = = == ae 
Brandtzaeg Ab Rechteck | 1°025 246 | 3010 | 0°0499 | 36, 36a 
1935 Ca - | 1020 | 193 | 2990 | 0:0502 —- 
Ea a | 0°962 | 126 | 2930 | 0°0512 
Schmidt- A, _Rechteck | 0°972 | 175 | 6010 | 0°0029 36, 36a | Torstahl 60 
Stahlwerke B, _ / 0°998 |} 190 | 7390 | 0°:0077 | 
1950 as Rippen 1°015 218 | 7040 | 0°0047 31, 43 
1B | 35 | 0°928 175 | 6990  0°0117 


(Hingegangen am 23. Juni 1953.) 


Verdrehung bildsamer Metallstaibe iiber die FlieSgrenze. 
Von W. Starnberg, Wien. 
Mit 4 Textabbildungen. 


Zusammenfassung. Die Verdrehungsbeanspruchung von Metallstiben an der FlieB- und 
Bruchgrenze einschlieBlich des Verfestigungsbereiches wird mit Hilfe eines in der Grenzbedingung 
vorkommenden Parameters berechnet. Die Sonderfalle des FlieBens und des Sprédbruches werden 
mit Hilfe der speziellen Werte dieses Parameters erhalten. Die berechneten Werte zeigten gute 


2% A, Brandtzaeg: Hisenbetonfragen. Festschrift der Techn. Hochschule Drontheim, 1935; 


Eisenbetonfragen. Beton u. Eisen, H.,14 (1936). 
2% Magistrat der Stadt Wien: M.-A. 39, Stadt. Priif- und Versuchsanstalt, Zah! MA 39, 


M 318/49. 
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Ubereinstimmung mit den Versuchsergebnissen. Die Gleichungen fiir die Normal- und Schub- 
spannung und fiir den Gleitwinkel 2 « wurden aufgestellt. 


Summary. By means of a parameter occurring in the limiting condition, the torsional strain 
of metallic bars at the yielding and breaking points is computed. The special cases of flow and 
of brittleness fracture are obtained from special values of that parameter. The computed values 
agree with experiments. The equations for the normal and shearing stresses and for the slip 
angle 2 « are also established. 


Résumé. Au moyen d’un paramétre qui se trouve dans la condition aux limites, l’auteur 
détermine l’effort de torsion des barres métalliques aux limites d’écoulement et de rupture. Les 
cas particuliers de l’6coulement et de la rupture cassante résultent des valeurs spéciales de ce 
paramétre. Les valeurs obtenues sont de bon accord avec les essais. Enfin, l’auteur déduit les 
équations pour les efforts normal et de cisaillement et pour l’angle de glissement 2 a. 


I. Einleitung. 


Die Grenzen der Beanspruchung werden je nach der FormanderungsgréBe ge- 
schieden in: FlieB-, Verfestigungs- und Bruchgrenze. Die Spannungen grenzbean- 
spruchter Korper fiir je eine der vorerwihnten Grenzen wurden bisher immer fiir 

sich gesondert ermittelt. So z. B. untersucht Nadait 

(S. 117 und 121) die Biegebeanspruchung eines Stahl- 
% stabes getrennt fiir den ideal plastischen Fall und fiir 

den Fall mit Verfestigung. Aufgabe dieser Arbeit soll 

es sein, die Spannungen bei der Verdrehung von 

Metallstaiben fiir alle oberwahnten Grenzen zu 

berechnen, so daB sie alle mit Hilfe eines Parameters 

analytisch erfaBt werden. Die Fille, welche etwa den 

Arbeiten von Nadai!.? und Schaden® tiber Verdrehung 

zugrunde liegen, miissen sich dann als Sonderfalle er- 

geben. Nadai hat namlich eingehend die Verdrehung von 


C= Ci/fy 27 


0 ef Lapa e Staben bildsamer Metalle beim FlieBbeginn sowie bei iiber 
Abb. 1. Abschitzung der Zahl en Querschnitt voll ausgebildetem plastischem Zustand 
¢ = 64/9; untersucht. Schaden hingegen hat die Verdrehungs- 


festigkeit von Staben spréder Stoffe (Gesteine, Beton, 
GuBeisen usw.) berechnet. Um nun alle Stadien vom FlieBbeginn bis zum Bruch zu 
erfassen, wird hier die von Leon?‘ eingefiihrte Kennzahl 


C= Op/oz (1) 
als Parameter verwendet. Leon hat zwar durch diese Kennzahl nur das Verhaltnis 
zwischen Druck- und Zugfestigkeit festgelegt, damit aber gleichzeitig eine wichtige 
GroBe zur quantitativen Erfassung des Sprédigkeitsgrades von Stoffen geliefert. 
So werden z. B. mit c= 1 bildsame Metalle (Baustahl, Aluminium usw.), mit 
c= 3 +5 GuBeisen, mit c = 10 + 25 Beton und Gesteine usw. erfaBt (s. Hiitte 15, 
8. 765). Torre® hat die Kennzahl ¢ Gl. (1) dahingehend verallgemeinert, da8 sie 
bei ein- und demselben Stoff veranderlich angenommen wird. Diese Annahme folgt 
aus der Beobachtung der Spannungs-Dehnungslinien (Abb. 1) fiir Druck und Zug 
eines bestimmten Stoffes. Fallen diese zwei Linien zusammen, dann gilt fiir den 


1 A. Nadai: Der bildsame Zustand der Werkstoffe. Berlin: J. Springer. 1927. 

2 A. Nadai: Plastizitat und Erddruck. Handbuch der Physik, Bd. VI, 8.428. Berlin: 
J. Springer. 1928. Plastic Torsion, Transf. of the Americ. Soc. of mech. Ingineers, Appl. 
Mechanics, 8. 29. 1930. Theory of Flow and Fracture of Solids. New York: Mc Graw-Hill Book 
Comp. 1950. 

8 K. Schaden: Osterr. Ingenieur-Arch. 6, 12 (1951). 

4 A. Leon: Ingenieur-Arch. 4, 421 (19338). 

5 Hutte I: Ing.-Taschenb., 27. Aufl. Berlin: Wilh. Ernst & Sohn. 1941. 

6 C, Torre: Osterr. Ingenieur-Arch. 4, 174 (1950). 
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ganzen Bereich etwa von der Proportionalitaétsgrenze bis zur Bruchgrenze c = 1. 
Trennen sich aber diese zwei Linien etwa nach der Proportionalitaétsgrenze, dann ist 
die Zahl ¢ von dieser Grenze an bis zum Bruch verinderlich. Insbesondere beweegt 
sie sich bei bildsamen Metallen, welche hier ausschlieBlich betrachtet werden sollen, 
in den Grenzen 

ls ors oi bis Ss. (2) 


Die Versuche von Ro& und Hichinger’ zeigen dies sehr deutlich, daB die Spannungs- 
Dehnungslinien. auseinandergehen. 


Torre® bezeichnet zum Unterschied von Leon die Zahl ¢ (s. auch Abb. 1) mit 
© == 0 g/Oes (3) 

Fiir c = 1 befindet man sich etwa an der FlieBgrenze des betreffenden Stoffes, fiir 
¢ > 1 an einer signierten Formiinderungsgrenze «* nach Abb. 1. An der FlieBgrenze 
wird o, = o, (Quetschgrenze = FlieBgrenze bei Druckbeanspruchung) und o, = o, 
(Streckgrenze = FlieBgrenze bei Zugbeanspruchung), an der Bruchgrenze hingegen 
wird og = op (Druckfestigkeit) und o, = o, (Zugfestigkeit). Sohin sind also GI. (2) 
und (3) fiir ein- und dasselbe Material zu nehmen und nicht wie Leon 
gemeint hat, da’ mit jeder Zahl c¢ jeweils ein Material charakterisiert wird. 

Hiermit ist die Méglichkeit gegeben, die Beanspruchung der Metalle ein- 
heitlich, das heiBt mit einer einzigen Grenzbedingung bei Anderung 
des Parameters c, von der FlieBgrenze bis zum Bruch zu erfassen. 
In dieser Erfassung der verschiedenen Beanspruchungsgrenzen liegt der eigentliche 
Zweck dieser Arbeit, wobei die Verdrehung bildsamer Metallstibe nur als Beispiel 
dient, um die erwahnte Moglichkeit an einem einfachen Fall zu veranschaulichen. 
Damit wird erwiesen, da eine einheitliche Theorie der Grenzbeanspruchungen, die 
zunachst von Torre® untersucht wurde, durchaus méglich und auch zweckdienlich ist. 

Die in diesem Beitrag vorkommenden Spannungsgr6Ben werden immer auf den 
jeweiligen Querschnitt bezogen. Das bedeutet nach Ludwik® z. B. fiir den ein- 
achsigen Zug bei verschwindender Volumsinderung e = 0: 


Oz = Oz (1 + €z), (4) 


wobei e, die Bruchdehnung ist. Ist die Volumsanderung verschieden von Null 
(e + 0), dann folgt der gleiche Ausdruck nach Torre® (S. 188, Anm. 27) 


EX -&7 
07 = 1 a e f O70: (5) 


Das Material wird nicht als isotrop und homogen vorausgesetzt. Die anisotropischen 
Eigenschaften eines Materials kommen mathematisch durch Spannungs-Dehnungs- 
beziehungen zum Ausdruck. Da hier ausschlieBlich mit den Gleichgewichtsbedingungen 
und der Grenzbedingung gearbeitet wird, erstreckt sich die Giiltigkeit der so ge- 
wonnenen Gleichungen auch auf anisotrope Stoffe. Nach Mises'® und Hill’ unter- 
scheidet sich wohl die Grenzbedingung fiir isotrope und anisotrope Stoffe durch die 
in ihr auftretenden Koeffizienten. Da es sich hier aber um ein quasiisotropes Metall 
handelt, kann man sich damit begniigen, die Anisotropie durch die Zahl c zu erfassen. 


tia ie Rok und A. Eichinger: Bericht Nr. 14 der EMPA in Zirich, 1926. Bericht Nr. 34, 
Zurich, 1929. Bericht Nr. 172, Zurich, 1949. \ 

8 ©. Torre: Schweiz. Arch. angew. Wiss. Techn. 15, 116, 145 (1949); Osterr. Ingenieur- 
Arch. 4, 93 (1950). 

9 P. Ludwik: Elemente der Technologischen Mechanik, S. 18, Anm. 15. Berlin: J. Springer. 
1909. 

10 R, v. Mises: Z. angew. Math. Mechan. 8, 161 (1928). 

11 R. Hill: The Mathematical Theory of Plasticity. Oxford. 1950. 
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Das wiirde sich hierdurch auswirken, daf ein Metall verschiedene Festigkeiten o, und a, 
in verschiedenen oder mindestens in einigen Richtungen hat. Diese Kigenschaften 
treten meistens bei Gesteinen deutlich zutage, wahrend bei bildsamen Metallen fast 
keine diesbeztiglichen Unterschiede zu erwarten sind. 

Die Anregung zu dieser Arbeit verdanke ich Herrn Dozenten Dipl.-Ing. Dr. techn. 
C. Torre. 

II. Die Grenzbedingung. 

Die von Torre® vorgeschlagene und eingehend untersuchte Grenzbedingung fiir 
spréden Bruch und plastisches Verhalten bildsamer Metalle wird gegenstandlicher 
Arbeit zugrunde gelegt. Die Grenzbedingung wird nachfolgend im Sinne ihrer An- 
wendung noch eine Erweiterung erfahren. Sie diente bisher nur allgemeinen 
Festigkeitsuntersuchungen, wihrend sie hier zur Berechnung der Spannungen in 
einem tordierten Stab herangezogen wird. Sie lautet zunachst in Hauptnormal- 
spannungen ausgedriickt 


(a, — 04)? + (6; — 3)? + (63 — 0,)? + 2 (ce — 1) a, (0; + 62 + 53) = 2¢6/7, | 
MiG tae AM 81S ae a 
bzw. in Koordinaten des Spannungstensors 
(6, — 04) + (6y — 04)? + (6, — 02)? +6 (tay? + tye +2) + 
4+2(¢—1)a,(t, +0, +,) =2¢0,. (7) 
Die Grenzbedingung Gl. (6) stellt im Koordinatensystem der Hauptnormalspannungen 


ein Paraboloid dar, das fiir c > 0 (was hier tibrigens der Fall sein wird) zur negativen 
Kompressionsachse 


(6) 


01 = 0g = 6, <0 
offen ist. Fiir c = 1 geht die Grenzbedingung Gl. (6) in die bekannte FlieBbedingung 
von Mises!2 und Hencky? tiber (sog. Gestaltténderungstheorie). 


Fiir den Fall der reinen Verdrehungsbeanspru- 
chung (Verdrehung unter itblichen atmosphirischen 
Verhaltnissen und Zimmertemperatur) gilt 


On 0" on oe 0. 


‘4 


Dann lautet Gl. (7) fiir diesen Fall 


c 
z Tye +, a a 3° r Op" (8) 


In Abb. 2 ist ein Querschnitt des verdrehten Stabes 
gezeichnet mit der Stabachse parallel z-Achse. Die 
Querschnittsebene liegt in der Koordinatenebene (x, y). 
Man sieht aus Abb. 2, da die resultierende Schub- 
spannung 


Abb, 2. Spannungen im Stab- 
querschnitt. 


Trop == Teo iTe yes (9) 
Die Unterschiede t,, und t,, usw. sind hier bedeutungslos, da bekanntlich t, , = tT», 
ist. Aus Gl. (8) und (9) ergibt sich, daB bei der Grenzbeanspruchung die resultierende 
Schubspannung ,,, in jedem Punkt des Querschnittes Abb. 2 der Grenzbeanspruchung 


Bey / gleich sein mu - 
Treg = + OY ie (10) 


12 R.v. Mises: Nachr. Ges. Wiss. Gottingen, H. 4, 582 (19138). . 
18 H. Hencky: Z. angew. Math. Mechan. 4, 323 (1924). 
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Fiir eine Verdrehung, bei der die mittlere Hauptnormalspannung o, = 0 ist, gilt 
bekanntlich 
d= 0. 6; = —.6, =t,. (11) 


Dann wird nach Gl. (6) und (11) 


Cc 
= 46,)/%. (12) 
Aus Gl. (11) und (12) ist ersichtlich, daB die Grenzschubspannung 1, gleich der 
resultierenden Schubspannung 1,,, eines grenzbeanspruchten Querschnittes ist. Es 
ist wichtig, zu Gl. (12) hervorzuheben, daB die Grenzschubspannung von den Festig- 


keitskonstanten abhingt. Dies sind die Zuggrenzbeanspruchung o, und die Zahl c 
nach Gl. (3). 


Nadai’ rechnet mit Gl. (9) als Grenzbedingung, indem er fiir ideal plastische Stoffe 
T.. = & = konst. 
annimmt. Es bleibt offen, ob fiir 


k= oder k=—& 

3 oder 3 
genommen werden soll. Der zweite Ausdruck stimmt mit Gl. (10) fiir ¢ = 1 iiberein. 
Somit gilt nach Nadai eet RRS. (13) 


Die Grenzbedingungen Gl. (8) und (13) besagen, daB es sich hier um gleiche Grenz- 
bedingungen, aber mit verschiedenen Konstanten handelt, so daB die Berechnung 
der Verdrehungsfestigkeit in Anlehnung an die Ergebnisse von Nadai erfolgt. Es 
ist aber, wie zuvor schon erwahnt wurde, die erweiterte Bedeutung der Ergebnisse 
hervorzuheben. Insbesondere, da man auf diese einfache Art in Gebiete iiber die 
Plastizitatsgrenze kommt und da} der ganze Verformungsbereich von der Plastizitats- 
grenze bis zum Bruch mit Hilfe des Parameters c analytisch erfaBt wird. 


III. Die plastische Spannungsfunktion. 


Die im vorigen Abschnitt genannten Schubspannungen Tt,, und 7,, miissen auBer 
der Grenzbedingung Gl. (8) noch der Gleichgewichtsbedingung in der z-Richtung 
Cie 


é 

y 
geniigen. Die anderen zwei Gleichgewichtsbedingungen in w- und y-Richtung er- 
geben nur, wegen o, = 0, = 0, = Tz, = 0, daB sich t,, und t,, in der Richtung der 
Stabachse nicht andern diirfen, also ér,,/¢, = 0 und 0t,,/0, = 0. 

Es wird hier nach Nadai! (S. 95) die Spannungsfunktion F (a, y) eingefihrt, 
die die Gleichgewichtsbedingung Gl. (14) befriedigt, was sich aus 
oF oF 
“ey” Cy eee ee (15) 
ergibt. Sie kann geometrisch als eine tiber den Stabquerschnitt gespannte Flache 
angesehen werden. Setzt man Gl. (15) in den Ausdruck fiir die Grenzbedingung 
G1. (8) ein, dann mu8 die Spannungsfunktion F der Gleichung 

oF \2 OF \2 ¢ 
(Ze . (Fa) mcs eo 

geniigen. Es ist wichtig, Gl. (16) noch in Vektorform zu schreiben 


OT mz 
Ox 


= (14) 


T xz = 


(VF): (VF) =(VFP =z 02 
bzw. ‘ee 
| grad F | =o.) = (17) 
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Gl. (17) besagt, daB der Gradient von F, das heift die groBte Steigung der Flache F, 
konstant ist, waihrend bei der Beanspruchung im Gebiete der elastischen Form- 
anderungen diese Steigung jeweils gleich der resultierenden Schubspannung 1,,; an 
der betreffenden Stelle des Stabquerschnittes ist (Erster Satz von L. Prandtl). 

Es miissen noch die Randbedingungen lings des Stabquerschnittes 

d Cys oF / éF 
(i) “rte gana (18) 

erfiillt werden. Diese Gleichung sagt aus, da8 die resultierende Schubspannung Tres 
in jedem Randpunkt des Querschnittes tangential zur Randkurve y (x) gerichtet ist. 

Der zweite Satz von L. Prandtl besagt, da das Drehmoment M des Stabes gleich 
dem doppelten Inhalt des von der Spannungsfliche F (x, y) eingeschlossenen Raumes 
ist. Nach Abb. 2 und Gl. (15) betraigt das Drehmoment 


M =\\ (y+ tos — 2-14.) dx dy = \\(y + 2 F_) dedy. 


Der Beweis fiir M = 2 V aus obiger Gleichung findet sich z. B. bei A. und L. Foppl* 
(S. 68). 

Die plastische Spannungsfunktion ist bei voll ausgebildetem plastischem Zustand 
-ein Hiigel iiber dem Stabquerschnitt von konstanter Neigung, bei dem das Dreh- 
moment mit Hilfe des eben erwaihnten 2. Satzes von Prandtl leicht zu berechnen ist. 


IV. Berechnung des Drehmoments. 


Das Drehmoment eines Stabes kann mit Hilfe der Gl. (17) und des 2. Satzes von 
Prandtl (M = 2 V) leicht berechnet werden, wobei die Form der Spannungsflache 
_nur noch von der Form des Stabquerschnittes abhingig ist. Es wird hier fiir nach- 
folgende vier Querschnittsformen ermittelt. 


a) Kreisquerschnitt. Fir den Halbmesser r betragt die Hohe h des Spannungs- 
hiigels 


h=r tg B, 
wobei 6 den Winkel einer Erzeugenden des Kegels (das ist nimlich die Gestalt des 
Spannungshiigels iiber dem Kreisquerschnitt) mit der Querschnittsfliche angibt. 


Nach Gl. (17) ist nun tg6 = a, ae bzw. die Hohe h = r- o;,: ee 
Damit ergibt sich das Volumen des Kegels zu 
Vazarh =taro,|/<. 
Und nach dem 2. Satz von Prandtl (M = 2 V) folgt fiir das Drehmoment 


me 
M =F nro,)/<. (19) 
b) Rechteckquerschnitt. Seiten b >a. 
, - 

M = (86 —a)o,)/ <. (20) 

Fir den Sonderfall des Quadrates mit a = b lautet 

as lo 

Se a.|/ (21) 


c) Kreisringquerschnitt. Bei diesem Querschnitt ist die geneigte Fliche nur 
iiber die Ringdicke zu nehmen, wiihrend iiber dem Innenring vom Halbmesser 7 


“ A. und L. Féppl: Drang und Zwang, Bd. II, 3. Aufl. Verlag R. Oldenbourg. 1944. 
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eine ebene Flaiche zu spannen ist. Der Spannungshiigel ist somit ein Kegelstumpf 
(s. etwa Foppl', S. 79ff.). 
yet 


M = 32(R —1)o,)/ <. (22) 


d) Gleichseitiges Dreieck. Seitenlinge a. 


as ; 
M = 554.) > (23) 


In diesen Gleichungen sind zwei schon veréffentlichte Sonderfalle von Nadai! und 
Schaden® enthalten. Mit c = 1 ergibt sich der Fall der plastischen Verdrehung eines 
ideal plastischen Stoffes, und zwar nach der Grenzbedingung von Mises! und 
Hencky!. Diesen Fall hat Nadai untersucht. Mit c = 3 ergeben sich die Ausdriicke 
der Verdrehungsfestigkeit spréder Stoffe, die Schaden untersucht hat. 


V. Vergleich mit Versuchsergebnissen. 


Nachfolgend werden diese theoretischen Ergebnisse mit den aus Versuchen er- 
haltenen verglichen. RoS und Eichinger! haben Verdrehungsversuche an Rohren 
und vollen Zylinderstiiben aus Siemens-Martin-Stahl durchgefiihrt. Aus dieser Ver- 
suchsreihe sind die folgenden Daten entnommen. 

a) Hohle Stabe nach Abb. 13 und 22 von Ro’ und Kichinger. 

x) Bei einer Verzerrung von 3°/9) ergibt sich nach Abb. 13 

6g = 2717 kg/cm?, o, = 2590 kg/em?, 
daraus folgt 
¢ =='G,/0, = 1049. 


Der AuBenhalbmesser des Rohres betrug R = 1'5cm, der Innenhalbmesser + = 1:25 cm. 
Damit errechnet sich nach Gl. (22) das Drehmoment 


aes 4860 kgem. 


M = 2x (1°5 — 1-258) 2590|/ 


Nach Abb. 22 von RoS und Eichinger wurde gemessen 

M,; = 4500 kgem, Unterschied zum theoretischen Ergebnis — 13%. 

My, = 4460 ,, , i : " Qf jaeataes 

Die gute Ubereinstimmung der beiden Ergebnisse von Rechnung und Messung 
ist darauf zuriickzufiihren, da sich schon bei geringer Verzerrung von 3°/, der 
plastische Zustand iiber den ganzen, verhaltnismabig diinnen Kreisringquerschnitt 
ausgebreitet hat. Bei vollen Staben und geringer Verzerrung ist die Ubereinstimmung 
hingegen noch nicht vorhanden, wie aus nachfolgender lit. b hervorgeht. Eine dies- 
beziigliche Nachpriifung, etwa nach der Atzmethode von A. Fry?*, liegt nicht vor. 

8) Bei einer Verzerrung von 18°/)) ergibt sich nach Abb. 13 


oq = 3240 kg/em?, 9, = 3110 kg/cm’, 


daraus folgt 
¢ = ¢,/0, = 1:042. 


Nach Gl.. (22) betragt das Drehmoment 
1°042 


M =n (1'58 — 1-259) si10]/ — 5460 kgem. 


Nach Abb. 22 von Ro& und Hichinger wurde gemessen 


15 M. Ro& und A. Hichinger: Versuche zur Klarung der Bruchgefahr. Bericht Nr. 14 der 
EMPA, Ziirich, 1926. Siehe auch Bericht Nr. 172 der EMPA, Ziirich, 1949. 
16 A, Fry: Kruppsche Mh., Juli 1921; Stahl u. Eisen (1921). 


Ingenieur-Archiv VII, 4. 21 
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My; = 4960 kgem, Unterschied zum theoretischen Ergebnis — 10%. 

Die etwas groBere Abweichung vom gemessenen Wert ist bedingt durch die zeich- 
nerisch offenbar nicht mehr so exakte Darstellung im betrachteten Bereich. 

b) Volle Stabe nach Abb. 12 und 29 von Ro’ und Eichinger. Mit Riicksicht 
auf die Verdrehung, die in Abb. 29 ohne Verfestigungsbereich angegeben ist, wird 
eine Verzerrung von nur 2°/,, in Betracht gezogen. Die entsprechende Druck- bzw. 
Zugbeanspruchung betragt 

og = 2690 kg/cm’, 0, = 2630 kg/cm’, 
daraus folgt 
6 =,04/6;, == 1023: 


Der Halbmesser des Stabes T9 betrug r = 1°5 cm, des Stabes T6 r = 1:°0cm. Nach 


1:023 
Gl. (19) ergibt sich das Drehmoment fiir den Stab T9 M = 2 mt 1°58: 2630]/ : Ss 
= 10850 kgem, fiir den Stab T6 M =+ 2 1-08- 2630]/ r0?s _ 3220 kgem. Aus 


Abb. 29 von Ros und Eichinger wurde gemessen 

Mr = 8470 kgem, Unterschied zum theoretischen Ergebnis — 281%, 

Mr, = 2450 ,, , ; i ne Bot geet eee 

Diese gréBeren Unterschiede ergeben sich erwartungsgema8 dadurch, da sich 
der vollplastische Zustand, der dieser Arbeit zugrunde liegt, wegen der geringen 
Verzerrung von nur 2°/), noch nicht ganz iiber den Versuchsstab ausgebreitet hat. 
DaB die berechneten Werte gréBer sind als die gemessenen, weist auf die Richtigkeit 
der Theorie hin. 


c) Zum SchluB sollen die gefundenen Gleichungen noch mit neueren Versuchs- 
ergebnissen von Ro§ und Hichinger verglichen werden (s. Anm. 15, Bericht Nr. 172). 
In Tabelle VII von Ro’ und Eichinger sind Versuchsdaten an GuBeisen EK 50 
angegeben. Die Druck- bzw. Zugfestigkeit betragt 
Op 19 = 8330 kg/em?, oz, = 2680 kg/cm? 
und damit 
¢ = oplozg = 311. 
Nach Gl. (12) ergibt sich die Verdrehungsfestigkeit fiir ¢ = 3 


ty = +0, Vz = + o, = 2730 kg/em*. 
Nach Tabelle VII von Ro§S und Kichinger betrigt 
Ty = 2540 kg/cm?, Unterschied zum theoretischen Ergebnis — 7°5%, 


was eine befriedigende Ubereinstimmung ergibt. 

In gleicher Weise ist auch die Tabelle I von Ro’ und Hichinger’ aufschluBreich. 
Ks wird zunichst die Bruchgrenze der auf Verdrehung beanspruchten Stiibe unter- 
sucht und nur jene Zugfestigkeiten oz mit den Verdrehungsfestigkeiten ty verglichen, 
bei welchen den groBen Bruchdehnungen grofe Schiebungen entsprechen. So erhalt 
man fiir Flu8stahl in der Spalte ty/oz folgende Unterschiede: 


In der Zeile Nr. 1 + 5%, 


GS OP WwW 
— 
w 
=< 


— 5% usw. 
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Zu diesen Zahlenwerten, von welchen die meisten brauchbare Ubereinstimmung 


zeigen, ist zu bemerken, da sie ohne den Beiwert + nach Gl. (12) ermittelt wurden. 


Fir Chrom-Nickelstahl entnimmt man Tabelle I, Zeile Nr. 23, in der Spalte op 
die Werte 


Op = 8340 kg/em?®, oz = 7900 kg/cm, 
daraus folet 


C = dp/oz = 1°056. 
Nach Gl. (12) ergibt sich die Verdrehungsfestigkeit zu 


ees 7900]/ we = + 4680 kg/om*. 


Gemessen wurde tp = 4500 kg/cm?, also mit — 4°, Unterschied zum theoretischen 
Ergebnis. 


Ebenfalls fiir Cr-Ni-Stahl entnimmt man Zeile Nr. 29 in der Spalte op die Werte 
Op = 7650 kg/em?, of = 7500 kg/cm, 


daraus folgt 
¢ = op/oz = 1°02. 
Nach Gl. (12) ergibt sich die Verdrehungsfestigkeit zu 


i dat’ Be 7500]/ 12 — + 4380 kg/om?, 


Gemessen wurde tp = 4400 kg/cm?, also mit + 0°5°% Unterschied zum theoretischen 
Ergebnis. 

Die Ubereinstimmungen zeigen, da& die Grenzbedingung Gl. (6) 
faihig ist, mit Hilfe der neu eingefiihrten Zahl c =o,/o, als Parameter, 
die Verdrehungsgrenzbeanspruchungen von der Plastizitatsgrenze bis 
zur Bruchgrenze zu erfassen. 


VI. Der Grenzwinkel und die Spannungen in der Grenzflache. 


Den Grenzwinkel nennt man allgemein Gleit- bzw. Bruchwinkel, die Grenzfliche 
allgemein Gleit- bzw. Bruchfliche. Nach Torre!” kann der spitze Grenzwinkel 2 « 
zwischen zwei Grenzflichen, die Normalspannung o und die Schubspannung 7 in 
diesen Flachen berechnet werden. Wenn die Grenzbedingung als Funktion der Haupt- 
normalspannung 


0, = 0; (G3) (24) 
gegeben ist, dann lauten die vorerwaéhnten GréBen"” 
tg? « = do,/do; = o,' (25) 
und 
= a, wa ena ek (26, 27) 
Ferner lautet der mathematische Ausdruck fiir die Belastungsart 
Og = 101 +9 Gs; mig. (28) 


Unter Belastungsart versteht man, wenn eine Art der Beanspruchung (Verdrehung, 
Zug, Druck usw.) sowohl unter atmospharischen Verhaltnissen durchgefiihrt wird, 
als auch unter erhdéhtem allseitig gleichem Druck. Zum Beispiel Verdrehung unter 


17 C@. Torre: Osterr. Ingenieur-Arch. 1, 36 und 316 (1946). Auf Seite 336 dieser Zeitschrift 
wird darauf hingewiesen, welche physikelische Bedeutung der Winkel 2 « und die Spannungen o 


und t haben. 
21* 


308 W. Starnberg: Verdrehung bildsamer Metallstiibe itber die FlieBgrenze. 


atmosphirischem Druck und unter allseitigem Druck. Das wird Verdrehungs- 
belastungsart genannt, die sich nach Gl. (28) zu 


1 
r=q=F5 
ergibt. Dann betrigt die mittlere Hauptnormalspannung nach Gl. (28) 
O, = (1 + 4). (29) 
Herrschen nur die tiblichen atmosphirischen Verhiltnisse, welche durch Verschwinden 


des allseitigen gleichen Druckes (o4 p = 0) gekennzeichnet sind, so betraigt dann 


02 = 04p = V (30) 
bzw. nach Gl. (29) 


01; = — G3 = Ty. 


Das verdeutlicht am besten die Mohrsche Darstellung gemiB Abb. 3. Ist o, = 
= o4p <0, dann erhalt man in Abb. 3 die drei Spannungskreise links vom Ursprung, 


Abb. 3. Zur Verdrehungsbelastungsart bei o, = o 


‘Ap = 9 bzw. 6, = 6,)< 0. 


die den im Ursprung gezeichneten ahnlich oder gleich sind. Es sei noch erwahnt, 
daB sich nach Gl. (28) die Druckbelastungsart mit 7 = 1, hingegen die Zugbelastungs- 
art mit r = 0 ergibt. 

Durch Einsetzen von Gl. (28) in Gl. (6) erhalt man 

(o, — 03)? (1 — rq) + (¢ — Ipe, [o, (1 +r) + 03 (1 + 9)] =c0", r+q=1. (31) 
Gl. (28) stellt im Koordinatensystem (03, ¢3, 0,) eine durch die Gerade o, = o, = o3 
gehende Ebene dar. Durch eine Anderung der Parameter r und q dreht sich auch 
diese Ebene um die erwahnte Gerade herum. Wie schon gesagt, stellt Gl. (6) im 
gleichen Koordinatensystem ein Rotationsparaboloid mit der Hauptachse o, = 
= 0, = 0; dar. Man sieht somit, da die Ebene Gl. (28) das Paraboloid Gl. (6) in 
Meridianen schneidet. Die Projektion aller dieser Meridiane, die sich durch Anderung 
der Parameter rv und q ergeben, auf die Koordinatenebene (c,, 03) liefert eben Gl. (31). 

Um die Gl. (25) bis (27) auswerten zu kénnen, ist die Kenntnis der ersten Ab- 
leitung 0,’ = do,/do, erforderlich, die sich aus Gl. (31) ergibt zu 
»__ doy 2 (0; +9) (lg) — (¢— 1) a, (1 +g) 


OL day = BXa aed) Cheerios) Gp liek bikie Tea (32) 
Wie erwihnt, ergibt sich die Verdrehungsbelastungsart fiir r= q = > bzw. fiir 
02 = ee = 0. Die reine Verdrehung folgt dann aus o, = 0 baw. nach GI. (11) und (12) 


aus 6, = —0,=1,= + a,|/ Daraus ergibt sich nach Gl. (25) und (32) 


A. Slibar und K. Desoyer: Zur Schwingungstilgung bei Sternmotoren. 309 


Ue" og =o a ; (33) 


Aus Gi. (26), (27) und (33) bestimmen sich nun die Spannungen in der Gleit- bzw. 
Bruchflaiche zu . 


oT 1 a 
oe a Ce o.|/ a Ol ee ge (34, 35) 
Man sieht, daf fiir ¢ = 1 nach Gl. (33) tg? x = 1 baw. der Grenzwinkel (zwischen 


zwel Grenzflachen) 2« = 90° ist. Bei dem vorliegenden idealplastischen Stoff be- 
tragen die Spannungen in der Grenzfliche nach Gl. (34), (35) 


0, 
y3- 


Das ist die bekannte maximale Schubspannung nach der sog. Gestalt&énderungstheorie. 


o = 0, T= + 


ee 
oad iB 
Abb. 4. Gleitlinien bei verschiedenen Werten c = 1, 2, 3. 


Ein zweiter Sonderfall ergibt sich fiir c= 3. Hier betraigt der Bruchwinkel 
tg? « = 0 bzw. 2x=0. Dieses Ergebnis besagt, daB es sich hier um den Trenn- 
bruch handelt. Die Spannungen lauten dann erwartungsgemaB 


C— Oe T= (es 


was auch der Beanspruchung beim Trennbruch entspricht. Abb. 4a und 4b ver- 
anschaulichen diesen Zustand. Die Gleitlinien I, II in Abb. 4a sind bekanntlich 
Hauptschubspannungstrajektorien, wahrend die Bruchflache in Abb. 4b zur Zug- 
spannung o, senkrecht steht, so dal o = o,, das heiBt, da die Normalspannung in 
der Gleitfliche gleich der Hauptnormalspannung ist. 

Ein allgemeiner Fall ergibt sich fiir c = 2 bzw. nach Gl. (33) 2 « = 52°, s. Abb. 4c. 
Dort befindet sich die Lage der Gleitlinien zwischen dem idealplastischen Fall 
(c= 1, 2« = 90°) und dem spréden Fall (c = 3, 2« = 0). 


(Hingegangen am 2. Juli 1953.) 


Zur Schwingungstilgung bei Sternmotoren. 
Von A. Slibar und K. Desoyer. 
Mit 3 Textabbildungen. 


Zusammenfassung. In der vorliegenden Arbeit werden die Verdrehungen einer elastischen 
Welle untersucht, die mit Pendeltilgern ausgestattete rotierende Massen verbindet. Diese Massen 
stehen im allgemeinsten Fall unter der Einwirkung von Stérmomenten beliebiger Kreisfrequenzen. 
Die Tilger werden dabei als physikalische Pendel eingefiihrt, die Welle wird als elastisches 
Kontinuum behandelt. 

Im weiteren wird die Spezialisierung fiir den haufigsten Anwendungsfall, den Sternmotor 
mit Luftschraube, vorgenommen und die nach Tilgung des Einflusses einer bestimmten ausge- 
wihlten Stérfrequenz verbleibende Torsion der Welle untersucht. Ein beigefiigtes Diagramm 
erlaubt die Bestimmung der Wirkungen der tibrigen Stérfrequenzen auf die Welle, 


310 A. Slibar und K. Desoyer: 


Summary. The paper under review deals with the deformations caused by the twisting of 
an elastic shaft, connecting rotating masses, provided with oscillating absorbers. In the most 
general case these masses are acted upon by moments of derangement of cyclic frequences of 
any kind. In this connexion the absorbers are introduced as physical pendulums, whereas the 
shaft is treated as a continuous, elastic body. 

Further, the author concentrates on the mode of application, most frequently occurring, Viz. 
the star motor with the air propeller, whereby the torsion of the shaft is examined remaining 
after the absorption of the influence of a specified, selected derangement frequency. A graph 
adjoined permits of determining the effects on the shaft of the remaining frequences of derangement. 


Résumé, Cet ouvrage traite du probléme des déformations d’un arbre élastique dues aux 
torsions, arbre dont on se sert pour joindre des masses rotatoires munies de dispositifs & absorption 
pendulaire. Dans le cas le plus général ces masses se trouvent sous l’influence de moments de 
dérangement (perturbation) de fréquences cycliques quelconques. Les dispositifs d’absorption 
sont introduits sous la forme de pendules physiques, tandis the l’arbre est traité comme représentant 
une continuité élastique. 

Puis auteur se concentre sur le cas d’application le plus fréquent, c’est a dire sur le moteur 
a étoile avec hélice, pour examiner la torsion de l’arbre restant aprés l’absorption de l’influence 
d’une fréquence choisie de dérangement specifique. Le diagramme joint au traité permet de 
déterminer les effects sur l’arbre des fréquences de dérangement restantes. 


I. Einleitung. 


Zur Begrenzung der Torsion der Kurbelwelle bei Systemen Sternmotor-Luft- 
schraube werden hiufig Pendelschwingungstilger verwendet!. Nach den einfiihrenden 
Arbeiten von E. 8S. Taylor?, O. Kraemer?, W. Schick!, A. Kimmel und I. Lutz- 
weiler® ist bei dem als mathematisches Pendel aufgefaBten Tilger zur Beseitigung 
des EKinflusses der Grundfrequenz der erregenden St6rmomente die als Abstimm- 


bedingung bekannte Beziehung z 
Oe ait (1) 


einzuhalten. Darin bedeuten 2 die Kreisfrequenz des erregenden Stérmomentes, 
LI die Entfernung des Anlenkpunktes des Pendels von der Drehachse der rotierenden 
Masse, w, die mittlere Winkelgeschwindigkeit der umlaufenden Massen und / die 
Lange des Pendels. Die Beziehung (1) nimmt die einfachere Form 

L 

TF =HP On 
an, wenn die Kreisfrequenz 2 des erregenden Momentes als ein k-faches der mittleren 
Winkelgeschwindigkeit w, eingeftihrt wird. Erst die nichtlinearisierte Behandlung 
des Tilgers an einer Drehmasse liefert, wie von den Verfassern gezeigt wurde®, Aus- 
sagen tiber die GréBe der giinstigsten Masse des Pendelkérpers. 

In der vorliegenden Arbeit soll die Behandlung von Systemen mit Fliehkraft- 
tilgern bei Beriicksichtigung der als elastisches Kontinuum aufgefaBten Welle durch- 
gefiihrt werden. Wa&ahrend bei der nichtlinearisierten Behandlung des Tilgers die 
Forderung nach modglichst gleichférmigem Lauf des gesamten Systems unter dem 
Kinflu8 einer einzigen Erregerfrequenz im Vordergrund stand, wird hier weitgehende 
Verkleinerung der Verdrehungsbeanspruchung der Kurbelwelle angestrebt. Von 
méglichst allgemeinen Bedingungen ausgehend, soll dabei im Verlaufe der Rechnung 
auf den Sternmotor mit Luftschraube spezialisiert werden. 


! Zum Beispiel H. Schrén: Die Dynamik der Verbrennungskraftmaschine. Wien: Springer- 
ved ieee — K. Haug: Die Drehschwingungen in Kolbenmaschinen. Berlin: Springer- 

2 E.S. Taylor: J. Soc. automot. Engr. 88, 81 (1936). 

* O. Kraemer: Z. Ver. dtsch. Ing. 82, 1297 (1938); 83, 901 (1939); Motortechn. Z. 1, 3 (1939) 

4 W. Schick: Ingenieur-Arch. 10, 303 (1939). ; 

* A. Kimmel und I. Lutzweiler: Ingenieur-Arch. 12, 100 (1941). 

° K. Desoyer und A. Slibar: Ingenieur-Arch. 21, 208 (1953); 22 (1954). 
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Il. Die gegebene Anordnung und deren Bewegungsgleichungen. 


Im Sinne obiger Bemerkung legen wir ein System nach Abb. 1 zugrunde. 

Zwei rotierende Massen sind durch eine elastische Welle verbunden. An den Massen 
seien zwei als physikalische Pendel autgefaBte Tilger angelenkt. 

Wir verwenden die folgenden Bezeichnungen: 
vi v2 --- Tragheitsmomente der rotierenden Massen Fl CUR ae 
.. Halbmesser der Welle, 
. Lange der Welle, 
.. Gleitmodul des Wellenwerkstoffes, 

@ ... Massendichte des Wellenwerkstoffes, 

Ly, L,... Anlenkradien der Tilger ,,1‘‘ und ,,2“, 


oS 


My, mM, ... Massen der Tilgerkérper, 
81, 8g... Entfernungen der Tilgerschwerpunkte von deren Anlenkpunkten, 
tg ts2--- Tragheitsradien der Tilgerkérper beziiglich ihrer Schwerpunkte, 
Y1, 2 -.. Drehwinkel der rotierenden Massen, 


9,, 0, ... relative Ausschlagwinkel der Pendelkérper, 
y (z, t) ... Drehwinkel der Welle in bezug auf ein raumfestes Koordinatensystem. 
Alle Winkel werden in dem aus Abb. | ersichtlichen Sinne positiv gezahlt. 

Auf die rotierenden Massen ,,1‘‘ bzw. ,,2‘‘ wirken iuBere Momente, die durch 
abgebrochene Fourier-Reihen 


DP 
M,= M,,+)>) (My, sink gy, + M,,* cos kg), 
k=1 


(2) 
D 
bzw. M,= M,, + >) (M,, sink vy, + M,,* cosk @,) 
baa 


dargestellt seien’. Wir zahlen 
beide Momente im Sinne der 
wachsenden Winkel g, und 9g, 
positiv. Fiir den mit y (z, t) 
bezeichneten Drehwinkel der 
Wellenquerschnitte gilt 


oy 1 dy 
: ez? u* ot (3) 
mit ile 
= |/— 4 
u | ; (4) 


und den Randbedingungen 


und Pip) es ) Abb. 1. Schema der behandelten Anordnung. 
y(t) =2(t). — 5) 

Die Lagrangeschen Bewegungsgleichungen fiir die beiden rotierenden Massen mit 

den angelenkten Pendeln lauten: 


a (Spt) — ae = a + ZOOS), (6a) 

‘i . (Fs ee pies (6b) 
un q 

ae (Spt) — et = a — 39 (ZI, (7a) 

ae eet (7) 


3 Der Einflus der Schwerkraft wird wie tiblich nicht beriicksichtigt. 
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Darin bedeutet 7’, bzw. 7’, die kinetische Energie der betreffenden Masse mit ange- 
lenktem Pendel. In den Gl. (6a) und (7a) stellt der jeweils letzte Term das von der 
Welle auf die Masse ,,1“ bzw. ,,2“° tibertragene Drehmoment dar. 

Fiihrt man in die Gl. (3) zur Ermittlung der stehenden Wellen den Ansatz 


; Raga y (2, t) = Zig) vethoet (8) 
ein, so ergibt sich 
BL ke? we? 
. Gets ci, ae Z (9) 
mit der Losung 
y (z,t) = >) (A, sink w,t + B,cosk wot) (sin oes z+C, cos oe :) Lwot+xz. (10) 
k=1 


Die beiden angefiigten partikuliren Lésungen beriicksichtigen die konstante mittlere 
Rotation der Anordnung und die rein statische Verdrehung der Welle unter der Be- 
lastung mit den konstanten Anteilen der Momente. x ist eine noch zu ermittelnde 
Konstante. Die fiir die beiden letzten Terme der Gl. (6a) bzw. (7a) bendtigten Aus- 
driicke ergeben sich aus (10) zu 


é ony? . 
(4) =» a (A, sim kat = B; cos k wy t) a Ke; 


U U 


| 
Lan 
| 


(+) =» £00 (4, sin k wt + B, cos kort) (cos Ee pune sete 7 pix 
‘a 


Die kinetische Energie 7, der Masse ,,1‘‘ mit angelenktem Pendel lautet: 


Ce 2 m if a . 
T= Pt het, + ts, (G1 + (12) 


Aus den Koordinaten xs, und ys, des Pendelschwerpunktes in dem nicht mitrotierenden 
Koordinatensystem der Abb. 1 


tg, = L, cos p, + 8, cos (py, + 44), 
Ys, = L, sin g, +8, sin (p, + 94) 
erhalt man das Quadrat der Absolutgeschwindigkeit vs, des Pendelschwerpunktes 
Ug? = &g2 + Ys? = (L? + 8? + 2 L, 8,008 9] ge +s? b? Bz 


WHR ee 13 
+ 25 9,8, + 2L, 8, 9,9; cos I. St 
Damit folgt aus (12) und (13) 
Ce vere . . 
T= wit = [L,? + 8? + 2 L, 8, cos 8] 2 + _ sP OP + ms, 9, + (14) 


ay 3 eee ae 
+ my, Ly 8 9, 9, cos J, + ae ts (Pi + A). 


Die kinetische Energie 7’, der Masse ,,2“ mit angelenktem Pendel erhailt man, wenn 
man in (14) den Index 1 gegen den Index 2 austauscht. 

Durch Kinsetzen von 7', bzw. 7’, und deren Ableitungen sowie der Ausdriicke (2) 
und (11) in die Beziehungen (6a), (6b), (7a) und (7b) erhalt man 


(31 + m, Dy? + 2m, Ly s, cos 9,) P + Mm, (87? + ds”) (P, 3 34) = 
— 2m, L, 8, 9, 9, sin 9, — m, L, s, o. sin 3, + 


oe Pp 
+ mL, 8,9, cos 0; = My) +3’ (M,, sink y, + M,,* coskg,) + (15a) 
ke 


co 
rn k ; 
+3 Gat! '—™ (A, sin kwyt + B, cos k wot) +x 
&=1 


> 
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(sy? + 43,2) (G1 + d,) + Ly 8, (p, cos 0, + p,? sin 9,) = 0, (15b) 
(Se + m, Li + 2m, L, 82 cos 34) Ps + my (83% 4- Ug,”) (%, + 3.) 7 
— 2m, Ly 8. 9,9, 8in 0, — ms Ly 8, 9,2 sin 9, + 


2 P 
+ My Ly 8,92 cos 9, = My, +>) (M,, sink o, + M,,* cos k gp) — (16a) 
k—1 
It — k : fF 
= Fat ay Agsin k wot + B,, cos k wy t) |cos ie 1— C,, sin aes i +- “h 
=1 
(83° + tg,?) (Pz + 3,) + Lg 82 (Pz cos B, + —,? sin 0.) = 0. (16 b) 


Wir setzen zur Linearisierung obiger Gleichungen voraus, daB die relativen Aus- 


schlagwinkel #, und #, klein bleiben. Die obigen Gleichungen nehmen dann die 
Form an: 


(S1 + my Ly? + 2m, Ly 81) py + my (8? + tg.) (Pi + d,) + my, Ly 8, d, =...,  (15a’) 
(8° + t52) (@, +O) + Lis: @ + 92 6) = 0, (15b’) 
(So + me Le? + 2m, Le 8) Py + my (8.2 + ¥s.°) (P2 + om) +m 5 Lig 89 J, =. .ay  (16a’) 
(82° + ts,2) Ba + Bs) + Less Ba + G2? Ps) = 0. (16b’) 


III. Auswertung der Bewegungsgleichungen. 


Zerlegt man die Winkelgeschwindigkeiten w, = g, und w, = ¢, in den konstanten 
Anteil w) und die Stérung 4, bzw. 4a,, schreibt also 


0, = @, + Aa, 
1 0 1 | (17) 


WM, = Wy + As, j 


dann sind fiir stationiren Lauf der Anordnung wegen w, = 0 die GréBen w, bzw. a, 
durch Aw, bzw. Aw, zu ersetzen. 
Wir fiihren die dimensionslosen Groen 


ik i Lies Se eS ee ron pe) 
‘are = opr Nis == OF x, = Sy x Oy? aa v0 
Pita weee ie ek so ceniiet * (k=1,2 ) l 
2 of vko Se vk > 4>-+++ DP), (18) 
It 
3 Gut 


——_— =A,, wt=t (y=1 oder 2) 
Sy Wo 

ein und kennzeichnen die Ableitung nach rt durch einen Strich. Mit der Zerlegung (17) 

und den Abkiirzungen (18) gehen die Gl. (15a’), (15b’), (16a’) und (16b’) tiber in: 


\ Aw tr 
(1 +O, + O, 0; Ay + 201 04) > + 101 (A +1) 4, = 


» re - 
= K,, +>) (Kixsinkr + K,,* coskt) + A, ee (A, sink rt + B, coskt) + “| 
k=1 k=1 
(15a’’) 
Aw,’ Ao” + A " 
ae, coat (15b”’) 


Aw,’ yy 
(1 + O, + Os oy dy + 2 Oy 02) =I + O00 (ta + 1)% a 


Paes Pp 
= K,, +> (Ks, sin be -- Ky,* cos.k7) — 16a") 
k=1 
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fore) k ey, 
WA y 4% (4, sin kt + By coskt) (cos a i — C;,sin = } mae 
k=1 
Aw,! dy Oo” + Oy i 
Wo its _ Aes ai 1 ; oe 


wobei wie iiblich die GréBe Aw, in der Kombination (w, + 4a,) gegen w, gestrichen 
wurde. Setzt man (15b’’) in (15a’’) und (16b”) in (16a”) ein, so erhalt man zwei 
Differentialgleichungen fiir #, und 9#,: 


a, 9," + 68, = Ky + A, +2 |(K “fee dly, 
at 


- A,) sin kr + 


(19a) 
ote (x, Pe Jigs B,] cos kel, : 
a0,’ + 6,3, = Ko — Apex + 
ae Dykes = pe i (cos se L— C,sin ee i)| sinkt + (19b) 
k=1 
+ Ess — ig Bi, (cos ee 1—C,sin Hate i)| cos k | 
mit den Abkiirzungen 
1 
a, — (6; Oy ae AS (1 + —,)] : A+ 
1 (20) 


b, = — [14+ 9, + 9,4, (A, + 2)]- (vy = 1 oder 2). 


ied 


Die fiir den uns hier ausschlieBlich interessierenden stationaren Vorgang giiltigen 
partikularen Lésungen lauten 


Saat k wo “ k 
ri Kyo + Ai " ’ Kyy + as ~ Ay, 2 Ky,* + Ay ‘“ By, 3 (21a) 
— . _ ha sin kt se cos kt], 
co . 
# eee Oe 2 A; (c, ay eee i 
za Pas A ay nee Ae U u : 
2= B : _—oo a sinkr + 
k= ‘ 
. ko k wy k aw tb) 
K,;,* + Ay—— B,,(C, sin l os i 
a by ae dy - COS k € 4 


Setzt man (21a) und (21b) in (15b’’) und (16b”) ein, so erhalt man 


Aw,’ Ye See 
oS oad Ae | by 

co 5 (22a) 

— 1—k*®A k , a 

Sp Fema | (x: ko A, — A,} sin kt — (x, — +- A et B,} cos k ‘I 
und * 
Aag !: | Sap ete 
One Ag +1 | by ‘i 
ey iP) ‘ ka ok k 
+ Spee {[Kee + Aomget An (Cxsin “22-1 — cos £2" 1)| sin be + 

) i | 

k ee: 
a. Ege + A, — B,, (Cesin — 1 — cos aa Joos ex} (22b) 


Zur Schwingungstilgung bei Sternmotoren. 315 


Zur Erfiillung der Randbedingungen (5) bilden wir gemaf (10) 


Y2=0 = = 2 (A, sink w)t + B, cos k wy t) Cy, + aot, (23a) 
ee Stal esica uti “81 + O,c0s~2" 1) + wgt + xl, 
k=1 
(23b) 


woraus sich durch Differentiation nach ¢ die Winkelgeschwindigkeiten wm, und w, 
ergeben und daraus durch Ableiten nach t und Division durch w, 


oy’ —— Ao, w 7 
. Sa ig (A, sin kr + B, cos kt) C, (24a) 
und 
, as , co 
= = — = — »’ 2 (A, sinkt + B, cos kr) (sin “ So he» cos “2. 20.1), (24b) 
0 M ts 


Wir setzen die Ausdriicke (22a) und (24a) gleich, ebenso die Ausdriicke (22b) und 
(24b), und fiihren den Koetfizientenvergleich durch. Der Vergleich der Konstanten 
liefert 


und 


Kin Ay x = 0 
(25) 


ein Ergebnis, das man auch sofort aus den Gl. (22a) und (22b) schlieSen kénnte, 
da Aw,’ und Aa,’ fiir stationiren Lauf kein konstantes Glied besitzen kénnen. Da 
fiir stationiren Lauf die Beziehung 


My, + My = 0 (26) 
gelten muB, so ergeben die beiden Gl. (25) mit (18): 
— ee (27) 


Die Konstante x stellt den Verdrehungswinkel der Lingeneinheit der Welle bei stati- 
scher Belastung durch die konstanten Endmomente M,, und M,, = — M,, dar. 


Der Vergleich der Koeffizienten der Winkelfunktionen liefert die Beziehungen: 


1 pes (2h = 
Ay eu eta, maha ue 
1 ay ay i 
hes ee Sy ; -| Ka + th ie {0s sin = cos i) | = (28) 
csr eaeh ise Giese 20 1) 


und 


1 1—#,, : se ie 2 
Gel oka, (Kas wet se Morel 


1 1— kA, bo ( ko, 
j+1 b,—Ra, [ Kee" aie cea Noh fess 
oe. 7), 


= k? B,, (sin = ew hk veda 
Aus den beiden Gl. (28) bestimmen sich A, und A,-C;, aus den beiden Gl. (29) 
ebenso B, und B,,: C;. 


= ae (29) 
J 
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IV. Die Verdrehung der Welle. 


Der klein zu haltende Ausdruck (vy, — 9), der ein Ma8 fiir die Verdrehungs- 
beanspruchung der Welle darstellt, ergibt sich aus (23a) und (23b): 


p 
 — 91 =>) (I, sin ko yt +1,* cosk wot) + #1 (30) 
mit ; 
gene es ee c.(1 By ees i) 
: (31) 
k 
Peeps rp, o,(1 rege i). 


Setzt man die aus (28) und (29) ermittelten Konstanten A,, A; C;, und B,, Beez 


in (31) ein, so folgen die Amplituden des gesamten relativen Verdrehungswinkels 
der Welle zu: 


bef bh, ) 
=H) Gi, 4 1).(b, == & ayy 
aoe ne 


SOG k wo ie m. 
(i + 1) @,— Bay) (1 — cos u ) aj Ms Kaut-nde ton) 


5 a LS 1—k* A, WRG ct oe | 
— isin Pt Gye, Bay Ke ay Kel} | (32) 


Tr 1 { 1—k? A, : 
gates (Ay + 1) (6; — k? ay) 


Lee pn ko ko 
(A, + 1) ec (1 ens a i) - (A, K,,* — A, Ky,*) — 
_ ko 1— kA, 1 — kh, } 
— k* sin aa i| Ged) (b, — k? a,) Kix (A, + 1) (6, — K? ay) K.x* J 


) 
mit P ' 
1— WA Oo. | 2 tobe gee 
D> ti eel. ae 


1—k*® A, Ka . ka eat 1— kA, : k wp ka, | 
—, £06, Fay 28a em atl Ee lees 11G,=sey ft a ee 


+ ke sin =e i]. (ee ees (33) 


Aus (32) ist folgendes zu ersehen: wirkt an der Masse ,,1‘‘ nur eine Stérfrequenz k’ 
und an der Masse ,,2“ nur eine Stdrfrequenz k’’, so werden die Amplituden J\,., [y*, 
Ly, und J‘,* zu Null, wenn man die beiden Abstimmbedingungen 


]—k? A, =0.. und T—4'"24,—206 
einhalt®. 


V. Spezialisierung fiir den Sternmotor. 


Von Interesse fiir die Anwendung wird der einfacher zu tibersehende Spezialfall 
sein, der durch die Verbindung eines Sternmotors mit Tilger und einer anzutreibenden 


Luftschraube unter Zwischenschaltung einer elastischen Welle entsteht. Fiir diesen 
Fall vereinfachen sich die Gl. (32) unter Beachtung von 


Koy = Kee 6 (se 1 So Sears 


iee—— (0) (34) 
‘ ey 
te =0 und — A 
Ss (+ 1, i a,), 
§ In diesem Fall ist: K,, = K,,* = 0 fir k + 0, + k’, 


Ky, = Ky,* = 0 fir k +0, +k’. 
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zu 
= ME Gh Seay) k wo k wo pa LAOS 
: y= D’ (A, + 1) (6, — ka) (1 cos — 7 ge fg ft ei j (35) 
mit 
3s if 1— k2 Ay : k wo “i k wo i k wy k wo 
ie (A, iu) (b, —k? @) Ay U [ 4: U sin U l + k? cos a == 


(36) 


k k : 
—k [ 4. ° cos ie 1 — k? sin we i}. 


Uu 


Aus der Gl. (35) folgt sofort die Amplitude /’,* durch Ersatz von K,; durch K,,*. 
Durch Erfiillen der bekannten, aus Gl. (35) ersichtlichen ersten Abstimmbedingung 


83" + ty? 1 


1—k?A~,=0, das heiBt a te 


(37) 


lassen sich die durch die Erregerfrequenz k hervorgerufenen Torsionen der Welle be- 
seitigen, da dadurch sowohl J}; als auch /;* zu Null werden. Wir setzen nun fiir die 


weiteren Betrachtungen voraus, da die Bedingung (37) fiir eine ausgewahlte Ordnung k 
erfiillt sei und beschaftigen uns im weiteren mit den durch die restlichen Erreger- 
frequenzen k + k hervorgerufenen Torsionen der Welle. Im betrachteten Falle werden 
die bezogenen Amplituden der noch auftretenden Relativverdrehungen der Welle 


PO Tyth) l 


a Moen a SEE ciey = 38 
i ee Vea, X, + ¥, (1 + @ (1 + o,Z,)] (38) 
mi 
A, |, Beales en even eane ee 7 
xX bedi u u u 
Tae | Wo k wo ci k a, 2 
A, — |1—cos l) + ksin 1 
u u u 
ko ko ko ee, 
As © cos 0 1 — k? sin — l 
Y= o r = ko , ee 
A, “2. (1—cos 1) + kein | 
u u 
und 
7 we(ke + 2)41 


ee — 2 


Bei den GroBen I',, I,.* und Z,, soll der obere Index die Voraussetzung der Erfiillung 
der Abstimmbedingung (37) ausdriicken. 

Die Amplitude W © der von der Erregung mit der Ordnung k verbleibenden 
Relativausschlage der Welle laBt sich damit in der Form 


W = Ve @ 4 Pee) (41) 
schreiben, woraus mit Beniitzung von Gl. (38) folgt 
(Xz + ¥x) 1 
W © ae —, (42) 
co Peer Ray ee sr @,:(1+ 0,2, ) 


Die linke Seite von Gl. (42) 148t sich nun in einem Schichtendiagramm (Abb. 2) als 
Funktion der noch frei wihlbaren Parameter o,, 9, und der dimensionslosen Para- 
meter X, und Y,, welche die Daten der Anordnung, die mittlere Winkelgeschwindig- 
-keit w, und die zu untersuchende Ordnung & enthalten, darstellen. Die Bezifferung 
an den Schichtenlinien der Abb. 2 bedeutet den jeweiligen Wert der GréBe nach 


Pat. (42): 
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Im Rahmen der konstruktiven Moglichkeiten kann nun die giinstigste Kombination 
von o, und @, fiir eine gewiinschte Betriebsdrehzahl und die betrachtete Ordnung k 


Abb. 2. Diagramm zur Ermittlung der bezogenen Amplitude. 


aus der Abb. 2 entnommen werden. Uberdies kann man nach durchgefiihrter Wahl 
von o, und @, den Einflu8 aller weiteren Storfrequenzen ebenfalls aus der Abb. 2 
entnehmen. 


VI. Der bifilar angelenkte Tilger am Sternmotor. 


Aus konstruktiven Riicksichten werden Pendeltilger haufig bifilar nach dem in 
Abb. 3 dargestellten Schema angelenkt. Bei dieser Art der Ausfiihrung ist die Relativ- 
bewegung der Pendelmasse gegeniiber der Hauptmasse eine reine Translation. Der 
Schwerpunkt S, der Pendelmasse beschreibt gegeniiber der rotierenden Haupt- 
masse ,,1“ eine Kreisbahn vom Radius s, um den Punkt A,. Wir bezeichnen jetzt 
mit LZ, den Abstand OA, und den Tragheitsradius des Pendelkorpers beziiglich seines 
Schwerpunktes mit is*. Die tibrigen Bezeichnungen behalten wir bei. 


Der Vergleich der kinetischen Energie 7,* des Systems 
Masse ,,1°‘ mit Pendel 


. i 2 
Ty* = (Si + Mm tg,*?) oe + oe Us? (43) 


mit der entsprechenden Beziehung (12) zeigt, daB im vor- 
hergehenden Rechnungsgang an Stelle von $, der Ausdruck 
(31 + my %s5,**) einzusetzen und die GréBe is, in der bis- 
herigen Rechnung durch Null zu ersetzen ist. Die erste - 
Abstimmbedingung nach Gl. (37) geht dann iiber in die 
Form, wie sie fiir das mathematische Pendel bekannt ist: 


Abb, 3. Schema des bi- 70 Le 
filar angelenkten Tilgers. cles (44) 
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Da in diesem Falle die GréGen A, und o, identisch werden, geht bei Verwendung 
eines bifilar angelenkten Tilgers eine zur Erzielung optimaler Wirkung verfiigbate 
Grobe verloren. Dies ersieht man daraus, daf in (42) die GroBe o, — A, bereits durch 
die Abstimmbedingung (44) vorgeschrieben ist. Der Einflu8 der weiteren Frequenzen 
der wirkenden Stérmomente ist wie vorhin aus der Abb. 2 zu entnehmen. 


(EHingegangen am 22. September 1953.) 


Endliche Deformationen diinner Plattenstreifen 
mit freien Liaingsrandern. 


Von F. Miiller-Magyari, Wien. 
Mit 2 Textabbildungen. 


Zusammenfassung. Analog den Elastikaproblemen eines inkompressiblen Stabes werden 
die entsprechenden Differentialgleichungen fiir die endlichen, aber dehnungslosen Verformungen 
eines diinnen Plattenstreifens (Band) abgeleitet. 


Summary. By analogy with the well-known elastica-problems of a thin rod, the author 
develops the differential-equations for finite, but inextensional deformations of a thin strip, 
the two opposite borders of which are free. 


Résumé. Par analogie avec les problémes connus d’une barre incompressible, l’auteur 
développe les équations différentielles relatives & une bande de faible épaisseur, avec deux bords 
opposés libres, pour des déformations inextensionnelles, mais de grandeur finie. 


I. Einleitung. Problemstellung. 


Seit den Anfangen der Schaffung einer Elastizititstheorie diinner Platten, also 
beginnend mit den Arbeiten von Lagrange, Clebsch, Kirchhoff, E. und 
F. Cosserat usw., um nur einige zu nennen, wei man auch schon, daB die Mittel- 
flache einer diinnen Platte unter Belastung im Bereiche der endlichen Deformationen 
eine Flache bildet, welche sich nur wenig von einer abwickelbaren unterscheidet. 
Wie diese etwas vage Aussage ,,nur wenig‘‘ mathematisch prazise formuliert wird, 
soll hier vorlaiufig nicht zur Debatte stehen, tatsachlich besitzt die diinne Platte 


eine relativ groBe Dehnsteifigkeit 
Dees as (1) 


(# Elastizitaétsmodul, ¢ Wandstarke, » Querkontraktionszahl) gegeniiber ihrer Biege- 
steifigkeit 


me a ary (2) 
welche der dritten Potenz der sehr kleinen Wandstirke proportional ist. Es ist daher 
plausibel, da8 die Platte einer ihr aufgezwungenen Verformung vorwiegend durch 
Biegung und nicht durch Dehnung der Mittelflache ausweicht, vorausgesetzt natiir- 
lich, daB nicht die geometrischen und mechanischen Randbedingungen die reine 
Biegung verhindern, was aber leider, zumindestens in gewissen Zonen, sehr oft der 
Fall ist. 

In dem in der vorliegenden Arbeit behandelten speziellen Fall des diinnen Platten- 
streifens mit freien Langsrandern (Band) nach Abb. 1, welcher an den Querraéndern 
jeweils eine starre Leiste aufweisen soll, 1aBt sich jedoch die reine Biegetheorie, wie 
gezeigt wird, mit sehr guter Annaéherung verwenden. Der Plattenstreifen bzw. Band 
wird die Gestalt einer abwickelbaren Fliche annehmen, wobei die Anfangs- und 
Enderzeugenden durch die Lage der starren Querleisten gegeben sind. Durch jeden 


Punkt der Mittellinie geht eine Erzeugende, welche mit ihr den Winkel > + « (s) 


320 F. Miiller-Magyari: 


einschlieBt, wobei « = «(s) eine erst zu ermittelnde Funktion der Bogenlange s 
der Mittellinie bedeutet, deren Randwerte, wie gerade erwahnt, gegeben sind. Nach- 
dem die Gestalt der abwickelbaren Mittelfliche, einer sogenannten Torse, bereits 
durch zwei Funktionen einer Veranderlichen bestimmt ist (meistens wihlt man dafir 
in der Differentialgeometrie ray ame und Torsion der Gratlinie der Torse), kommt 
es also nach Auffindung von « = « (s) nur mehr auf die Ermittlung einer weiteren 


trerer Rand 


Abb. 1. Der diinnwandige Plattenstreifen mit freien Lingsrandern. Koordinaten und Bezeichnungen. 


Funktion an. Als giinstig wird sich bei unserem Problem die Kriimmung der Mittel- 
linie erweisen, ebenfalls als Funktion der Bogenlinge s der Mittellinie. 


Die Aufstellung der beiden, im allgemeinen gekoppelten, gew6hnlichen Differential- 
gleichungen fiir die sogenannten bestimmenden Funktionen der abwickelbaren Ver- 
formungsflache, natiirlich mit Parametern, welche die Abhiingigkeit von der Be- 
lastung ausdriicken, ist das Ziel der Arbeit. 


Beim analogen Stabproblem, namlich Inkompressibilitét, das heift dehnungslose 
Schwerpunktsfaser, vorausgesetzt, sind bekanntlich drei Funktionen gesucht, z. B. 
zwei Verschiebungen und die Verdrehung des Stabquerschnittes, oder man kann die 
instantanen, Winkelgeschwindigkeiten. des begleitenden Dreibeines, also ~ die’ 
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Krimmungen x; als Quasikoordinaten einfiihren!.2 und erhalt dann sehr elegant 
mit Hilfe der Variationsrechnung unter Beriicksichtigung gewisser Vertauschungs- 
relationen die drei Gleichgewichts-Differentialgleichungen in der Kirchhoffschen 
Darstellung und weiterhin alle gewiinschten Aussagen. In unserem Falle des diinnen 
Plattenstreifens, des Bandes, welches sozusagen einen Stab mit einem sehr lang- 
gestreckten Rechteckquerschnitt darstellt, Breite 26S ¢, verlieren wir also einen 
Freiheitsgrad der Verformung, es sind nur mehr zwei gesuchte Funktionen da, die 
Methode der Quasikoordinaten aber soll sinngemiiB auch hier angewandt werden. 


Il. Einige geometrische Beziehungen. 


Die Resultate dieses Kapitels sind natiirlich jedem Geometer bekannt, im Zu- 
sammenhang mit der mechanischen Bedeutung und Formulierung des Problems und 
zur Kintiihrung der notwendigen Bezeichnungen und Abkiirzungen seien sie hier kurz 
wiederholt. 


Die Abbildung und Abwicklung zweier Flachen aufeinander, nicht nur jener vom 
Kriimmungsma Null, wie in unserem Falle der urspriinglich ebenen Platte, ist schon 
seit eh und je ein wichtiges und interessantes Sondergebiet der Differentialgeometrie 
und weist sehr viele Querverbindungen zu Nachbargebieten und -disziplinen auf 
(man vergleiche z. B. den Artikel von A. Voss in der Enzyklopidie der Mathe- 
matischen Wissenschaften, II] D 6a). In dem vorliegenden Falle des diinnen 
Plattenstreifens mit freien Liaingsrindern sind die geometrischen Beziehungen sehr 
einfach und fast anschaulich und elementar ableitbar. 


Wir setzen die Gleichung der Torse in der Gestalt 
Boas (5) A G(s). a 108 (3) 


an, 2; sind die Koordinaten eines Flachenpunktes beziiglich eines raumfesten 
Koordinatensystems, die Bogenlinge s der Mittellinie und der Parameter / sind die 
beiden Flaichenparameter und die sechs Funktionen x, = 2; (s) und y; = y; (s) sind 
vorliufig unbekannt. Die 2;(s) bestimmen die Leitkurve 4=0 der Torse, welche 
offensichtlich eine geoditische Kurve der verformten Mittelflache sein wird, da sie 
aus einer geraden Linie hervorgegangen ist, die Funktionen y; (s) wiederum legen 
die Richtungen der Erzeugenden im Raume fest. 


Denken wir uns die verformte abwickelbare Mittelflache wieder in die urspriing- 
liche Plattenebene zuriickgeftihrt, so geht jede Erzeugende in eine Gerade tiber, welche 


die Mittellinie unter dem Winkel > + « (s) trifft. Sind nun die freien Langsrander 


x, = +6 parallel zur Mittellinie, wie in unserem gezeichneten Beispiel, so treffen auch 
die Erzeugenden mit demselben Winkel am freien Lingsrand auf und man erkennt, 
daB die Annahme der abwickelbaren Verformungsflache in einer gewissen Randzone 
sicher nicht stimmen kann, denn dort an den freien Langsraindern miissen die 
Biegemomente senkrecht zum Rand verschwinden und das ist nunmehr im all- 
gemeinen infolge der Kriimmungsverhiltnisse nicht der Fall. Diese St6rung hat wohl 
einen EinfluB auf die drtlichen Spannungen, die charakteristische Form der Deforma- 
tion im gesamten Gebiet aber wird davon nur wenig beriihrt (Randstérungsproblem). 


Denken wir uns ferner im unverformten Zustand an jedem Punkt der Mittellinie 
2—0 ein normiertes Dreibein parallel zu den raumfesten Koordinaten x; angeheftet, 


i $2 Funk: Stabilitatstheorie bei Staben unter Druck und Drillung. Osterr. Ingenieur-Arch. 1, 
H. 1/2, 2—14. ; z 

ie Raher: Allgemeine Stabilitatsbedingung ftir krumme Stabe. Osterr. Ingenieur-Arch. 6, 
H. 3, 236—246 (1952). 


Ingenieur-Archiv VII, 4. 22 
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@) 
so wird es durch die Verformung in ein normiertes Dreibein &; (j = 1, 2, 3) tiber- 
gehen und die Koordinaten dieses Dreibeines sind die Richtungskosinusse in bezug 
aul 2; 
k6rperfestes Koordinatensystem 
(1) (2) (3) 


&; &; Sy 
x lo x Oo (4) 
raumfestes 1 2 2 
Koordinaten- Ls py B 9 ps 
system 
vs Me V2 V3 


Zwischen ihnen bestehen die bekannten sechs Relationen und sie lassen sich in der 
(1) 
iiblichen Weise durch die drei Eulerschen Winkel 3, y, y ausdriicken. &,; liegt in 
(2) 
Binormalenrichtung der Mittellinie, €; in der Hauptnormalenrichtung, welche auch 


(3) 
gleichzeitig mit der Flachennormale zusammenfallt und &; ist schlieBlich der Tangenten- 
vektor an die verformte Mittellinie 
Uy = 3; % = fs; Xs = Ys. a iB (5) 
Der Vektor y; der Erzeugenden der Torse ist komplanar mit &, und é; 


(1) (3) 
y; = cos «:&, — sin a: &,. (6) 


Die Krimmungen (instantanen Winkelgeschwindigkeiten des begleitenden Drei- 
beines, bezogen auf das korperfeste Koordinatensystem) x, x, und die Torsion t = x, 
sind 


Os 
P eeay Hias Lael Sale © Oe Bs t Be qnh?? 


(3) 
1) 
< dé, f ) Oks 


Xo == GE a 08 4 Oy Si nare: 6 acety (7) 
1 
dé. (2) Ou 
Pepa eT Ys or ae 
“3 ds Si os > eee >) 


wobei aber in unserem Falle die beiden wichtigen Beziehungen gelten 


Hoe Os, wee (S + x] Se cae das heiBt x; = — x,- tg a, (8) 


3 


welche man sofort aus (7) unter Beriicksichtigung der Frenetschen Formeln gewinnen 
kann, bzw. durch Differentiation von (6). 


Die Hauptkriimmung (Flichenkriimmung) x, in einem Punkt der Mittellinie ist 


4 —— % 1 
0 “cos? (9) 


und die Linge der Erzeugenden, gemessen von der Mittellinie bis zum Beriihrungs- 
punkt mit der Gratlinie 
COS & 


h= Ta (10) 


fag 


SchlieBlich lautet noch der Ausdruck fiir das Flachenelement 
df = (h — a): da-da = (h — A) «+ ds- da = (cos « — Aa’): ds- dd. (11) 


Damit sind alle fiir das Weitere notwendigen geometrischen Beziehungen und Formeln 
angeschrieben. 
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Il. Formulierung des Variationsproblems. 


Die bestimmenden gewéhnlichen Differentialgleichungen fiir die beiden gesuchten 
Funktionen werden mit Hilfe der Variationsrechnung aufgestellt. Fiir die gesamte 
potentielle Energie des diinnen Plattenstreifens unter beliebiger auBerer Belastung 
gilt 

TT gesamt = II, + V -» Extr. (12) 


IT; ist die innere Formanderungsenergie des Plattenstreifens, V das Potential der 
Belastung. Bekanntlich ist V gleich der negativen Arbeit der aiuBeren Kriifte 


V2 TT, (13) 


und driickt sich formelmifig durch Flachen- oder Linienintegrale oder Einzelprodukte 
aus. So ist :-z. B. im Falle der Flachenbelastungen p; = p; (%13; 23) 
X, = Acos x ee. 5 
pied: Di = Pi (X13 Xs) = pi (85 A; &), (14) 
v3 = s — Asin « 
() aie 
huss \ p,; u, df, (15) 


wobei wu; die Verschiebungen eines Flichenpunktes aus der Ruhelage bedeuten 
i= ty — ty = &, + A(y, —Cos.x), 
a= S. oe Lo = % +AYs, : (16) 
ts = fp — ey = (ay — 8) 4 Ay, — sin &): 
Somit wird 
q) L b/cos “ 
= |p, (83.4; «) U; (8; A; Ys; «) (cos « — Aa’) ds da. (17) 
s=0 4=—Dd/cos a“ 
Fiihrt man hier die Integration nach / aus, so erhalt man schlieBlich fiir das Potential V 


einen Integralausdruck 
(1) (1) 
Vie Ve (Se at gtk) (18) 


Die Riickfiihrung auf ein einfaches Integral mit der Integrationsvariablen s ist unter 
Beniitzung der Transformationsformeln (14) immer moglich, wenngleich auch bei 
beliebigen nichtkonstanten Flachenbelastungen p; = p; (s; 4; x) die Ausfiihrung der 
Integration nach 4 auf Schwierigkeiten formaler Art stoBen kann. 


Einen wichtigen Spezialfall stellen angreifende Krafte und Momente an den 
starren Endleisten dar. Man kann sie nach den iiblichen Gesetzen der Mechanik 
starrer Korper auf den Anfangs- bzw. Endpunkt der Mittellinie 2 — 0 reduzieren 
und erhailt unter der Voraussetzung flacher Neigungen, das heift geringer Ver- 
drehungen 0, des begleitenden Dreibeines (Existenz des Potentials auch fiir Momente!) 


(2) 8=L e— 
if = Pee; +M,0,; ) (19) 
ee ae js =0 
also 
3 L 
= = ( [P, a, + Py 2, + Ps (4s — 8) + My + Myx, + Myuglds. (20) 


s=0 
Die innere Forminderungsenergie /7; des Plattenstreifens ist unter der Annahme 
der reinen Biegung a 
sa. | i df, (21) 
3 


if, = 


22% 
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wobei x die Hauptkriimmung in einem Punkt der verformten Mittelflache bedeutet 
(die zweite, Hauptkriimmung ist ja stets Null). Fir die Hauptkriimmungsradien & 
der Punkte einer Erzeugenden der Torse gilt die lineare Abnahme zur Gratlinie hin 


vi s 
R= R,(1— =) (22) 
und daraus folgt fiir die entsprechenden Kriimmungen x = 1/R 
: (23) 


lee a ey 
x) baw R, sind stets die Werte im Schnittpunkt der Erzeugenden mit der Mittellinie. 
Der Ausdruck fiir die innere Formanderungsenergie erhalt somit die Gestalt 


- L b/cos « i L Kb “ T= te: 

C i & COS* & COS* % 
Sess DH a ie SEN eres pany of 2. Ps ee SO il 
T= | 4) “bye Oo Od = oem ee 

8=0 A=—b/cos « s=0 ~~ eos? x 
(24) 
Das Variationsproblem lautet demnach 
b a’ 
( Kb cos? x cos? « ; 
dT sana =i V (s; Ui, 1 +++ V3; o) = | ak ee ae eR eS ds — Extr (25) 
se8 : cos? « 
mit den Nebenbedingungen 
i 7) 
Lee: et we ene $y = My COS 063 
/ OX, ; : ¢ 
t = Bs; Myst oe a ae x, = 0; (26) 
; OX, 5 
Ly = Ys} ces api nM: Hg = =, 1g oO = — Xa SIM % COR X. 


Weil «’/cos? « = (tg x)’ ist, empfiehlt es sich, statt «= «(s) die neue Funktion 
B (s) = tg « einzufiihren, ferner nur die Kriimmung ~, (s) (die Hauptkriimmung der 
verformten Mittellinie) als zweite gesuchte Funktion beizugehalten und man erhalt 
die verainderte Formulierung des Variationsproblems 


 (B; B’) 
L a 
* : b . 1 1+ bp’ 
TT gesamt = V (83 ¥55 Oy 2. Yg3B) | 3 ah  (L + BP SG In mays ds — Extr. 
ee (27) 
mit den Nebenbedingungen 
5 ==" 0 
wie oben (26) a Bp ty OU. (28) 
Versieht man die Nebenbedingungen x,’ = 1; usw. mit Lagrangeschen Faktoren 


Ax; Ag; Ay und x, = 0 bzw. x, + Bx, = 0 mit A bzw. A, so lautet die erste Variation 
* 
BED). YOVn oc RCO ap ap’ Kb ax a,’ 8a 
de —E +] 2 x (pp de F pr ia) : 2 a ra i, (Sr — at) + 
OX, op de i oy Ox — [| ox Ox op 
| dn ( de “a t Ay ( oa a raz 4 Al Ge TB aey ibe )Jas. 
(29) 


Beachten wir jetzt die bekannte Verkniipfung zwischen Geschwindigkeit und Winkel- 
geschwindigkeit 


OKs 
De = (41%, — Hq %) USW. (30) 
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und daher analog durch Vertauschung von s und ¢ 


0X3 . v 
A Pee (% X%. — 2%, X,) USW., (31) 


wobei %,, %, *; als ,,Winkelvariationen‘‘ bezeichnet werden kénnen [analog den 
Formeln (7)], ferner die sogenannten Vertauschungsrelationen 


Oxy Ox, v ¥ 
Sen hel 8 8 Me ee, 
Ox, OX < v 
Os 2 eat M3 Hy — Hy 43, 
Or On 
Fs 3 4 a ae ws 9c 
De Be 1 Ha — Ha 4, (32) 


so wird die erste Variation (29) nach Ausfiihrung der Produktintegration 


M3 r 

Fe tS tor Bh + (Kd mp +A) + Am + Ai 4 

L 
re Be eee Bee es 

8=0 7 

+i ]— K b> 2% + A, as + dp Ba ny ae An, — 45) + 

+g | — K b+ 2, %59 — Ae a — Ap Bi — ty — + Aloe, — Boe) | + 

+g | K D+ x x29 AX, ti | AB x] + 
eee | ee aeons heath 


Durch Nullsetzen der Klammerausdriicke (einschlieBlich der Glieder in aa welche 


ae cs stehen, er- 
hilt man in der bekannten Art und Weise 1. die Differentialgleichungen, welche 
die Gleichgewichtslagen des Plattenstreifens bestimmen, 2. die natiirlichen Rand- 
bedingungen fiir s = 0 und s = L, soweit nicht dort Randwerte vorgeschrieben sind 
und daher die Variationen dort gleich Null sind. Wie man oben sieht, sind es, ab- 
gesehen von den meist trivialen Gleichungen fiir A, usw., vier gewohnliche, aber 
gekoppelte und nichtlineare Differentialgleichungen fiir die unbekannten Funktionen £; 
x,; A; A. Beim Anschreiben mu8 man natiirlich noch die Nebenbedingungen (28) 
beriicksichtigen. ° 

Unter Umstiinden ist es ohne weiteres méglich, A und A zu eliminieren und man 
erhalt wirklich nur zwei bestimmende Differentialgleichungen fiir 6 und x,, wie 
bereits am Anfang dieser Arbeit angedeutet wurde. Dies soll an einem einfachen 
Beispiel illustriert werden. 


bei den als willkiirliche Variationen anzusehenden GroBen 2 OS, 


IV. Einzelkraft und -moment an den Querrindern. 


Die Krafte und Momente werden nach dem Anfangs- und Endpunkt der Mittel- 
linie, also nach s = 0 und s = L, gebracht und nach den drei raumfesten Koordinaten- 
achsen zerlegt gedacht. Sie miissen natiirlich die 4uBeren trivialen Gleichgewichts- 
bedingungen erfiillen. Das Potential wurde bereits in (20) angeschrieben. Die Rand- 
werte des Winkels « = « (s) sind durch die Stellung der Querleisten vorgeschrieben, 
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op 


dort verschwindet also die Variation =- und es ergeben sich aus (33) unter Beriick- 
sichtigung von (20) folgende Differentialgleichungen und Randbedingungen : 
Kb KO! a FF 
pia Ve oe Sages VE ee 
iid aires 
— K 6-7 (4 ¢) i Aa 2 + Ap Be + dy Y2 — AB x, — Ge (AB) + M, Bx, = 9, 


. dA > 
+ Kb: puppy — de 1 — Ap By — Ay V1 ds + A (1 + 6?) x, — M, Bx, — Mzx, = 9, 


Let = f= Mey == 0): (34) 
Gi Te as . ; 
USi= Tees ae as ale (35) 
fur s= 0-und-s = 1: 
hes Po 2 Gp P= 8 et PS (36) 
Kb-x,0 +A, —M.= 6 A= Ig Ae (37) 


L=20 a 


Abb. 2. (a) Verformung des diinnwandigen rechteckigen Plattenstreifens unter reiner Drehmomenten- 
beanspruchung und (6) Abhiangigkeit des Winkels der Erzeugenden von der Liangskoordinate s. 


Betrachten wir z. B. als Spezialfall nur mehr ein angreifendes reines Drehmoment M,, 
wie es in Abb. 2 dargestellt ist. Fiir s = 0 ist « (0) = 0 und daher auch f = tg « = 0, 
in der Mitte des Plattenstreifens s = L/2 = 1 ist « vorlaiufig unbekannt, a (1) = «,, 
das heiBt 6 = B,, aber aus Symmetriegriinden ist (’ (1) = 0. Sofort erkennt man 


aus (34) bis (37 
(34) (37) A, = 4g =, =O: (38) 
Aw M 
AS Mee Se (39) 
Schreibt man diese Gl. (39) um, so erhalt man fiir die Flachenkriimmung 
te = 4, (cos « = — Bed Salata ik 
: ; » ei ee ves _1 + d/h cos o (40) 
1 — d/h cos « 


und man erkennt, da sie nichts anderes als die Gleichgewi tsbedingung zwischen 
den inneren Schnittmomenten lings einer Erzeugenden des Plattenstreifens und der 
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Komponente M,- sin « des iuBeren Momentes darstellt.. Die einzige iibrigbleibende 
Differentialgleichung ist 


B ithe BP dhe 
P 4 2 gi V8 — B as (Se 5) =0, (41) 


welche die Eulersche Differentialgleichung des Variationsproblems 


l 1 
Qs (1) — 0; (0) = 0 — @, (0) = \%, d= —| Bx, ds = 
s=0 s=0 
En Be aacrth 
= Kb \ gp & = KE \ FB: A) ds > Ext. (42) 
= 10) s=0 


: : 1 
ist (Dis auf den* Faktor 3): das heiBt, die zweite bestimmende Gleichung fiir die 


Winkelfunktion £ = f(s) ergibt sich aus der Forderung, daB der Verdrehungs- — 
winkel @, am Ende méglichst groB wird. 


Nachdem der Integrand nicht von der Integrationsvariablen s abhangt, gelingt 
die Integration durch eine Quadratur, es ist 
, p* 2 b? B® 
C= Lea a 
ae le Sra yeaa Cage 
(22h 
wobei c eine Integrationskonstante bedeutet. Sie ergibt sich aus der Bedingung 
foe 0 furs ¢ 20 


(43) 


I By? 


Se ee pee: 44) 
Man trennt also in (43) und erhalt 
1 1 Lise 6 Pai eps)? 
a (aeee — 1) oe = pea ) 


eine Differentialgleichung 1. Ordnung, welche man sehr einfach mit einem kombi- 
nierten numerisch-graphischen Verfahren aufldsen kann. Man beginnt bei s = 0 
mit 6 = 0 und bf’ = 6 =1, erhalt dadurch ein As/b = Ao = Ap/B und schreitet 
so in mehreren Intervallen bis 6 = 6,. Zu einem angenommenen f, ergibt sich dann 
das zugehorige 1/b. Fiir kleine 6 6’ <1, also in der Nahe von s = /, kann man zur 
Ergainzung und Kontrolle auch mit Potenzreihen in (s — 1)? arbeiten. Die linke 
Seite von (45) wird dann einfach 


1 1 1+obpf . b? Bp? 
alee pe 1— 6p’ er], 3 6 (46) 


In Abb. 2 ist der Verlauf von 6 = f(s) fiir ein gewahltes 6, = 0°8 gezeichnet 
(x, = 38°6°). Es wurde eine Naherung mit nur neun Intervallen in f# gewahlt, das 
zugehérige 1/b erhalt man mit 1:33 und der Verdrehwinkel 0, am Anfang baw. Ende 
(halber gesamter Verdrehwinkel auf die Lange L = 21) berechnet sich zu 


l 1°33 


REE ae or BB UB Sas GR 
OD eae Vira —T5F Dap hns- 340 = ee Oi080.M (47) 
a= STP o=0 


Dieses Beispiel wurde bereits von L. Kirste* behandelt, welcher zuerst einmal 
ebenfalls die Gl. (40) benititzt statt (41) usw., jedoch eine etwas andersgeartete 


oa L. Kirste: Elastic deformation of thin plates into developable surfaces. Fourth Congress, 
1953. Internationale Vereinigung fiir Briickenbau und Hochbau. 
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Naherungsbedingung aufstellt, namlich jene der verschwindenden Resultierenden 
der inneren Querkrifte lings einer Erzeugenden im Plattenstreifen. Naturgemaf 
erhalten wir infolge der Forderung (42), aus denen unsere zweite Differential- 
gleichung hervorgegangen ist, einen grdBeren Verdrehwinkel, der Unterschied ist 
jedoch unbedeutend. 

Zum AbschluB sei noch der Grenzfall 1/b —> co angeschrieben: der Winkel « (s) 
ist dann konstant iiber die gesamte Lange gleich 45°, das hei®t 6 = 1 und f’ = 0. 
Die Funktion y in (27) und weiterhin wird in diesem Grenzfall gy = 8 und die 
Kriimmung ~, der Mittellinie ist ebenfalls konstant 

Mee ot 
ara” goo vaeses 

Die Mittellinie nimmt dann die Gestalt einer Schraubenlinie an. Die Verwindung 0, 
(Verdrehung pro Langeneinheit) ist ebenfalls konstant und nach (47) 
Me} 

SOD aes 

lauter Resultate, welche man ohne weiteres auch elementar verifizieren kann. Hs sei 
hier noch bemerkt, wie bereits in Anm. 3 gezeigt wurde, daB die Resultate fiir Langen 
1/b > 2°5 bis 3 schon nicht mehr viel von diesem Grenzfall abweichen (geringer Stér- 
einflu8B der gelagerten Querrander). 


(48) 


6,-0.4 


(49) 


(Eingegangen am 29. September 1953.) 


Der Zykloidenraster. 
Von Dr. R. Bereis, Wien. 
Mit 3 Textabbildungen. 


Zusammenfassung. Es wird em auerst einfaches Zeichenverfahren fiir Zykloiden und Pseudo- 
zykloiden (inklusive deren Ausartungen) angegeben, das auf einer graphischen Integration beruht. 


Summary. An extremely simple drawing method for cycloids and pseudo-cycloids (including 
their degenerations) is explained, which is based on a graphic integration. 


Résumé, L’auteur explique une méthode de dessin excessivement simple pour cycloides 
et pseudo-cycloides (y-compris leurs dégénérations), méthode qui se base sur une intégration 
graphique. 

Einleitung. 

Durch Spitzenkreis ky (O, ro) und Scheitelkreis k, (O, 7,) ist eine Drehschar 
von Zykloiden! festgelegt, die durch jeden reellen Punkt P des Kreisringes (ko, k,) 
zwei reelle Kurven sendet. 7, ist der Fixkreisradius; der Rollkreisradius hat den 


] “ : ; 
Wert r= 4 (71 — 7). Je nach GroBen- und Realitatsverhaltnissen von (ro, 7) baw. 


(79,71) handelt es sich dabei um verschiedene Arten von Zykloiden bzw. deren Aus- 
artungen, wie nachstehende Tabelle zeigt: 


a) Zykloiden mit reellem Fix- und Rollkreis. 


Polare Speerkoordinaten 


1 Epizyklordersx.ccasmii orien hoes = 0 h = r, cos (5 p+ ¢) 

— - nei a _ - | " 

2 Hy pozykloide). 2. ox ccnp am : OF = h = 7, cos (3 yg + o} 
bee he's | AE REL es ssi: Seite ‘ae | : 

3 ‘Kreisevolvente ............. Te OSes Ty ree h=r,(y + ¢) 


* R. Blum: Cykloiden u. Cykloidalen usw., Progr. Wilhelms-Realsch. Stuttgart (1902). — 
R. Bereis: Uber die Béschungslinien auf Drehquadriken. Mh. Math. 56, 344—351 (1952). 
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b) Pseudozykloiden. 


cous LoS a es Se ee ee ee Se ee OS ee ee 


| 
} Eiyperzykloide: ...200...05.. eee eg See eat nc | ki =@.cn (5 yp + 7 
‘aby ee | 
ate. |. 2 
| 
2 POMREAVELOIIG ccs cpt ion a a's sans | pete = @s ti = 0 | h = bsh (F yg + ‘| 
| 
3 | Logarithmische Spirale....... | Grenzfall i h=ezkp+e Tar. 


In oben angefiihrte Tabelle wurden auch die Gleichungen der einzelnen Dreh- 
scharen in Speerkoordinaten? aufgenommen. Bei dieser Darstellungsart wird 
bekanntlich die orientierte Kurventangente durch ihren orientierten Abstand Ot = h 
von O und dessen Polarwinkel p festgelegt. 

Ubt man auf eine solche Drehschar von Zykloiden noch alle méglichen Streckungen 
von der Mitte O aus, so erhalt man eine zweiparametrige Ahnlichkeitsschar, die im 
allgemeinen durch jeden Punkt der Ebene co! Kurven hindurchschickt. : 


Zweck vorliegender Note ist es, auf ein einfaches Zeichenverfahren fiir Zykloiden 
hinzuweisen, wobei von einer solechen Kurve die Mitte O, ein Linienelement (Ff) 
und das Radienverhaltnis als bekannt vorausgesetzt wird. Dieses Verfahren beruht 
ebenso wie ein schon friiher verdffentlichtes? auf einer graphischen Integration. 


Der Zykloidenraster. 


Man hebe auf einem Millimeterpapier eine Grundlinie x hervor und trage sodann 
zwei einander auf x schneidende Gerade a und b ein, deren Steigungen der Bedingung 
teas te Pa 7? tr, 2 ry 
genugen. Sodann lege man ein durchsichtiges Pauspapier, auf dem die Mitte O und 
ein Linienelement (P,¢) der zu zeichnenden Zykloide markiert sind, derart auf den 
»4ykloidenraster’, dag O auf x und P auf die Gerade a zu liegen kommen, 
wahrend # in die z-Richtung weist. Jener Punkt K von b, der mit P auf einer zu x 
normalen Rasterlinie liegt, ist die Kriimmungsmitte fiir P. Nun kann man ein Stiick 
des Kriimmungskreises zeichnen. Auf diesem wahle man sodann einen Nachbarpunkt 
und verfahre mit ihm in gleicher Weise. So fortfahrend, erhalt man die gewiinschte 
Zykloide auf dem Pauspapier durch aneinandergereihte Bogenstiicke von Kriimmungs- 

kreisen approximiert. 


Abb. 1 zeigt die Anordnung des Rasters bei Zykloiden mit reellem Polkreispaar, 
Abb. 2 bei Pseudozykloiden. Bei den Grenzfallen — Kreisevolvente und logarithmische 
Spirale — kann man sowohl der gréferen Genauigkeit, als auch der leichteren Hand- 
habung zuliebe das angegebene Verfahren etwas modifizieren. Bei der Konstruktion 
einer Kreisevolvente kann man das Pauspapier um O drehbar machen, ist doch in 
diesem Fall O K = 7, fiir eine bestimmte Kurve eine konstante GroBe. Beim Zeichnen 
einer logarithmischen Spirale legt man giinstigerweise den auf der x-Achse frei wahl- 
baren Scheitel S nach O, worauf das Pauspapier wieder um O drehbar auf der Unter- 
lage befestigt werden kann. 


2 Die Gleichung einer Kurve in polaren Speerkoordinaten ist im wesentlichen mit der 
magischen Gleichung einer Kurve identisch, nur daB bei der magischen Gleichung der 
Polarwinkel der Kurventangente statt des Polarwinkels des Ursprungslotes Verwendung findet. 
Siehe W. P. Hiern: On the magical equation to the Tangent of a curve. Quart. J. of Math. 6 
(1864). — Die Gleichungen der Zykloiden in polaren Speerkoordinaten wurden von Fr. Fa- 
bricius-Bjerre bekanntgegeben: Uber zykloidale Kurven in der Ebene und im Raum. Dan. 
Mat. Fys. Medd. 26, Nr. 9 (1951). 

3 R. Bereis: A. a. O. 
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Wahlt man zur Beschreibung einer ebenen Kurve polare Speerkoordinaten, so 
lautet, unter Zugrundelegung der iiblichen Vorstellung der GauBschen Zahlenebene, 
die komplexe Kurvengleichung 


z=(h+h’'1)eé?, 
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wahrend ihre Evolute durch 

z= (h' +h" ilie’? 
festgelegt ist. Man erhilt somit den Abstand y, der Kriim- 
mungsmitte A eines Kurvenpunktes P von jenem Polstrahl «, 
der zur Tangente ¢ von P parallelliuft, durch Multiplikation 
des Lotes y= Px mit dem Faktor h’’:h (Abb. 3), also 


hk” 
VA oy aa 
Bei einer Ahnlichkeitsschar von Zykloiden nimmt der 
Faktor h’’:h den konstanten Wert —r,2:7,2 an. Somit be- 
steht zwischen den Triigern aller parallelen Linienelemente der Ahnlichkeitsschar und 


ihren Kriimmungsmitten eine perspektive Affinitit, die bei dem oben geschilderten 
Zeichenverfahren ausgewertet wird. 


(Hingegangen am 25. September 1953.) 


Kin Beitrag zum Maschinenrechnen: 
Die Berechnung vielstelliger Quotienten nach dem Aufbauverfahren!. 


Von K. Holecek, Wien. 
Mit 4 Textabbildungen. 


Zusammenfassung. Sind von einem Quotienten mehr Stellen erforderlich, als auf der Rechen- 
maschine in einem Arbeitsgang erreichbar sind, dann braucht man — falls das subtraktive Ver- 
fahren angewendet wird — nur den jeweiligen Divisionsrest als neuen Dividenden zu betrachten 
und kann durch Weiterfiihrung der Division beliebig viele Quotientenstellen ermitteln. Wird die 
Division nach der fiir die Praxis einfacheren Aufbaumethode durchgefiihrt, dann scheint die 
Weiterfihrung dadurch behindert, daB bei diesem Verfahren die Reste nicht unmittelbar zur 
Verfiigung stehen. Es wird gezeigt, da8 man auch ohne Kenntnis der Reste durch Hinhalten 
eines ganz einfachen Scliemas beliebig viele Stellen des Quotienten ermitteln kann. 


Summary. Supposing you need more digits of a quotient than available on the calculator 
at one shift, you will have to take the remainder as the new dividend — in case the substractive 
method is used — and to carry on the division in order to find as many places of the quotient 
as wished. Whenever you carry on the division through the multiplying method, which is much 
simpler in practice, you will find that there is a certain handicap in going on with the dividing 
process, as the remainders do not appear on the calculator. It may be seen, however, that through 
the latter method, without knowing the remainders, any number of places can be ascertained, 
provided a certain and very simple routine is observed. 


Résumé. A-t-on besoin d’un plus grand nombre de chiffres au quotient que celui que l’on 
pourrait obtenir en un seul procédé au moyen de la machine 4 calculer, il suffit — au cas ot Ja 
méthode soustractive serait appliquée — de considérer le reste de la division comme un nouveau 
dividende et l’on peut en continuant la division obtenir un nombre quelconque de chiffres au 
quotient. Si la division est faite & base de la méthode division par multiplication qui est plus 
simple dans la pratique, la continuation semble alors étre empéchée du fait que par ce procédé 
les restes ne sont pas immédiatement disponibles. Il est démontré que l’on peut, méme si lon 
n’a pas connaissance des restes, obtenir un nombre quelconque de chiffres au quotient en se 
servant d’un schéma tout a fait simple. 


Kinleitung. 


Maschinenrechner, die mit Handkurbelmaschinen arbeiten, verwenden fiir die 
Division mit Vorliebe das Aufbauverfahren, das heift sie stellen den Divisor im 
Einstellwerk (EW) ein und kurbeln bei geeigneter Wagenversetzung bis der Dividend 


1 Diese Untersuchung ‘wurde durch eine Vermutung des Rechenmaschinenkonstrukteurs 
Curt Herzstark angeregt. 


332 K. Holecek: 


im Hauptzahlwerk (HW) aufgebaut ist, wobei das Umdrehungszihlwerk (UW) den 
Quotienten zeigt. Der Vorgang ist also eine Multiplikation und entspricht der 
Divisionsprobe. Das hat gegeniiber dem subtraktiven oder Abbauverfahren den 
Vorteil, daB nach Fertigstellung der Rechnung Dividend, Divisor und Quotient in 
dieser Reihenfolge im HW, EW und UW stehen und kontrolliert werden k6nnen. 
Uberdies erspart man sich das Einstellen, Ubertragen und Léschen des Dividenden 
sowie das Léschen der Eins aus dem UW. Wenn der Dividend von einem friiheren 
Rechnungsgang schon im HW stebt, ist natiirlich das subtraktive Verfahren ginstiger, 
allerdings nur dann, wenn der Rechner schon mit Vorbedacht dafiir gesorgt hat, dab 
der Dividend auch geniigend weit links im HW steht. 


Die Anzahl der Quotientenstellen, die in einem Zuge ermittelt werden konnen. 


Diese Anzahl, die naturgemi® fiir beide Verfahren gleich sein muf, ist nicht 
immer gleich der Stellenzahl des UW. Sie kann geringer sein, wenn der Divisor so 
lang ist, daB man nur einen Teil des UW fiir die Berechnung von Quotientenstellen 


hy 
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Abb. 3. Abb. 4. 


Abb. 1 bis 4. Schema der Zahlwerke zu Beginn der Division fiir verschieden lange Divisoren. 

HW Hauptzahlwerk, EW Einstellwerk, @ @ © ecingestellter Divisor, © @ @ Dividend, e e e un- 

besetzte Stellen, z, Stellenzahl des Dividenden, z, Stellenzahl des Divisors, 2, Stellenzahl des 

Quotienten, nz (= 8) Stellenzahl des Einstellwerkes, ny (= 6) Stellenzahl des Umdrehungs- 
zaihlwerkes, ny (= 11) Stellenzahl des Hauptzihlwerkes. 


ausniitzen kann. Die Abb. 1 bis 4, die schematisch die Einstellung fiir je eine Division 
a: 6 darstellen, zeigen, in welcher Weise die Anzahl der in einem Zuge errechenbaren 
Quotientenstellen von den Stellenzahlen der Zaihlwerke und des Divisors abhiingt. 
Ks bedeuten: z, die Stellenzahl des Divisors, z, jene des Quotienten und nx, ny, ny 
der Reihe nach die Stellenzahlen von EW, UW und HW. Fiir die Bilder wurden 
8 x 6 x 11 Stellen gewahlt, die Stellenzahlen der Kleinstrechenmaschine Curta. 
Abb. | zeigt am Beispiel eines dreistelligen Divisors den einfachsten Fall, der 
so viele Quotientenstellen liefert, als UW-Stellen vorhanden sind (in der Abb. also 6). 
Hier und auch in den folgenden Abbildungen wurde ein solcher Divisor b gewahlt, 
der groBer ist als die gleich viele Ziffern umfassende Zahl a,, die vom Dividenden a, 
von dessen héchster Stelle beginnend, abgetrennt zu denken ist (z. B. a = 312153 
b = 3181, a, = 3121). 
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Die volle Ausniitzung des UW ist nur méglich, solange der Divisor in dem zur 
Verfiigung stehenden Teil des EW Platz hat. Wie Abb. 2, die den Grenzfall zeigt, 
zu entnehmen ist, kénnen nur dann mit GewiBheit n, Stellen vom Quotienten er- 
mittelt werden, wenn fiir die Stellenzahl z, des Divisors 6 gilt: 


Zp S Ny— Ny. (1) 


Falls a, = b (s. oben), wird die Schranke um 1 groBer. 
Abb. 3 zeigt am Beispiel eines 7stelligen Divisors einen Fall, bei dem man das 
UW nicht mehr voll ausniitzen kann. Es ist dann 


Zq = My — % (2) 
oder gegebenenfalls um 1 mehr, wenn a, = b. 

Abb. 4 zeigt die sicherlich immer mégliche Mindestanzahl (z,),,j, von Quotienten- 
stellen. Sie ergibt sich, wenn der Divisor das EW vollstindig ausfiillt, also wenn 
2 = Ny. Es ist (2,)min = MH — Me- 


Die Berechnung weiterer Quotientenstellen bei Anwendung des subtraktiven Verfahrens. 


' Zur Ermittlung weiterer Quotientenstellen braucht man nichts anderes zu tun, 
als den Divisionsrest als neuen Dividenden zu betrachten und mit dem gleichen Divisor 
wieder eine Division durchzufthren. So bekommt man eine weitere Zifferngruppe 
des Quotienten. Der Vorgang kann natiirlich beliebig oft wiederholt werden. 


Es empfiehlt sich, die Reste immer um gleich viele Stellen nach 
vor zu nehmen, niamlich um ebenso viele, als man Quotientenstellen 
in einem Arbeitsgang entwickeln kann?. Dann haben alle Teilquotienten 
genau gleich viele Ziffern, die man ohneweiters aneinanderreihen kann, und es werden 
Unsicherheiten oder Irrtiimer vermieden, falls einmal ein Teilquotient mit Null be- 
ginnt (z. B. 740: 329 = 2°24924 012158. .). 


In der Tat gilt fiir die Division a: 6, wenn q; und r; den 1-ten Teilquotienten bzw. 
den i-ten Rest bedeuten: 


a r 
in i oe an (3;) 
r r 

Gq =e ; ox " bt (32) 
qs = 10" ‘ + , (33) 

Ty — Ty 
17. = MA, ; b (3,—1) 
a 10” T-1 v (3 ) 
YW b v 


Daraus ergibt sich: 


7 = Gi +10 F geet gs Fo RARE goed 


ry 
eNO gtd 2" aes (4) 


— Da nies die Art der Reste natiirlich nicht voraussehen kann, ist fir diese Stellenzahl schon 
am Beginn der Rechnung auf alle Falle der kleinere Wert zu nehmen, also jener, der sich bei 


a, < b ergibt, auch wenn tatsichlich aq = b ist. ; 
"3 Fiir “aie Stellenzahl z, der Teilquotienten steht hier der Einfachheit halber n. 
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Die Formel (4) driickt aber nichts anderes aus, als daf} man die Ziffern der Teil- 
quotienten nur aneinanderzureihen braucht, sofern die Reste um genau » (=25,23) 
Stellen vorriicken, wenn sie die Rolle eines neuen Dividenden tibernehmen. 


Der Dezimalpunkt wird beim ersten Arbeitsgang bestimmt. Bei den weiteren 
Teilquotienten mu8 man sich darum gar nicht mehr kiimmern, da man ja nur die 
einzelnen Ziffernfolgen braucht. 


Der ganze Vorgang ist also denkbar einfach, doch wird es bei der Durchfithrung 
oft als lastig empfunden, da man bei Maschinen ohne Wirteleinstellung jedesmal 
eine Reihe von Schritten auszufiihren hat, um den Rest nach vorne zu bringen und 
den Divisor wieder bereit zu haben. Bei der Weiterfiihrung der Aufbaudivision, die 
im folgenden gezeigt wird, entfallen in der Regel alle diese Schritte, bis auf eine ge- 
meinsame Zaihlwerksloschung. 


Die praktische Ermittlung weiterer Quotientenstellen bei der Aufbaudivision. 


Hier scheint zunachst die Schwierigkeit zu bestehen, da man bei einer Aufbau- 
division a: 6 den Divisionsrest 7 nicht unmittelbar zur Verfiigung hat. Denn entweder 
bildet man im HW eine méglichst knapp unter dem Dividenden a liegende Zahl a’, 
dann ist r = a — a’, oder man iiberschreitet a méglichst wenig, etwa um w, und bildet 
im HW a=a-+u, dann ist r= b—u. Es ware also in beiden Fallen eine kleine 
Zwischenrechnung nétig, um den Rest 7 zu bekommen. Nun kann man aber, wie 
die folgende Uberlegung zeigt, mit dem UberschuB selbst das Auslangen finden und 
braucht den Rest erst gar nicht zu ermitteln. Das ist insofern ein Vorteil, als man den 
UberschuB w in der Regel direkt ablesen kann, wenn der Dividend nicht gerade sehr 
lang ist. 


Kurbelt man auf einen Uberschu8 u, im HW, dann ist 9g, = a + “3 , wenn all- 


gemein q; die Anzeige des UW bedeutet, die sich bei der Bildung eines méglichst kleinen 
Uberschusses ergibt. Es ist gy, = 9g, + 1 — die niedrigste Stelle als Einerstelle voraus- 
gesetzt — und u,;=—b —1,. 


Die Gl. (3,) kann daher geschrieben werden: 


—U, To 


le 
q; = 10° FA b 


oder 


; U Ce 
Fs = LO” (10" ; , (5) 


Gl. (5) bedeutet aber nichts anderes, als da man das Zehnerkomplement von q, 
erhalt, wenn man den Uberschu8 uw, um n Stellen nach links verschoben im HW 
aufbaut und dabei wieder einen Uberschu8 erzeugt. Dieser neuerliche Uberschuf 
(iiber w,) ist gemiB Formel (5) der zweite Rest r,. Die Ziffern von q, = q. erhalt man 
bequem lesbar, wenn man einfach das UW auf Gegenlaufigkeit schaltet*. 


Da der zweite Divisionsschritt bereits den richtigen Rest r, liefert, kann damit 
so wie mit dem Dividenden a zu Beginn der Rechnung verfahren werden. Die Ver- 
haltnisse wiederholen sich also nach je zwei Schritten. 


Fur die Fortsetzung der Aufbaudivision a:b ergibt sich somit die folgende Ver- 
fahrensregel : 


* Bei Maschinen mit automatischer Steuerung der Laufrichtung des UW (gemaB der ersten 
Kurbeldrehung nach der Léschung des UW) kann man die Umschaltung dadurch erzielen, da 
man eine subtraktive Hilfsdrehung macht, die durch eine nachfolgende additive fiir die Zahl- 
werke unwirksam gemacht wird. Das UW ist aber dann bereits ‘auf Gegenlaufigkeit geschaltet. 


Die Berechnung vielstelliger Quotienten nach dem Aufbauverfahren. 335 


1. Der Dividend a wird bis zum kleinstméglichen UberschuB wu, 
im HW aufgebaut, Der (n-stellige) Wert g,, der dann im UW steht, 
ist an seiner niedrigsten Stelle um eine Einheit zu verringern, womit 
sich der erste Teilquotient q, ergibt. 


2. Das UW wird umgesteuert und nunmehr (nach Léschung beider 
Zihlwerke) der UberschuB uw, um genau n Stellen nach links ver- 
schoben aufgebaut, wobei neuerdings ein Uberschu8B zu erzeugen ist. 
Dieser letztere UberschuB ist bereits 7,. Das UW enthalt die richtige 
Ziffernfolge des zweiten Teilquotienten q,. 

Man kann beliebig weiter fortsetzen, immer unter Beachtung des 
obigen Zweierturnus. 

Gl. (4) kann man nun durch q, und q, ausdriicken, wobei es der einfacheren 
Schreibweise halber zweckmiBig ist, die niedrigste Stelle von q, als Einerstelle an- 
zunehmen und den Ausdruck zu diesem Zwecke, falls nétig, mit einer geeigneten 
Zehnerpotenz 10* zu multiplizieren: 


nt = 2n 
er Gel) ae LO gai > [10 * (g, 01) 10 Pg |i 
es + F107 0) (9 = 1) 0 OD 9] 10> © Dae (6) 
Vii. i ke. yageles 


Eine mehr als zweifache Fiillung des UW wird in der Rechenpraxis freilich nur 
selten vorkommen. Ein solecher Ausnahmefall ist z. B. das Quadratwurzelverfahren 
nach Collatz’. 

Zum Abschlu8 médgen die Verhaltnisse noch an einem Zahlenbeispiel demonstriert 
werden. 


Beispiel: 740: 329 = ? (nach dem Aufbauverfahren). 


1. 329 im EW an den Stellen 1 bis 3 einstellen, Wagen an die 6. Stelle, UW gleich- 
laufig. 740 im HW mit mdglichst kleinem Uberschuf8 aufbauen: 740:00325. Das 
UW zeigt 2°24925 = q, 

q, = 2°24924 (= erster Teilquotient). 


Bei diesem Rechnungsgang Bestimmung des Stellenwertes von g,. Das HW zeigt als 
UberschuB die Ziffernfolge wu, = 325. 


2. HW, UW loschen, UW aut gegenlaufig umschalten, Wagen an die 6. Stelle, 
im HW die Ziffernfolge 325 mit méglichst kleinem Uberschu8 aufbauen®: 325000018. 
Das UW zeigt 

do = Yo = 012158 (= zweiter Teilquotient). 


Das HW zeigt als UberschuB die Ziffernfolge wu, = 18 (= zweiter Rest!). 


3. HW, UW loschen, UW wieder auf gleichlaufig schalten. Der Wagen darf jetzt 
nur an die 5. Stelle kommen, weil die Eins von 18 an die 8. Stelle gehort. Im HW 18 
mit méglichst kleinem UberschuB aufbauen: 18000248. Das UW zeigt q, = 054712, 


dz = 054711 (= dritter Teilquotient). 
5 CoHatz: Z. angew. Math. Mechan. 16, 59/60 (1936). 
6 Da der UberschuB in diesem Beispiel stets um sechs Stellen weiter links im HW aufzu- 
bauen ist, mu8 die Ziffernfolge 325 an die Stellen 7, 8, 9 zu stehen kommen. Wiurde man sich 
verleiten lassen, den nach der ersten Umdrehung an den Stellen 6 bis 8 erscheinenden Wert 329 
durch subtraktive Umdrehungen in der folgenden (5.) Wagenstellung an 325 heranzubringen, 
dann bliebe die héchste Stelle des UW unentwickelt und man bekime daher nur 5 brauchbare 
Ziffern des Teilquotienten, was zu Irrtitmern fiihren kénnte. 
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4. HW, UW loschen, UW auf gegenlaufig umschalten, Wagen an die 6. Stelle, 
im HW 248 mit méglichst kleinem Uberschu8 aufbauen: 248000200. UW zeigt 


V4 = Va = 246200 (= vierter Teilquotient) usw. 
Ergebnis: 740: 329 = 2°24924 012158 054711 246200 


i, 2. 3: 4. Teilquotient. 
(Hingegangen am 13. Oktober 1953.) 


Bergman’s Linear Integral Operator Method in the Theory 
of Compressible Fluid Flow. 


B. Super- and Transonic Flow. 
By M. Z. E. Krzywoblocki, University of Illinois. 


Summary. In this second part of his work on Bergman’s linear integral operator method 
the author discusses supersonic and transonic flow. At first, the operator obtained by the use 
of Riemann’s function is thoroughly discussed as well as the question of the domain of convergence. 
In the section on transonic flow the application of the integral operators of the first and second 
kind in the case of the exact and the simplified compressibility equation is shown. 


Zusammenfassung. In diesem zweiten Teil seines Werkes bespricht der Verfasser die 
Uberschallstrémung und die Transschallstromung, und zwar auf Grund der Linear-!Integral- 
Operator-Methode von Bergman. Zunachst wird der Operator genau erértert, welcher sich aus 
der Verwendung der Riemannschen Funktion ergibt, und dann wird auch die Frage des Konvergenz- 
bereiches behandelt. In dem Abschnitt tiber die Transschallstr6mung wird die Anwendung des 
Integraloperators der ersten und zweiten Art im Falle der genauen und der vereinfachten 
Kompressibilitatsgleichung dargelegt. 


Résumé. Dans la seconde partie de son ceuvre l’auteur discute le courant super-sonique 
et le courant trans-sonique selon la méthode de Bergman concernant |’opérateur linéaire-intégral. 
D’abord V’auteur traite en détail le probléme de l’opérateur obtenu & l’aide de la fonction Riemann, 
ainsi que la question du domaine de la convergence. Dans la section qui a pour sujet le courant 
trans-sonique, il démontre l’application de Popérateur intégral du premier et du deuxiéme genre 
dans le cas de l’équation de compressibilité exacte et simplifiée. 


Part II. 
Supersonic Flow. 
1. Some remarks on different types of equations. 


Let us discuss briefly various simple cases of partial differential equations and 
the general solutions of these equations. The following three equations will be 
considered 


1 5 = é 
7G (Pua + Poo) = Fy, = 9, 4=U+Vs, Z3=u—-719, {3 5054) 
Pog — Vyy,=90, o 4% ,0=0, pe etyvaz ey (3.1.2) 
Woe +4(1—pw)¥Y,,—2%, =0 (3.1.3) 


In equation (3.1.2) by the transformation 2 + y = u,x — y= 9, Y (2, NS Y (u, 9), 
equation ¥,, — ¥,,—0, changes to 4%,» = 0. The first equation is of elliptic 
type, the second of hyperbolic and the third of mixed type. In the first case, the 
general solution is given by 


¥ = f(s) +9@), (3.1.4) 
and in the second one, by 


Y = f(u) +9 (9), (3.1.5) 
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where / and g are arbitrary (sufficiently many times differentiable) functions of one 
variable. As « and # are real variables, it is seen that in (3.1.4) there is an arbitrary 
function of one complex variable and in equation (3.1.5) two arbitrary functions 
of one real variable. The former remark follows, since in order to obtain real 


solutions, g must be chosen the conjugate to f; that is g=f(u—i9). By a 


transformation A, = (1 — yz)", equation (3.1.3) can be reduced to the form 
Was +-F.2, = 0 and the general solution is 
P=f(td—wa):—d8) +96 —p)e +8). (3.1.6) 


2. Differential equation in the supersonic flow. 
Introduce in (1.1.4) the variables log g and 6 instead of g and 6, i.e, Y= 
= 7 Piogg and ®, = q-!D,.,. Then we obtain: 
Op 0 Mee ~“Ptogg = — (1 — M) 9? Yo. (3.2.1) 


Differentiating the first equation (3.2.1) with respect to log q and the second one 
with respect to 6, and equating leads to 


ea els M*) Poo Se Fog a)iog 4 = 0. (3.2.2) 
For convenience, an auxiliary variable H = H (q) is introduced, given by 
2 q 
xed or H = \ od (logq). (3.2.3) 


0 y ‘ 


Since @/é (log g) = o 0/0H , we obtain from (3.2.2) the expression 


We now introduce the quantity A [compare with (1.2.1)], defined by 
dafdH = 0-1 (1 — M2, or da/dg = q-1 (1 — M2)’. (32 25) 


Then since M? is given by (1.9.19), we can find the value of 2 by a purely formal 
computation. In the case of a subsonic flow, we obtain the equation (1.3.13), and 
in the case of a supersonic flow, the expression 


4 (M) = + é{aretan [(M? — 1)'/2] — h-l arctan [h (M2? — 1)/:]}. (3.2.6) 


If we introduce the new variable A defined by 


A=ti, (o..2 7) 
we get 7 
A = h-' arctan [h (M2 — 1)"2] — arctan [(M? — 1)*/:]. (3728) 
In this case equation (3.2.4) becomes 
Van Pog + 4h. Fa = 9, (3.2.9) 
y, = 229 yee — 1h. (3.2.9a) 


3. Solution of the differential equation. 


The solution of the differential equation (3.2.9) follows an analogous pattern 
to that one applied in the case of a subsonic flow. As has been indicated by Berg- 
man}, there exist various operators for transforming solutions of one equation into 


1 §. Bergman: On two-dimensional supersonic flows. N.A,C.A., T.N. No, 1875 (1949). 
Ingenieur-Archiv VII, 4. 23 
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solutions of another equation. Often for different purposes it is convenient to use 
different operators or different forms of the same operator. In the present section, 
two operators are derived. They transform solutions of the relatively simple hyper- 
bolic equation 

Pia — Poo = 92% (3.3.1) 
i.e. a pair of differentiable functions g (3) and h (7), 3—=A+0,7=A— 8, of one 
variable 3 or 7, into solutions of the compressibility equations. As soon as the 
function WY (A, 6), satisfying equation (3.2.9) is obtained, the corresponding flow 
in the physical plane can be obtained exactly in the same manner as in the subsonic 
case. In contrast to the subsonic case, flows which in the logarithmic plane satisfy 
equation (3.3.1) have not been studied. Relatively little is known concerning how 
to choose g and h in order to obtain in the physical plane a flow around a prescribed 
obstacle, although there is reason to expect that this investigation should be simpler 
than in the subsonic case. On the other hand, it seems that this method of attack 
would not be of considerable interest in applications since most flows in which we 
are interested are mixed or supersonic with shocks. These combinations make necessary 
considerable modifications in the approach. 


4. Operator obtained by the use of Riemann’s function. 


One method of generating the stream-function of a supersonic flow from 
expressions g and fh [solutions of (3.3.1)] is by the use of Riemann’s function. To 
this end, let us transform equation (3.2.9). If in (3.2.9) instead of A and 6, the 
variables 3 and 7, defined in section III.3, are introduced, and, instead of WY, the 
“reduced stream-function”’ 


au 
pe — Wexp - ie ar} (3.4.1) 


a 


(where a is an arbitrary constant) is considered, then a formal computation shows 
that Y* satisfies the equation 


Poet ¥1Q 7 4) O° =o (3.4.2) 
dN 
ea N.2 Ua ee 
(k + 1) [5(k+1) 12k i“ ne ier 
5(k yt 
mare. Y tag 2 BO UF ope toe OT 8) ee ee 


Obviously it is allowable to take for the lower limit @ of the integral in equation (3.4.1) 

an arbitrary point, since replacing a by a*, a* + a, means only that ¥Y* is multiplied 
a* 

by a constant, namely exp | N,(t,)dt,]. The so-called “Riemann Integration 
4 i 

method’’? essentially permits the derivation of an integration formula which expresses 


the solution of the hyperbolic type of equation in a clear way by means of the initial 
data. Consider the most general form of a linear, hyperbolic differential equation 
of the second order in two independent variables: 


Gey tau, +bu, cu =0. (3.4.4) 


Suppose that we know that a solution of (3.4.4) exists. The “characteristic curves” 
of the equation (3.4.4) are the lines x = const., and y = const. Assume two 


is R. Courant and D. Hilbert: Methoden der Mathematischen Physik, Vol. I and II. Berlin: 
Julius Springer, Berlin (1931 and 1937). 
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arbitrary characteristic lines x — const. and two y = const. The so-obtained 
parallelogram has four vertices: A and B on one diagonal and P (3, 7), D (a, Yo) 
on another. Assume that we know the initial value of the function u along the 
“initial curve’ ADB. By means of the theory of characteristics we can determine 
also the value of uw, along AD and the value of u, along BD. This is the so-called 
“characteristic initial value problem”. 

The function v, depending upon the arguments x, y and the parameters 3, 1, 
satisfying the following conditions 


Vay — (UV), — (bv), +cv=9), (3.4.5a) 

Vz (2X, H; 3, 4) = b (x, y) v(x, 7; 3, 4) on AP, (3.4.5b) 
O, (4, 9; 3. 0) = 4,4) 0G, 9; 3,4) on BP, (3.4.5¢) 
D2: 2) ==) Ob (PF), (3.4.5d) 


is called the “Riemann function” belonging to the differential equation (3.4.4). 
The function w satisfying (3.4.4) is then given by the expression 


A B 
u(x, y) = u(D)v(D) + \ v (uy + au)dy + \v(u, + 6u) da. (3.4.6) 
D D 
Expression (3.4.6) represents the classical Riemann formula. Consequently, by 
means of the initial value on the curve 4 DB we are able to determine the function 
u (x,y). Let us apply Riemann’s formula to (3.4.2), i.e., let R (3, 7; 4, 7) be the 
Riemann function of equation (3.4.2). Then 


1 


Y* (3, 4) = Y* (30, No) FB (3, 73 30> 70) + | PF (40 Y) B (30: ¥3 30> Yo) dy + 
"No 
3 
= \ aes (2, No) R (x, Nos 30> o) dx. (3.4.7) 
30 
Since for ¥* (39, 7) we may write an arbitrary prescribed constant A and for Y* (2, 7) 
and ¥*(3,,y¥) we may write two arbitrary continuous functions g (x) and h (y) with 
the condition that g (39) = A (jo) = A, the formula 


rr (3, 7) = A R (3; UE 30> No) ae \ J x (2) R (x; "o> 30> No) dx + 


a0 


LU] 
+ \ hy (y) B (30; Y5 30. No) Cy, (3.4.8) 
no 
represents an integral operator which transforms two arbitrary functions g and h 
of a real variable into a solution of equation (3.4.2). Thus according to (3.4.8) the 
values of ¥* (3, 7) are determined within a rectangle having vertices at 39, 1/9; 3, "0; 
400 733. %, by prescribing its values along the line segments: y = 1%, 39 S% S 3, 
aNOLa =o, Yq = Gory. ‘It 1s possible to show that the Riemann function of 
equation (3.4.2) can be written in the form 


R (3, 43 30 %o) = ae ee at (3: 73 30> %o)> (3.4.9) 
0 
where 
R (3, 13 30 No) = 1, (3.4.10a) 
37 : 
R® (3,3 30 70) = \ | Fi (% + y) dx dy, (3.4.10b) 


30 "0 
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5” 

R® (3,73 30) No) = \\ Fi (a + y) BR (2x, Y3 $0. No) dx dy, (3.4.10¢) 
Zo %o . 
an 

R (3,73 30,0) = \) Fi (a + y) RO» (x,y; Bo: No) dx dy. (3.4.10d) 
30 "0 


The tabulation of R = R (3, 3 30, %9) is rather complicated since R is a function 
of four variables. Below, for the illustrative purpose, the second term R® of the 
expansion (3.4.9) is computed. 


5. Evaluation of the function R® (3,7; 4, 7%). 


In this section, the second term R” (3,7; 39,7) in the series for the Riemann 


function will be evaluated. The function F, is a function of the variable A = ava 
and therefore 
4 
| Py (a + y) da =I (5 + y) —T (39 + y), (3.5.1) 
where : 
b 
VP, (x) dt =F (4) — T (), (3.5. 1a) 
to 
and 
ieee ” 
\dy|\ Fy (x + y)dx =\ (P(g + y) —T (39 +y) dy = 
No do "o 
= FO (3 +4) —F(3 + jo) — FP (3p +0) +P (Go +10), (8-5-2) 
where 
ty 
| £® (x) dr = T® (t,) — T (ty). (3.5.2a) 


to 
For A =0, i.e, M=1 [see (3.2.8)], the quantities N,(M) = N,(1), (3.2.9a) 
F, (B) = F, (0), (3.4.3) or F, (A) = F, (0) tend to infinity and this fact causes a 
certain amount of inconvenience in the tabulation of the functions / and J'@), 
According to (3.4.3), for k = 1°4, 


F, (B) = 0°45 B-* + 0-63 B+ — 0-21 B2 — 0°51 — 012 BY, (3.5.3) 
and ’ ; ; 
A = |/6 arctan|(|/6)~! B| — arctan B. (3.5.4) 


The values of M, B, F,, [' (2A), [@ (2A) are tabulated in1, where an ordinary 
numerical procedure was applied. 


6. Second type of integral operator. 


In analogy to the subsonic case, where formula (1.9.2) yields a representation 
for solutions of equation (1.9.1) in terms of an arbitrary analytic function of one 
complex variable, there is derived in this section a representation for solutions of 
equation (3.2.9) in terms of two arbitrary differentiable functions of a real variable. 
Obviously, instead of equation (8.2.9), equation (3.4.2) for the reduced stream- 
function Y%* may be considered. The desired representation can be obtained as a 
consequence of the following theorem: 

Let H” (rt), t= 2A, denote functions defined by the following recurrence 
relations 


EY (t) = — F, (2), or EO @ = — | FiGp adr (3.6.1) 
0 
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EBOtY (7) — — EO (¢ ) — F, (t) E™ (x), (3.6. 1a) 
or . 
Hee+0 (x) = — | [BY (1) + Fy (4) BE (x,)] dey, (3.6.1) 
0 ’ : 
where F(t) is the given function which must satisfy the inequality 
2 2 
Riv a for’'0 7 = 4, (3.6.2) 
Ito T| 


where uw and tT)» are conveniently chosen constants, and the derivatives of F, must 
satisfy the inequalities 


d” F, (ct) — 2m + I)! 4 tp ph Nuts 

ROD Ee eg ae Ey, ee ave 
e., for t real and positive. Then the series 
VE 8s 9) =F (3) eB a) fO1G3), (3-6-3) 
1 
and 
V® (3,9) =g (n) +>) B™ (3 + 7) g™ (x), (3.6.4) 
1 


where / and g are two arbitrary differentiable functions of a real variable and 


pret = | fo (j) 431, f(g) =F, 


(3.6.5) 


ge* (ny) = \ 9 (9g) dn, gH) = 9, 


HPs Oop 


represent the solutions of (3.4.2). The proof of this theorem is carried out under 
the assumption that the series (3.6.3) and (3.6.4), as well as those for their first 
derivatives, converge uniformly. This assumption permits one to interchange the 
operations of summation and differentiation. A formal computation yields 


Vad, a (Bf) + BO f{), (3.6.6) 
or with the use of (3.6.5), ie., ps = f@—, we obtain 
VD — f+ f BO 4 d7 fm (E” + Be+), (3.6.6a) 
il 
The second differentiation (with respect to 7) gives 
Vin = ft ED i fm (BS + Bet), (3.6.7) 
Also BS 
Pe VOU sf Hao: he ie Ee. (3.6.8) 
1 


Substituting (3.6.6a) and (3.6.8) into (3.4.2) yields: 


yy noe vo — {EP + F,| oy, {|B oe ee a Bo), (3.6.9) 


Since += 2A = 3-47, we have EM (2A) = H(z) = BOE eke ehus 


the first bracket on the right-hand aide of (3.6.9) vanishes due to (3. 6.1) and the 
second one vanishes due to (3.6. 1a). 
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7. The proof of the uniform convergence. 


In section III.6, the assumption was made that the series (3.6.3) and (3.6.4), 
as well as their first derivatives, converge uniformly. Of course, it is sufficient to 
prove that the series and both their formal (i. e., term-by-term) derivatives converge 
uniformly. Assume that in a certain interval J there hold for certain functions A (t) 


and A (rt), and all their derivatives the inequalities 


d"A (zx) 
aan 


|A (x)| < A(n), ey, (3.7.1) 


for t contained in the interval J (which symbolically will be denoted by 1 « J), then, 
as is already known from the section I., the function A (rt) is a dominant of A (rt), 
which fact is symbolized by writing A S A, tel. If E™ (zx) is given by 


T 


EO (2) =\Fy (tq) dry, Fy (t) > Fi (0), (3.7.2) 
0 
then 
|Z (x)| 5 BY (x), ve I, (3.723) 
and also 
dB (x) | dB (t) : 
| dee aaa = = rare (3.7.4) 
Thus’ £™ (zt) is a dominant of HE (r), i.e., EY S E®, Suppose now, that 
B+) (x) = | HY, () + Fa (ty) BO (dey + H+” (2), (3.7.5) 


0 
where E™ (t)< E(t), 0<H™ (rt). Then it follows immediately from (3.6.1b) 
and (3.7.5) that 

|E@+) (7)| = H@+D (q), (3.7.6) 


| dB “- ih) 


Geer 
dt ? 


dt 


Is 
and by considering the corresponding derivatives of [He +F, E”| in comparison 
with |B + F, B| it follows that 

E™ +) (7) << B+ (7), Cie a 


Thus starting with the assumption that #™< H™, we have proved by induction 


that B" +) <H@+. Now by expression (3.6.2) we may assume that F, is dominated 
by the right-hand side of (3.6.2), i.e., 


: ee a = 
SpA aks Sapsyicese |: = F, (t). (3.7.8) 
Suppose now that H™ (zr) is given by 
a c 
H®” = —-> 0, (3.7.9) 


where the ¢,’s are some conveniently chosen positive constants, which will be 
determined later. Then, an explicit expression may be obtained for E,1 and namely 
with HO —1 and c, = 2mut (3.7.9) we have from (3.7.5) 


EY — \Fy (v1) dt, + HY = 2 22, (to) — F)- (377,10) 
( 
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Assume that 
EO) (t) = ¢y (tog — t)7*, (3.7.10a) 


and in order to obtain the recurrence relation holding between the c,,, write with 
the use of (3.7.9) 


t 


as 2 r ; — 

Be+ (=| AS Til) GO MD ahah OR ec age 7 Pe — 7) ty 
0 i (to9 73)" * 2 ' (To9 — 4)" * 2 seit Gna (Foo Z 

(3.75.10) 

Cnn 1 = Cy [M 12 ty (mn + -1)-1]. (arid 2) 


Relation (3.7.12) supplies the quantity c,—= 27>, if the assumption is made that 
€>= 1. Consider large values of x, which permit us to write 


2. too (% + 1)-? <1, (3.7<15) 


or 
Oe a5 oe Gg. hte ol) se, (0-4-1) (3.7.14) 
Sel (8 1) mlm sed prea ». Sregihi2e3n ome): 
Therefore we see that for every value of ~ we can put down 
Ces Bea) (357%, 15) 
where * =>c, is a conveniently chosen constant. 

The function 7 (3) is assumed to be differentiable and finite, which means that 
for all values of 3 less than a certain value 3,, say, the values of the function f (3) are 
less than a certain constant value y., say, which can be written symbolically in the 
form 


lf (@l <#2. for [3] Se. (3.7.16) 
From the expression (3.6.5) in conjunction with (3.7.16) we obtain 
42 3 
| Sed, PS, Spay --- (3.7.17) 
or in general: 
i (g)| SA for: ah Ste (3.7.17a) 
Consequently, substituting (3.7.10a), (3.7.15) and (3.7.17) into (3.6.3) we obtain 
; 21 u* |3| (n + 1)! y* |5" 
7) (3 ee clea! bars ts 
V (3; m<t{1 1! (t¢9—T) ets l Rit (re ee) I (3 7 1 ) 


The series inside the square bracket itself has a dominant, namely the series 
* n 
: GegkB Ua Nal tu. dae Ach At lal bf (3.8.18') 


which converges if 


Zook. & bao ee ey (3.7.19) 


The series inside the square bracket of equation (3.7.18) converges if 


A+ 6 


Se eee for Pl NC ee AA (3.7.20) 


Thus we proved that the series V has a dominant which converges and has all 
terms positive. According to the fundamental concepts of analysis, the series V 


is uniformly convergent. 
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8. The domain of convergence. 


The series V“ converges in the domain given by (3.7.20). Assume that the 
quantity t)) is negative: tT), <0. Then the domain of convergence is bounded by 


: 3 
the following two straight lines in the t (= 2 A), 0-system of coordinates: 31 + 
+0—%.=9, et 6 —t»=9. Similarly, we may prove that the series V®, 


given by (3.6.4), converges uniformly in the domain 


al I I fe | 
Too — 2 | ( 
or wees 
Fr a Dati ate (3.8. 1a) 
|T>o— TI 


which is bounded by two straight lines: ; t— 7 — t= 0,5 t 4 0 — to9 = 


The intersection of both domains determines the domain in which both series converge. 
The upper limit on the positive part of the (2 A)-axis is given by the quantity 
T) = 2 Ag, (i. e., vertical line). In a particular case we may assume Tt), = 0, which 
implies that the domain of convergence is a triangle bounded by the following three 
lines: + 6 = 0;rand + = t,. 

The last question which should be mentioned is the investigation of the behavior 
of the function F, (2 A). That is, we should show whether and where its derivatives 
satisfy the inequalities (3.6.2a). This question is too complicated to be discussed 
here, and the reader who is more advanced mathematically is referred to the original 
literature!, where proof is presented that the derivatives of F, satisfy the in- 
equalities (3.6.2a). 


9. Solution regular at d=r1r=0. 


The solution (3.6.3) and (3.6.4) and the theorem presented in section III. 6, 
cannot be applied directly to equation (3.4.2), as F, given by (3.4.3) has a pole 
for A = 0. That is, if M varies from 1 to oo, the quantity B= (M? — 1)’: varies 


from 0 to oo, and A, given by (3.2.8) varies from 0 to = (A-1 — 1)z. Thus for 


A = 0, the quantity B= 0 and the function Ff, —> co. In some cases, however, it 
is possible to overcome this difficulty by shifting the origin. Let a be a positive 


number a and 1 1 
=F =e eA, Rea Gn. (3.9.1) 


so that A* will mean A* = A — e t* — 71 —a. Equation (3.4.2) then assumes 


the form | Wee + Fo (3*, 7*) W = 0, (3.9.2) 
. 1 i 
where’ W (3*,7*) = Y* (27 = 7% y* +754), F,(t*) =F, (t* +a), so that F, 
is analytic for t* = 0. The corresponding changes in HH will be indicated by 
writing H” (rt) thus: 
BO (c*) = — | F, (t,*) dry* = — | F, (x,) dr, = B, (2), (3.9.3) 
0 a a 

Ben (ge) = — | [BO (ry*) + Pe (0) BO (2,*)|dry* = 

: ° (3.9.3a) 
= | [ee (t,) + F, (t) B(t,)\dr, = Hoty (t). 


a“ 
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The domain of convergence, remaining the same, referred to the t* system, is merely 
a shift of the original domain. It may be desirable to have explicit formulas for the #™ 


and their derivatives, as these may be employed in computation; if » is given the 
values n = 1, 2,3,... in (3.9.3a), then 


B® (rt) = —| F(t) dry, (3.9.4) 
HO (t) = Fy (t) — Fy (a) + 1/2 [B® (x), (3.9.5) 

E® (t) = F, (t) B® (t) — Fy, (t) + Fix (a) — 
i (3.9.6) 

—|\ FY (t) de, — P, (a) BY (x) + 1/6 [LE (1), 

where " 

BO (a) = F(a), (3.9.7) 
ky (t) = Fy. (t) — F, (t) B® (x), ete. (3.9.8) 


The expressions for H“ and E® (x), ES (x), are cited in?. Attention must be called 
to the fact that in those formulas instead of A, the quantity rt should be substituted. In 
A 


the same report Bergman cites the expressions for H® (A), F, , (A) and \ Pe (A) dA. 


Assume for the quantity in (3.9.1) the value a = (h-! — 1) a, so that 
P* = '2 A¥ 2 Aha — Ix (3.9.9) 


This means that the function F, (t*) has a singularity for A = 0 or r* = 2 A* = 
= — (A+ — 1)a, corresponding to an extreme case M -—> oo. This transformation 
may be very convenient in some cases and a more thorough discussion of it may be 
found in!. Let us mention that because the quantity A = 1/2 (kh — 1) a represents 
an upper bound for A, it is obvious that for the upper bound of the variable +, the 
quantity t) = (h-1 — 1)a may be chosen. That quantity t, should be used in the 
considerations presented in section III.6. It is obvious that PF, (t*) = PF, (2 A*) = 
= F, [2A — x (kh — 1)] = F, (2 A) is an analytic function of A which is certainly 


f “ 1 
regular in the circle of radius -, (kh — 1). 


10. Transformation to.the physical plane. 


In order to carry out the transition to the physical plane, the values of Y,* and 
Y,* are needed. Without going into too many details which may be found in the 


original literature®, we shall present here the final results: with H given by (3.2.3): 


— BA 1 : 
w*— BH (q 1) | ot ee it abs +3M,+M,| (3.10.1) 
Y,* = H[M, — My), (Be1029) 


where V =>’ W® (3*, 7*), 7 = 1,2, represents the solution of (3.9.2) and the 


values of M oe are given in the original papers of Bergman. (See the section discussing 
the tables available.) If the integration is varried out along a streamline, ¥* = const., 


3 §. Bergman: On supersonic and partially supersonic flows. N. A.C. A., T. N. No. 1096 
(1946). 
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then the corresponding values of « and y, that is, the image in the physical plane 
of the eae flow, are given by the formulas* 
= | (9 g Yy*)— [B? Y,*? — ¢ ¥,*?] cos 6 dq, (3.10.3) 
= (og? Ya") ee Ua a a ee ag (3.10.4) 


11. Alternative formulas. 


Sometimes the range of variability of the speed is comparatively small. In these 
instances, it is useful to replace F, by a constant, say F',. Equation (3.4.2) then 


becomes 
PR Pt = 0, Cage 


and its solution can be written in the form 


il 
, | 26 (3m) (1—#?) if (Se) dt 
eB: 1) =\{cos| Py Ile 2 | lt, 2 i} a—eyle’ 
1 


: (3.11.2) 


where f and g are two arbitrary, twice-continuously differentiable functions of one 
variable. When applying formula (3.11.2), the independent variables of functions / 


and g should be substituted by the expressions E 2fe r)| and 
Es Cie fa geal )|, respectively. 

In many instances, it is necessary to have representations of the solutions which 
are valid in the whole strip — 0 <98<o, 0<2A<2a(h+— 1), say. One 
possibility of obtaining such formulas consists in approximating F, (2 4’), 
2A’ =2A—(h+— 1)a, by a sequence of functions F,™ (2 A’), each of which 
satisfies the inequalities (3.6.2) and (3.6.2a), and which are chosen in such a manner 
that 


lim FY” (2 A’) = F, (2 4’), (Si1ice) 

m—> Co 
in the interval 0 <2 A’ <2(h1—1)a. Let Y*™ (32’, 7’) denote the solutions of 
Yee ee) (3) ey) PY is iene ae (3.11.4) 
where 3’ and 7’ are taken in the interval — 3 < 3’< 3, — 7% <y’'<7, and 
3 and 7 are values, corresponding to the quantity rt = 2 (kh — 1) a, equivalent 


to t) = 2 Ay in the section III.7. The quantity + is the upper bound for the 
variable +’ = 2 A’ < 1 = 2(h-!— 1)a. The functions Y*™ are also to satisfy 
the conditions: 


pein) (37 0) = P* (3/, 0), W*™ (0, 0!) = W* (0, 17’). (3.11.5) 
Then, in a manner similar to that in Part I, it can easily be shown that 
Lina, WF 8) (37 95?) eet ae (3.11.6) 
m—> CO 


in the given interval. A sequence of functions fF,” (2 A’) can be obtained in the 
following manner: Let 


Fy (2 A’) = om + oA, (SL ioey 
be the series development of /, around the point 2 A’ = 0, i. 2A=(h3—1)a 


This series obviously converges in the circle of radius A eS 1 l)a, ics 
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it is relatively easy to show that the series (3.6.3) and (3.6.4), obtained as solutions 
of equation (3.11.4), with F,™ given by equation (see Fig. 1) 


(m) (2 A’) = Wee en _ Atn 
He (2A) = I'(1+n/m) ae, eB: 
n=0 j 


will converge for every finite positive number m < co, in the strip. — co <7 '< 60; 
0 <2A <a(h*— 1). On the other hand, according to a classical theorem of the 
theory of functions, the relation (3.11.3) holds in this case, so that for the functions 
Y* (3', 7’) obtained in the foregoing manner, the relation (3.11.6) holds in the 
whole interval — 3 <3’ <3, — 7 <7! <7, 

Below, we shall give another form of the solution of equation (3.4.2). Suppose 
that F’,, is a function which possesses a continuous first derivative. Let 2,* (3, 4,4) 
and F£,* (3, 7, t) be solutions of 


(1 = i?) dea — {1 Sa + 2t3 LETT + F,, 3 a | = 0, (SV LESS) 
and 
(1 — #) B*,,,—t B*,, 42 ty [%s), + Fy Het} = 0) “ (3-19-10) 


respectively. Let #, and #, possess continuous second derivatives, and_ let 
H*,, (nt) and EH*,, (3¢)+ be finite for t—0. Then 


+1 


* (3,9) = | {Ei* (1 Oh ly38—8)| + BY G0.) fly 28) 
=5 


—— 


(1 — #)- de, 
(QalTsih) 
where /,, # = 1,2, are two arbitrary, twice-continuously differentiable functions 
of their respective arguments, is a solution of the equation 
pe, + F,, P* = 0. (3.11.12) 


The proof of this theorem is given in‘. Let F,,, (6) possess derivatives of all orders 
in the interval 6, <6 < fy, 0 < By < By < oo. If a constant c exists such that the 
inequalities 


ree ele J 

| oo er KAGE yay) (err 13) 

can be obtained, then there exist solutions #, and H, of (3.11.9) and (3.11.10), 

respectively, satisfying the conditions of the theorem cited above. Thus, we obtain 

a solution for ¥Y* in terms of two arbitrary, twice-differentiable functions. In® 

Bergman shows that (1.9.2) and (3.11.11) are different forms of the same operator, 

the former valid in the subsonic case, while the latter holds in the supersonic range. 

Similarly, as in the subsonic range, if Y, (A, 0), » = 1, 2, 3,..., represents a set 

of particular solutions of equation (3.2.9), and A,, arbitrary constants, any linear 
combination 


N 
EEA TAO ay: (3.11.14) 

v=] 
is also a solution of equation (3.2.9). This is so since the equations are linear and 
therefore the principle of superposition of solutions holds. By varying the constants A,, 
y = 1,2,..., flows around various shapes can be obtained. On the other hand, it 
is often necessary to determine constants A, to yield a flow which approximates 


4 §. Bergman: The approximation of function satisfying a linear partial differential equation. 
Duke Mathematical Journal, Vol. 6 (1940), pp. 537—561. 

5 §. Bergman: Methods for determination and computation of flow patterns of a compressible 
fluid. N: A.C. A., T.N. No. 1018 (1946). 
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that about a prescribed boundary curve, the equation of which is, say, F (x, y) = 0. 
In this case, the problem of how to determine the constants A, is somewhat 
complicated and unsolved up to the present. Some remarks on this subject are in?. 


Part IV. 
Transonie Flow. 
1. General remarks. 


The problem of transonic flow is perhaps the most difficult. We have to find 
solutions which must partly cover the subsonic, and partly the supersonic regions. 
The complexity of the problem is the origin of several methods of solutions. Below, 
we shall try to outline and discuss a few of them. In the second half of this part a 
more general discussion of the problem of transonic flow will be presented without 
deriving the proofs which a reader who is more advanced mathematically may find 
in Bergman’s original papers. 

First we shall present a continuation of the subsonic flow solution through the 
sonic line, according to the proposition of von Mises’. 


2. Behavior of the G, at m = 1. 
The equation (1.5.3a) provides a recursion formula for the determination of 
the G,. It may be written in the form 


A 
d ° 
Qui = oe +\fG@,da,  Gy=1. (4.2.1) 


If 'we compute this last integral, we obtain an integration constant each time, since 
the lower limit is arbitrary, i.e., G, would depend on » constants of integration. 
By means of reasoning identical with that in section I.9, we choose as the lower limit 
of the integrals the value 4 = — oo, corresponding to f= m= 0 (1.3.11) and 
(1.3.13). However, at 4 =0 or m= 1, all G, become infinite. It is therefore 
desirable to introduce other functions r,, in the place of the G,, which remain bounded 
(i. e., finite). To find such 7,, we have to study the behavior of the G, near m = 1, 
or, what is the same, the behavior of f (7) and G,, (7) in a neighborhood of 7 = 0. 
Assuming that 


Gat TOL OT ee Oe te ae ee (4.2.2): 
we introduce into the recursion formula for the G,, the function r, (7') defined by 
PUTS SoS ee am | (4.2.3), 
Since Gy=— 1, we have ¢, = 1>7,= 1. From (4.2.2). we obtain: 
dG 
ai wd | ne lg (4.2.4). 


where 7,’ denotes the derivative of r,, with respect to 7. According to (1.3.13) we 
have, after a few transformations, 


da 
ar =~ 27 [(k +1)— 2k T+ ke — 1) F3. (4.2.5) 
By multiplying the right-hand side of (1.3.14a) and (4.2.5), we obtain 
da 1 
B= fap = ay TTA TH) oe) 2G) ae ay 
é 
=—F(k+1)T+[l—2k{5(k + 1j}4T...], (4.2.0) 


Sev: von Mises and M. Schiffer: On Bergman’s integration method in two-dimensional 
compressible fluid flow. Advances in Applied Mechanics, Vol. I. New York: Academic Press 
Ine. (1948). 
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From 
dG reat dG, 2 
Grir= qe =\fG.04 = GP A+ | BG, af, (4.2.7) 
where sere : 
aT I 
A= = — 5 [+ 1) — 2k 7? 4+(k—1) TM) 7, (4.2.7a) 
we obtain, by using (4.2.3) and (4.2.6): 
De eC eb es Fe Oat 8) Sab gt eye 
fe 


+26, 7+...) [(b +1) — 2k + (k—1)™) —Z& +) e,( 241 — 
1 


=e Wee Lyra ae (1 4b, oe.) on dae, (4.2.8) 


whence 
Cn+1 E1035 T—(& 
ees ( +3) + bay (3 +,1) As | 


=T-69+9(5-3n(k +1) +4 


(3n 4 3) (k + 1)|+ (4.2.9) 


+ T-en4 tee 4 [5 Gn 2 Gu 1)]|6,—h| 30 +e Bn + alls. 


Comparing the factors of 7 @”** on both sides of (4.2.9), we obtain 


Sth = 5 8n(k+1) +35 (b+1)(8m +3) = 


=53n (k + 1){1 +5 [86 (nm + 1)]-}. (4.2.10) 
a recursion formula for the c, that yields finally 
ae E 3 (ke + yl’ nt vat +5 (36y (y+ 1)]-4¢. (4.2.11) 
Since c, = 1, we derive further from (4.2.10) | 
= srk t+), | Po (de 


which determines all c,. Comparing now the coefficients of 7~@”"*” and using 
(4.2.10), we find 


Onai{l +5 [86 n(n + 1)] JI) = by (1 — 28m) + 520 n— I) (ee 
— (2k(k —1)2 +k [6 n (3m + 1) (k +.1)]4, (4.2.13) 


which together with b, = 1 determines all 6, and justifies our particular choice of 
the functions r, (7). If n is large, we obtain from (4.2.10) 


1 


Cn pie ez n(k +4 ), (4.2.14) 

and . 
b, @% —'2nk(k +.1)>, tim paey eae (4.2.15) 
Taking k = 1:4, we obtain from (4.2.13): Go OL On 4, Op Fee 


b, = — 2°8,; etc. 


4=— 
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3. Recursion formula for the functions 7r,. 


By substituting (4.2.2) into (4.2.7) we obtain: 


. 
Opeth UL) Oe aCe Seen Sn T-Ger Dy) AG, ar rae) dx. 
1 
: (4.3.1) 
4 


Representing the function A (4.2.7a) by a polynomial 4 = 7-2 > oT? and 
v=0 


the function B (4.2.6) by a power series B= T~-4)’8,7", and substituting 
& PNY 
them into (4.3.1), we obtain 
4 fe [e-e) 
Cyt Op Pace = DO TUL Ty S80, ) ae TOP ON ee ee ee ee ery 
y= 0 e y=0 
In order to compute the r, (7) from this recursion formula, we have to keep in mind 
that r, = 1 and that we obtain from (4.2.10) for k = 1-4, the values c,/cy) = 1/2, 
C,/¢, = 77/20, etc. Then we obtain 
oe co T 
1 


en(=T |e) sede = 0*\ Bada (4.3.3) 


Using this result for the next recursive step, we find 


4 fis 
Fer(L) = Sa, 7 (Pry — 37.) + T¢(| Bar), (4.3.4) 
y=0 1 
where the second term on the right-hand side originates from the expression 
ay T x 
1 | oo? Bi(a)r; (edi TOS BiG Bihar de (4.3.4a) 
1 it 1 


by means of integrating by parts. The results, (4.3.3) and (4.3.4), have one 
remarkable feature. Since B (7'), taken in its original form (4.2.6) without developing 
into a power series, contains only even powers of 7’, it appears that r, (7) and r, (7) 
are free from logarithmic terms and can be developed into power series at 7’ = 0. 
In general, this will not be true. The power series for r, (7), when introduced into 
the integral in (4.3.2) in order to determine r; (7'), will produce a logarithmic term. 
This term, however, is multiplied by T°, so that r,(7’) has a regular development 
at the point 7’ = 0 up to the ninth order. Further application of the recursion formula 
will preserve this logarithmic term and will also produce new additional terms of 
the type 7” (log 7)“. It is easily verified that we have in general 


n— 2 
rn (T) = 3) P (P) log TY, n> 2, (4.3.5) 
y= 0 

where the P\”(7') are regular series of 7’. One may easily deduce recursion formulas 
for the P® by means of (4.3.2). Here, we wish to point out only that each P™ (7) 
with » > 0 has just those powers of 7’ which are larger than 8. We shall restrict 
ourselves to determining the first coefficients of PY? (7). Let us assume the 

development 
PD SD ee Oa (4.3.6) 


yes 
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From (4.3.2), we find for the coefficients 9 the following recursion formula as 
long as y $8, » <3n+3: 


ner Oy On? = S{(H— 3m) xp +B,.l — (30 +3)}-4} 0. (4.8.7) 
u=0 
Certain numerical values of the. coefficients o” are given in’, p. 269. 
Let us note that it is possible to apply another expansion, since in the neighborhood 
of T = 0 (or M = 1), T and TT can be developed in a series of ¢ = [— 3(& + 1) (A/2)] 
or ¢ = (— 3°6 4)* for k = 1:4. This choice yields T =e +... If these series are 
substituted in (1.3.14a) for 7 and 7-1, a Dette of f in the power series of «, 
which holds for M <1, is obtained in the neighborhood of M = 1. 


4. Methods of analytic continuation. 


In the region (M < M, <1), where M, is near 1, the series (1.9.2) converges 
very slowly, and it is therefore necessary to employ a large number of terms in order 
to obtain a good approximation for Y*. When applying equation (1.7.7), the 
number m must be chosen rather large. If this is the case, it is then expedient to 
replace the expansion (1.9.2) by (1.9.9). Theoretically, this is, however, not the 
only way of overcoming this difficulty, and in the following other means of doing 
so will be indicated: this alternative approach employs the method of analytic 
continuation. 

Let Y (q, 6) be determined in a domain, say H, and let ¥,, n =1,2,..., be a 
complete system of particular solutions, each WY, being determined in a domain G. 
Suppose that zi and G actually do overlap, and denote their common part by I. 


Further, let » a, Y, be the series expansion of Y in J. Frequently, a hike es 


n=1 n=1 
will converge outside of J, say in the domain H, — J, where H, is G or some part of 
co 
it. If, in addition, »’a,¥%, can be termwise differentiated twice in H,, it represents 


n=1 


the analytic continuation of WY in H, — I. ed UeMLy the domain H, in which 


the stream function can be represented in the form vy a, Y, covers a supersonic 
=1 
region as well, and consequently this method will then yield the flow in this latter 


region. In this manner, a method for determining a mixed flow may be obtained. 
Below, we shall present some such methods and at first let us cite an auxiliary 
lemma: 

Let pq, 9), @ S9 Sm, —L S90 SL, be an analytic function of two real 
variables, g and 6, and let 


» | (q) cos (77>) + b, (q) sin (77")], (4.4.1) 
v=0 
a, (9) = +L" iA (q, 0) dd, (4.4. 1a) 
oS =f 
“ 0 
a, (q) = L\ p (qs 8) cos ("7 ) a0, (4.4.1b) 
i aid sy ; 
: 0 
b, (q) = L-*\ p (q, 6) sin (7) dd, (4.4.10) 


is 
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be its Fourier development. The series (4.4.1) converges uniformly and can be 
differentiated termwise any finite number of times with respect to both q and 6. The 
proof of this lemma can be found in any book on Fourier series, and is also cited in’. 
Since every solution of an elliptic equation with analytic coefficients is an analytic 
function of two real variables, the result obtained can be applied to the case where 
p(q, 9) is the stream function ¥ (q, 6) of a subsonic flow. Thus 


‘pp [a (q) cos ( zee + b, (q) sin (23")|, (4.4.2) 
ey , 
I 
iy (q) = -y L-1\ ¥ (q, 8) 48, (4.4. 2a) 
; —L 
a, (q) = L\ ¥ (q, 8) cos (777) dd, (4.4.2b) 
ry pes: 1h ey ; 
ke 
b. (q) = L-1\ ¥(g, 8) sin ("7") a8, (4.4.20) 


—L 


can be differentiated termwise. If, now, following Chaplygin, we introduce instead 
of g, the variable 
t= @ (24,2), (4.4.3) 


the then equation (1.1.5) for Y assumes the form 
(27 (l— rt)? YW] +h — (264+ 1] [270 = 7) 20 —1) PY =0, “4G 


where 6B = (k — 1)-1. Using the representation of VY (qg, 6) by means of the Fourier 
series, given above, and differentiating termwise, gives the result: 


> {|e (ro at ae) [1 — (2684+ 1)t]- (ex(l—7)] 1(1—1)-*- 
» 2B | 008 as +| 5. (<a 18 | [l= (28 + 1) siete ee 


watb,|]_. vw8 | 
47 | ay vee (4.4.5) 


J -( cA Wek 


In this equation, the coefficients of cos as and sin (73), for .+ ==1, Zee 


must be equal to zero separately, and therefore the a, and 6, are each solutions of 
equations 


da, yo oe? 
Afra ra al zt {1 (2p 1) r} {oe (Ql =r) (hb Sz) Waa ae 
(4.4.6) 
db, »? a? 
mae ler {1— (28 4 1) t} {e(1 — A}: (1 — 7) nO, 
(4.4.7) 


Each of these equations can easily be transformed into a hypergeometiic equation 
and thus every solution of them may be written in the form 


(A, F, + B,F,*) eP/2@D, (4.4.8) 


7 §. Bergman: On two-dimensional flows of compressible fluids. N. A. C. A., T. N. No. 972 
(1945). ah Soe 
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where A, and B, are constants, and 


F, = F (%, B,;.— B; 1 4), (4.4.9) 
He eV city Oe" Bay xm cttorectts 5, Bos Sorts Boil ait); (4.4.10) 
F (x, B, y,t) being thé hypergeometric series. The symbols used are given by: 
1 
yw = (vL-1 +1), o = =[(z — 8) +4] (4.4.11) 


Lf» 2 Yo 
p=> (= 6) A, A, = ( sr) (2B +1) + | 
In order to determine the constants A,, B,, the following theorem is employed: Let 


PY (9,0) and Ye (7,9), [@ S¢ 5%, —L S80 SL] be solutions of an equation 
of elliptic type. If, along a line, say ¢ = q, 


PO (qo, 6) = P© (qo, 9), | (4.4.12) 
and 

Pde Meier Ve tree, (4.4. 12a) 

then in the whole domain [q, <q <q, —L <6 <1), 
POG 50) = eX GsG): (4.4.13) 
Now, suppose that the function g (Z), (1.9.2a), is regular in some domain H, + H,, 
Hi due. 9y y le 0S Dy) (4.4. 14a) 
Eps [Oh 3 Grn sl Se], (4.4.14b) 


which domains lie in [6? <3 2?, 4 <0]. Then, by the formula (1.9.2), it follows 
that WY is also regular in H, + H,. Suppose that Y has been evaluated in the 
domain H,, but it is desired to avoid the evaluation of Y in H, by means of (1.9.2), 
since this series converges very slowly in H,. To achieve that, represent the function Y 
in H, by a Fourier series, which is always possible since we evaluated VY in H;,: 


Ww) as | (q) cos (=7~| + b, (q) sin =") (4.4.15) 


v=0 
Assume that in H, we wish to represent Y by another Fourier series 


co 


wr?) ay [| se” (rz) cos ( zue + S® (zr) sin (75*)|, (4.4.16) 
yv=( 
SO) SAPP, Be Fee ie 1,2. (4.4. 16a) 


In order that this series will be the solution Y under consideration, the constants 
A® and B® must be determined so that 


iG Cah te eee yO CP) Oy geod ea (4.4.17a) 
S@ (z) d= = O» (Yo), te 0, 1, 2, PN) (4.4.17b) 
dS¥» (z) da, (4) | LAST 
dq g=a dq |\g=a’ (oe Lae) 
dS (x) db, (q) 
ar =- (4.4.17) 
dq |\g=a dq \g=u 


Since in [4 <¢ SQ, —L<9SL), W(g,6) is an analytic function, the 
series (4.4.16) and its derivatives converge uniformly and absolutely in this domain, 
and this series represent the solution Y under consideration in the region H,. 


Ingenieur-Archiv VII, 4. 24 
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Moreover, this series and its derivatives may also converge outside of H,, say in 
H,=([¢.<¢9 <4 —L<6<L). If H, lies partially outside of the domain 
[a2 <3 22, A <0], then the expression obtained gives the analytic continuation of 
the solution outside the domain of representation by the integral formulas (1.9.2). 
In particular, H,; may include some region which lies in M> 1. Of course, in 
practical application the summation will include only a finite number of terms N, 
say. Thus the procedure of continuation consists of determining a linear combination 
| A, y (M, 6), so that along a line M = M,< 1, this latter function and its 
derivative with respect to M coincide with the corresponding quantities for the stream 
function of the subsonic flow pattern which is to be continued. 

Another similar method is proposed by Bergman in ” and 8. 

Unfortunately, the power series for the function of (4.4.6) and (4.4.7) converge 
so slowly that in most cases this approach must be abandoned for computational 
purposes. On the other hand, if certain changes in this approach are made, numerical 
results may be obtained. Two such procedures will be discussed in the following®: 

(a) The solutions of (4.4.6) and (4.4.7) are expanded in the form of a power series, 
not, however, around t = 0, which is the case of the hypergeometric series, but 
rather around some conveniently chosen value of t which lies inside the interval 
under consideration (— 2 <0 L). For example, it will often be convenient to 
introduce for t the variable o = t — 0°15, say. Since the particular solutions are 
often considered only for a small range of variation of t, say 0°13 <1 < 0°20, the 
coefficients of (4.4.6) and (4.4.7) may be approximated by polynomials in o, so 
that if, say, L = 2/2, the equation approximating (4.4.6) and (4.4.7) is (for 
B= (kh — 1)71)=.5/2); 


(— o° a 0°5500 o? + 0°2325 6 — 0°0191) - + (1°5000 o2 ae 1°4500 o ate 0°1838) oo a ieee 
— (6:000 ¢ + 071000) % U = 0, (4.4.18) 


where U denotes a function approximating the actual solution of (4.4.6) and (4.4.7). 
As is known, every solution of (4.4.18) can be represented as a linear combination 
of two independent solutions, each of which can be written in the form of an infinite 
series in o, which at least converges for |o| < 0°15, a range of convergence which 
will be sufficient for most purposes. Solutions U, (oc; v), U,(c; 7) are chosen so that 


UU (05.9) = OF Oe 0a) a, (4.4.19a) 

U(0>>).= 1, ) U0" 9p) = 0. (4.4.19b) 

Tone Wiss a etG! se ecteee eee ae (4.4.20) 
U,(o;y)=1+b,0°4+6,08e4+..., (4.4.21) 

a;, b;,..., = 2,3,..., being functions of ». Substituting (4.4.20) and (4.4.21) 


into (4.4.18), we obtain a set of equations which give the values of a; and b; as 
polynomials in »? (3, p%*). Since it is often sufficient to consider the series (4.4.20) 
and (4.4.21) in a small interval, say |o| < 0°05, it will suffice to employ only a 
few terms of these series, and hence it is necessary to compute only a small number 
Gf G5, U5: 


(b) To employ the second method mentioned, equations (4.4.6) and (4.4.7) are 
approximated by the equations 


: du du Ol bie 
O15 _| "85 ae “OE . id 2D) 
(0'15 + 0) (0°85 — 0) $4 + (1-225 +150) 4 — » (Sisse) ¥ = 0, (4.4.22) 
BS Bergman: The hodograph method in the theory of compressible fluids. Supplement 
to ,,Fluid Dynamico‘ by R.v. Mises und K. O. Friedrichs, Brown University (1941—1942). 
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and two independent solutions of this equation, wu, (c; »), Uz, (a; v), are determined 
not by the power series (4.4.20) and (4.4.21), but by an approximation method. 
That is, we write 


d 
oe or Ag v Ao. (4.4.23) 


Then 


(1°225 + 1°5 o) (0°1 — 60) 
dv =| Far + o) mae eay | | | oe +o)? res ey |e) Ao. (4.4.24) 
Assume that Ao is contained in a small interval, say 0°001, and that for o = 0, we 


have u=0 and v= 1. Then, if (4.4.23) and (4.4.24) are employed, Au (oc) and 
Av (ao) may be determined at o = 0, and hence 


u (0°001) = uw (0) + Au (0), (4.4.25a) 
v (0°001) = v (0) + Av (0). (4.4.25b) 
From these values of u and v, Au (0°001) and Av (0-001), and hence wu (0-002), v (0°002) 
may be found and the entire procedure continued until an approximate curve is 
found for uw (c), 0 So < 0°15. In this manner, the desired integration is performed. 


Both methods described above are appropriate for the purpose of computation, but 


they are insufficient to give an insight into the behavior of the particular solutions 
when y — oo. 


5. Application of integral operators to transonic flow. Basic equations. 


Below, the application of integral operators to transonic flow will be discussed. 
In general, the possible method of continuation, which will be discussed, will consist 
of the determination of the limit values of the subsonic stream function and of its 
derivatives with respect to M on the line M = 1, and then of the consideration of 
the initial value problem for the compressibility equation in the supersonic portion 
of the flow, M > 1. Since it is advisable to use different variables in the region of 
transonic flow, we have to begin with the basic equations. Assuming that the thermo- 
dynamical conditions follow the isentropic law, and introducing as new variables 6 
and 


g d 1 
1 ey 
H=\¢ (dgiq) =\q* [1-5 &—Ve|®™ ag, (4.5.1) 
an Nn 
we obtain the following linear equation for the stream function: 
where 
1 tid 
M=q|l (ke @| 2 (4.5.3) 


It is noted that in considering purely subsonic flow in previous sections, the 
variable 2 (1.2.1) was used; and for purely supersonic flows, the variable A = 1 A. 
In studying transonic flow patterns, it is found to be more convenient to use the 
variables H, 6. A formal computation shows that the Taylor development of / (7) 
in the neighborhood of H = 0 is 


L(A) = [2/(k — 1)J@ H/F —Y{(— 2 A) — [(2k + 5)/(2& + 2))- 
ALO geet A0 a hho ea AY « E> acre (4.5.4) 


In purely subsonic flow, we applied formula (1.9.1) in previous sections. Introducing 
24 
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the complex variables Z = A + 76, PB EID, equation (1.9.1) transforms into 


4W-+4N (W,+ %) =0, (4.5.5) 
where the following transformation formulas may be helpful: 
0 Pope ag 0 1 [a ne 
oT ee (4 : aa) ea (+ t iz), (4.5.5a) 
e 0 7) é a fane é 
Boe. 1 ger Wee lz 7 = (4.5.5b) 


It should be noted that the interval — co <A <0 corresponds to the interval 
— oo <H <0. Inthe supersonic case, we applied equation (3.2.9). Equations (1.9.1) 
and (3.2.9) can be transformed. If WY is replaced by 


we — WRA, eR[o =N, (4.5.6) 
then Y* satisfies the equations 
Pe Page eee =O) and Pia ee eae (4.5.7) 


where 
peng ene +1) M1 (26 =) eek stat Me sae (4.5.8) 
The expansions of N and F in the neighborhood of 4 = 0 are(®, p. 861): 
N =(1/124) {1 ee "| (2 +52) 5 Fee 2| 
» [2-Me 3% (he - 1)-e(— ay] +... a (4.5.9) 
Bre 8/(80(222 4 MG. (4.5.10) 


In the vicinity of H = 0, | (H) may be replaced by the first term in its expansion. 
Using this value of / (#) in (4.5.2), we obtain the so-called “simplified’’ compressibility 
equation: 

— CH Yoo + Pay = 0, (4.5.11) 


where C is a constant. In considering transonic flows, the solutions of the simplified 
equation will give a fair approximation of the exact stream function. The expressions 
for N, F and H will, in this case, reduce to 


N= NYS Hay (4.5.12a) 
Pt PY = 5/6 (l (2245 (4.5.12b) 
A = H+ =] (35/2) (2k S11) 2B 2) ay. (4.5.12¢) 


In his important investigation, Tricomi!® studied the boundary value problem 
of equation (4.5.11) and showed that if we consider a finite domain located partly 
in the subsonic and partly in the supersonic region, then under certain conditions 
the boundary value problem has a unique solution. Frank!" considered questions 
allied with Tricomi’s investigations in the case of the exact compressibility equation. . 


® §. Bergman: Two-dimensional transonic flow patterns. Amer. J. Math. 70, pp. 856—891 
(1948): 

10 F. Tricomi: Sulle equazioni lineari alle derivato parziali di 2° ordine di tipo misto. Atti 
R. Accad. Naz. Linzei, Mem. della Cl. Sci. Fis. Mat. e. Nat., Ser. 5, Vol. 14, 133—247 (1923). 

1 F. I. Frankl: On the problems of Chaplygin for mixed sub- and supersonic flows. Bull. 
Acad. Sci. URSS, 8, 195—245 (1944) (in Russian). — On the problems of Chaplygin for mixed 
sub- and supersonic flows. Bull. Acad. Sci. URSS. Vol. 9, 121—143 (1945) (in Russian). Also: 
N. A.C. A., T.N. No. 1155 (1947). 
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We shall discuss below the integral operators (of two kinds). The integral operator 
of the first kind may be applied to the subsonic flow. It is useful to shift the origin 
to the point 2, 9) (point of reference of the operator). Every real solution of equation 
(4.5.5) can be represented in a neighborhood of the origin Z’ = 0 as the imaginary 
part of the right-hand side of (4.5.13) — 


pig) = R22) 9 (2) + aS goon Tey (2,2) 9 (2')], (4.5.13) 


L'(n + 1) 
with cox 
Vip ey Lai 3=Z' 
In (Z a \ \ Z ae fare 
Gea) i= \ \9( n)AZ,, dZ, = (n—1)! (Zac2g) 9 (3) 43, (4.5.13a) 
0 0 0 3=0 
where 
+1 
/ I / 9\1 
g(Z')=\f| 52 qa #)| def(. — 2 (4.5.13b) 
eS il 


with a suitably chosen function f (or g). The function g (Z’) is called the associate 
of p (g) with regard to the integral operator of the first kind, p. The function R is 
defined by the formula 


id Z’ oe 
R(Z',Z') = exp|— | N(Z +.Z,')dZ,'|, (4.5.14) 
0 i 


and the symbols Z’ and Z’ denote: 
Oe a eth a. GA = Aya =O — Oy, Ag 0. (45 uo) 
Another form of the operator of the first kind is 


ok ie l 
Pi) = |B (2.2,0 [52 —#)| ay — ey (4.5.16) 
k 


where f is a suitably chosen associate function, called the cotati function of the 
operator P,, and E, is equal to 


EA 7 t= RAL 7) B22 .t), (4.5.17) 
where H,* has the development 
E,*=1+)> 2s (2" P (Z’, a (4.5.18) 


The functions P™ satisfy the fellewins recurrence relations: 


(1) fe + 1) (n) 


PZ +2F=0, (2n +1)P ej ee p+ 2FP =0, Wass 1 28h ae 
(4.5.19) 
with 
BET OO... pes Notes 


By the above requirement, the P” and hence the generating functions H, (of the 
first kind) are uniquely determined. All the proofs—for example, the proof that the 
series (4.5.18) converges absolutely and uniformly in a sufficiently small neighborhood 
of the origin Z' — 0, Z’ = 0—can be found in®. Thus, assuming that the associate 
function is regular in a sufficiently large domain, we may prove that Be applying to 
it the integral operator of the first kind, we obtain a solution of (4.5.5) defined in 
a sufficiently small neighborhood of the origin. As the next aa one can 
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show that if f is regular in the entire domain M <1, this solution can be 
continued throughout the whole subsonic region, i.e., the representation holds for 
E[M <1, — «<9 < oo), where the symbol H#[..] denotes a domain. There 
exists still another form of the intetral operator of the first kind, namely 

Zz 


p(g)= R(Z’,Z'){g9,(Z') —\ \ F9,dZ,dZ, + 


00 
VALE NSB IA, 
saad i tu aA dZ,dZ,+...}. (4.5.20) 
00 0 0 


Making some modifications, we can extend the definition of the operator P, so that 
it can be applied to functions of one variable A + 6 and A — 6, respectively, thus 
generating solutions of equation (3.2.9) in the almost entire supersonic region. The 


expression 

Priel a 8)] + Ps Gh 30) ], (4.5.21) 
where /, and f, are two linearly independent functions, will represent a (possible) 
stream function of a supersonic flow pattern, i.e., for H[M >1, — co <6 < ool. 


The integral operator of the first kind, convenient though it is for many purposes, 
has the disadvantage that it does not represent solutions of the compressibility 
equation in the neighborhood of the sonic line. Furthermore, it has the disadvantage 
from the practical point of view that the P™ (Z’, Z') are functions of two variables Z’ 
and Z’, or 4 and 6, which makes tabulation of the values of P very time-consuming. 


6. Integral operator of the second kind in the case of the simplified 
compressibility equation. 

The two disadvantages of the operator of the first kind can be removed by the 
use of another operator—to be termed “operator of the second kind’’—for which 
the P™ are functions of one variable only, and which yields a representation of the 
stream function in the neighborhood of the sonic line. This operator has a number 
of other distinctive features, which will best be elucidated by the detailed discussion 
of the so-called “‘simplified’’ compressibility equation, i. e., where N = N+ = (12 4) 
in equation (1.9.1), or, alternatively, where F = F'+ = 5/(144 4?) in equation (4.5.8). 
At first, let us represent the generating function in the case of the simplified 
compressibility equation. It is of the form 


vag ees (4.6.1) 

where H is given by (1.9.7); it can be given by the formula 
H (2:4) = 8S, (— 2.4) [CE 2 ay], (4.6.2) 

with . 

p [(— 2 aye] = 14+ 8,(— 2ae + 8,(— 2Aseh+..., (4.6.3) 
So= 9(2k + 1)/(6k —6) 3—"Ye (k a Hea a (4.6.3a) 
S, = (1/10) (3/4) (& + 1)-“s (2% + 5), (4.6.3b) 
S, = — (1/1400) (3/4)! (k + 1)—*2 (6442 + 70k + 75), ete. (4. a 


In the case of the simplified equation, we have N+ = (12 A)-, p [(— 2 A)s] = 
H+ (24) =S8,)(— 24)-*, F+ = 5/(144 22). One can show that in this case, ie 
generating function H can be given in the form 
E+ (Z',Z',t) = Et (A, 0, t) = A, S, (— 2 a)-* F 1/6, 5/6, 1/2, X) + 
+ B, Sy (— 24)-** [— # (A + 2 6) JieF (2/8, 4/3, 3/2, X 
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for |X| = |— @ (A +7 4)/(— 2 a)| <1, (4.6.4a) 
and 
E+(Z',Z’, t) = E+ (A, 0,t) = A, 8, [(— #) (A + 4 0)—"* F(1/6, 2/3, 1/3, X) + 
+ BS, (— 2 Ay“ [(— #) (A + 4 6)]-* F (5/6, 4/3, 5/3, X), 
for |X| = |— 2a/— #@ (A +7 6)| <1. (4.6.4b) 


One can prove in this case of the simplified equation that the imaginary part 
of P,* (f), 


+1 
Py () =| B*(4,0,) [541 — ey fae — ey, (4.6.5) 
si! ee = 
=A+76, (4.6.5a) 
is a solution of (1.9.1) with N = N+—1/(12A4), ie., 
Wee nn bPse iy}: (4.6.6) 


The solution obtained is valid in the large, i.e., in a finite domain. In particular, 
if the function 

eae ae 

g(Z)= | f|y2a—elaya — ey, (4.6.7) 

t=—1' 
is regular in a region B which possesses in its interior a number of branch points, 
each of finite order, then the obtained solution (4.6.6) is defined in B and possesses 
branch points at the same points and of the same order as (4.6.7). We assume here 
that the only singularities of g in B are branch points. In applying the integral 
operator method in the case where g has poles or logarithmic singularities, certain 
modifications are needed. If the function (4.6.7) is defined for all values 4 < 0, 
— co < @ < oo, then the solution Y will be defined in the same region. 

The expressions (4.6.4) obtained by means of the explained procedure, will not 
necessarily be analytical continuations of each other. Since, however, the hyper- 
geometric equation has only two linearly independent solutions, the constants, A,, 
B,, and A,, B, can always be adjusted as to make these two solutions analytical 
continuations of each other. 

The integral representation (4.6.5) can immediately be generalized to the 
supersonic case where it will produce, in an analogous manner, solutions of (3.2.9) 
with N, = N+ (4.5.12a). Indeed, replacing 4 by the variable 2 A, it may be verified 
relatively easily that W (i A, 6) satisfies equation (3.2.9) with N, = 1/(12 A). 
Repeating the procedure which led to the generating function (4.6.4) in the subsonic 
case, we now obtain the generating function 


E+ =[a, S,/(2 A)!*] F (1/6, 5/6, 1/2, X) + [b, Sp (A + 6)"2t/(2 A)"*] F (2/8, 4/3, 3/2, X), 
for) x) = |? (LE-— 6)/(2 A) =< 1. (4.6.8) 


As is known, a hypergeometric series in z defines a function which is analytic when 
|z| <1. This function has a branch point at z = 1 and if a cut (i. e., an impassable 
barrier) is made from + 1 to + oo along the real axis, the function is analytic and 
one-valued throughout the cut plane. In our case, where |X| is substituted for z, 
if F («,, By, y1, X) denotes only the hypergeometric series and not the hypergeometric 
function, (4.6.8) has to be replaced for |X| = |#@(A + 6)/(2 A)| > 1, by 


E+ = [a, Sof(@ (A + 9))"*] F (1/6, 2/8, 1/3, X) + 
+ [by So (2 A)/(@ (A + 4))""] F (5/6, 4/3, 5/3, X), (4.6.9) 


where the constants a,,b, are expressable in terms of a, }. 
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If 6 =A and #=1, there arise certain difficulties, since the hypergeometric 
functions in (4.6.8) will then become singular. By the transformation formulas 
of the hypergeometric function, (4.6.8) may be written in the neighborhood 
of [(A + 6) (2 A)+]# = 1, in the form 


E+ = [ay So[(2 Ay] F (1/6, 5/6, 3/2, X) + 
+ b, Se [2 AK2 A — #(A + 6)" F (1/3, — 1/8, 1/2, X), (4.6.10) 


with now 


XS. Pas Hie. 


In order to avoid the complications which arise from the fact that the second term 
of (4.6.10) has a singularity for A = 0, # = 1, we shall therefore take b, = 0. ‘The 
function H+ will accordingly be of the form 


E+ = [az So/(2 A)'*] F (1/6, 5/6, 3/2, X). (4.6.11) 


We note further that all these considerations can be repeated with A replaced by 
— A. Our operator will therefore yield two independent types of solutions, depending 
on whether the argument of the associate function is taken as 4 -++ 0, or A — 6. The 
exact conditions under which our operator can generate solutions of the compressibility 
equation in the supersonic case are given in the following theorem: 

Suppose /,(8), s = 1, 2, are real functions of the real variable and everywhere 
differentiable with the possible exception of 8 = 0; suppose further that in a fixed 
neighborhood of 8 = 0, f, can be approximated to any prescribed degree of accuracy 
by the expressions of the form 


AVS ; 
SACD OR a, lanl Ki, Sean (4.6.12) 
n=1 
then 
Y (A, 6) = Ry (fs). Re: (fe), (4.6.13) 


+1 
i ' 1 1 
Re (f,) =| Bt (4,—(— 1 8,0 f(g (4—(— 8) (1— 8) dey —Ays, (4.6.14) 
i=—1 

Sil, 
represents a solution of the compressibility equation, which is defined for any A > 0 
and can be interpreted as a stream function of a (possible) supersonic flow pattern. 
The proof of this theorem is given in *. In (4.6.14), the function H+ is defined by 

(4.6.8) or (4.6.9) or (4.6.11), depending upon the value of |X). 
Before’we proceed to investigate the behavior of our solutions on the sonic line, 
we note that (4.6.4a) is not the only solution which depends on one variable. The 
generating function (4.6.4) for the operator of the second kind may be replaced by 


E+ (A, 6, t) = A, Sy (— 2 d)-* F (1/6, 5/6, 1/2, X), (4.6.15) 
XS (Ad, 4 Oe Ay, 


Besides its greater generality, the generating function (4.6.15) (with J, + 0) has a 
number of additional features, which make it superior to (4.6.4) in many cases. 
In (4.6.4), the point 2 = 0, 6 = 0, is a singularity since by approaching this point 
in a suitable manner, the argument X of the hypergeometric function can be given 
an arbitrary value. In (4.6.15), such a singularity does not exist, since 4 and 
A— i, +176 cannot vanish simultaneously if 2, +0. Another singularity, which 
can be removed by using (4.6.15) with J) + 0, occurs in the supersonic case. For 
? = 1 and A = 6, the second term of the right-hand side of (4.6.10) becomes singular. 
In order to allow for this case, we had to assume that b, = 0, thus somewhat restricting 
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‘the generality of the solutions we could obtain. Setting 4 = 7A, it is seen that in 
the supersonic case the argument of the hypergeometric function in (4.6.15) becomes 
X=f(A+7A, + 6)/(2 A). (4.6.16) 
For t= + 1, we have X + 1 for real values of A, 6; the singularity in question is 
therefore removed and the constant b, in (4.6.10) may now take any arbitrary value. 
Thus, by the use of integral operators of the second kind with the generating 
functions (4.6.4), (4.6.8), (4.6.9), as well as (4.6.15), we obtain solutions 
YW = I m P3 (f) which are defined in the subsonic and supersonic regions respectively. 
Let us repeat that if f (4.6.7) is defined in a domain G, in the subsonic region and 
has, as its only singularities in G,, branch points of finite order (but not poles or 
logarithmic singularities), then the generated function is again defined in G,. In 
particular, it has branch points at the same points and of the same order as g. If g 
is a twice-differentiable function of one real variable in a domain G, in the supersonic 
region, the generated function will be a solution of (4.5.5). Thus, by this procedure 
we can generate solutions of (4.5.2) which are defined in adjacent domains, one in 
the subsonic, the other in the supersonic region. Our object is to show that these 
solutions can be continued across the sonic line. This fact is seen immediately if 
we introduce a new variable s 


Se) “for A= 0; 


a 4.6.17 
s= — As, for A> 0. ( 

We note that in the simplified case we have: 
§ = — 2-3-7 [2k — 1)J@-H/BE—| A, (4.6.18) 


The variable s = s (M+) (where M~ denotes the fact that the simplified case is taken 
into account) considered as a function of Mach number M~*, possesses the property 
that s(1) = 0, and that 


ds (M+)/dM+ = — 2 (3k + 3)~*s M+ + 0(1 — Mt), (4.6.19) 


is non-vanishing and bounded (finite) in a sufficiently small neighborhood of M+ = 1. 
In (4.6.19), the symbol 0 (1 — M+?) denotes the terms of order (1 — M*?). More 
details concerning this particular item may be found in the original paper®. Thus 
our solutions can be continued across the sonic line, since by expressing both of them 
(one on each side of the sonic line) in terms of the one variable s, we may easily show 
that the generated function Y (0, s) is an analytic function of 6 and s (and therefore 
of 6 and M+) at the point 0=0, s<=0. 


7. Integral operator of the second kind in the case of the “exact” 
compressibility equation. 
As was shown in Part I, in the subsonic region, equation (4.5.5) has a solution 


in the form 
Y —Im([P (f), (4.7.1) 


where 
Pif)=\ E4201 (52 (1— 0) aa — ey (4.7.2) 
The generating function of the integral operator P (f), which produces solutions of 


this equation, is : 
EB =H (2A) E*, (4.7.3) 
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and the function H* can be represented by a series 


co 


E*=14 3) (@Z)Q™ (2'A), (4.7.4) 
n=1 
where the Qs are solutions of the system: 
(2%. + 1). © +07 424 FO = @ (4.7.5) 
A 
QY = — 4 | Fda, (4.7.5.8) 


[for relation between A, A and M see (1.3.13) and (3.2.8)] the function F’ being 
given by equation (4.5.8). The above series converges for |Z|< 2 |A|, or for 
62 < 3 42 and was thoroughly discussed in Part I. It is therefore desirable to obtain 
solutions which are defined in |Z| > 2 |A| or 62> 3/7. Similarly, in the supersonic 
region a solution was presented which is convergent for — 34<6< A, if only 
the series V“ is taken into account, or for — A< 0< A, if the sum of the series 
V® + V® ig considered (see Section III.8). It is desirable, therefore, to obtain 
solutions which are defined for — 3 A > 6 and for 6 > A. At first, we shall discuss 
the subsonic region. We shall present the following theorem without the proof which 
the reader may find in °: 


Let 
oe n— > + As) +9) 
(0? =S On weg ETE, ya, (4.7.6) 
v=() 


be a set of functions which are connected by the relations 


Qi? +4Fq” = 0, (4.7.7a) 
2 (1 = ox Gi? Ae gts rene Nags ae ae egy, (4.7.7b) 
Te ers 6 cae lae 
Then each of the functions 
Be) — S79 /(— eZ)" etch (4.7.8) 
n=0 


multiplied by H is a generating function of the integral operator P (f) and each of 
these series converges in H# [2 |A| < |Z|]. 

In order to obtain a continuation of a given stream function to the supersonic 
region, replace A by the variable s = (— A)’s. Then the generating function in the 
subsonic case may be written as 


E® (4, 6, t) = H (2 2) E* (A, 6,4) = 


Co 69D ile SEN 
= 268, p(s) Sy Cee eng), Be (4.7.9) 
n=0 ¥=0 
for 7 = 18) 
[8] Sos, 2. 18 NZ i aoe te ee 
where 


Z = (—syh +76, (4.7.9a) 
and the domain of regularity is a circle of radius so, with the center at the origin. 
By replacing 2 by 1 A we see that (— A)’ is changed to — A’. Thus we obtain 
as the generating function in the supersonic case, 


co Cc 


, ee: : = 1 
HO (A, 6, t) = 2-6 By (a) —1S" YT #08" (— BIO 9] PZB CO 
n=) v=0 


(4.7.10) 
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for 
Silay, Del allem tT bong: peated 
where 
Z=i(—syh+i0. (4.7.10a) 
As mentioned above, the proof is given in °. 
Let us now restrict ourselves to a neighborhood of the point A= 0, 0= 6,4, 
% +0, lying entirely in the domain D = E [|s|< 4p, |3|<3-‘ |6["]. We may 
show that if f is regular in D, then both the generating function and the associate 
function, and therefore the expression P [f (— s‘2 +7 6)] in the subsonic case, may 
be expanded in integral powers of (@ — 6,) and s'!. Similarly, in the supersonic 
case, the expression P [f (i (— s)*2 +4 6)| may be expanded in integral powers of | 
(6 — 6.) and (— s)". One may also show that both functions P (in sub- and super- 
sonic regions) are determined and equal to each other for s = 0. In the next step, 
we may show that if / is regular in D, then the solution P [f (— s‘2 + 76)] is an 
analytic function of s and (@ — 6) at (0, 6,). Thus 


P if (— #2 +2 6)| =» f, (0 = 6) 8. CBee 
n=0 
Since s is unchanged by the result of putting 7 A for A, we obtain 
Pif i (—s +7 6)] =>» Ff, (0 — 9) 8”. (4.7.12) 
n=0 


The expressions P[f(— s":+76)] and P[f(¢(—s)/: +76)], equal functions of s 
and (# — §,), are analytic continuations of each other across the sonic line. Thus, 
assuming that the associate function is regular in a sufficiently large domain, and 
applying the integral operator of the second kind, we obtain solutions of a 
compressibility equation defined in four adjacent domains, namely 


D,=E(M <1, 6>32|A(M)|]+2[(M +1, 06 >A(M)], 
D,=E(M <1, 6 <3 |1(M)|], 

Da eM aol iO era MOE) EME Oe 8 AM), 
D,=E(M +1, —3A(M) <6 <A(M)). 


The solutions defined in D, and D, were presented in the present section, while those 
defined in D, and D, were derived in Parts I and IIJ, respectively. 

In the simplified case, using the theory of hypergeometric equations, it was 
possible to combine these representations into one, which yields solutions of (4.5.2) 
defined in the whole (M, — 6)-plane. In the exact compressibility equation, the 
problem remains of combining these four representations into one. This problem can be 
attacked by using the integral operator of the first kind in addition to that of the 
second kind, and, in analogy to the simplified case, developing a theory of differential 
equations with singular coefficients, which would furnish us with information 
corresponding to that used in the simplified case. More information concerning this 
item may be found in?. 

We should like to add here a remark regarding the general question of analytic 
continuation of a solution ¥ (Z, Z) of a linear partial differential equation. If W 
is given in two different domains, say B, and B,, by different representations, say 
in B, by the integral operator of the first kind in the form 


Y= WW, =p (lg (Z)1 +p lh (2), (4.7.13) 
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[see (4.5.20)], and in B, by another operator 
+1 
w— ¥,=|E(Z,Z,1)1 (52 qo P)) (dt/(. — 8"), (4.7.14) 
— 
[see (4.5.16)], (not necessarily of the first kind), and the origin lies in the inter- 
section of both domains B, and B,, which may be denoted mathematically by a 
symbol B,M B,, then the problem of the analytic continuation of ¥Y, into the 


co 
domain B, is equivalent to the determination of g and h from a given if =) a2 


n=0 
One may show” that setting Z = 0, and Z = 0, respectively, in the relation 
+1 
eS — 1 : 
plg(Z)] +p (h(Z)]=\BZZ,Hf(5 20 —#) aya — ey, (4.7.18) 
t=—1 
we obtain the identities 
co n —s 
D> 2 Se a, = 9(Z) + £0, Z) h(0), (4.7.16) 
n=0 v=0 
Xoo > Zr = R(0, Z)g (0) +h (4), (4.7.07) 
n=0 
where 1 
o@) = | BO), (1 — ey o> dt, (4.7.18) 
t=—1 
; 
72 = | BO () (1 — eyed, (4.7.19) 
t=—1 
#(Z, 0, t) = S’EO() 2", (4.7.20) 
n=0 
E (0, Z, t) =) B® (t) 2. (4720 


It is thus seen that the analytic continuation of Y, into B, and all the problems 
arising from it, such as the determination of the singularities, etc., are reduced to 
similar problems in the theory of functions of one complex variable which are given 
by ‘their power series expansions. 

One more item may be mentioned briefly, and this is the determination of the 
associate function in terms of the given stream function. There are possible cases 
in which the stream function is given in the form of a power series; say, and we want 
to have the associate function. In many instances, the stream function VY is given 
in some different form, say the values of Y and 0¥/éM are-given on a line M = const. 
If these quantities are analytic functions of 6, then from these data it is possible to 
determine the associate functions. Problems of this type are sometimes of importance. 
In *° one may find the procedure for the derivation of the formula expressing f in 
terms of the values Y = y, (0) and 0¥/@M = y,(6) on the sonic line, M = 1. 


8. Combination of the integral operator and Chaplygin’s solution. 


We now want to indicate another procedure which aims at constructing “mixed’’ 
flows. In applying this method, we use in addition to the integral operator, 


8S. Bergman: Certain classes of analytic functions of two real variables and their properties. 
Trans. Amer. Math. Soc. 57, Vol. 299—331 (1945). 
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Chaplygin’s solutions. Suppose that g (Z4) is a complex potential which in the 
(incompressible) physical plane yields a.flow past a closed curve. For the sake of 
simplicity, let us assume that the flow is symmetric, and the complex potential (in 
the logarithmic plane) is a two-valued function, possessing as its only singularity a 
branch point of the second order at the point 4, corresponding to the flow at 
infinity (see the section on singularities for more details on this subject). Employing 
the operator of the first kind, we obtain a solution of equation Y) = Im [p (g)| 
in the subsonic region. For the illustrative purposes, let us assume that this solution 
will not be determined at first in its entire domain of definition, but only in a portion, 
say D, of this domain, namely in the portion situated in [A < A,], 4) < 0, say, so 
that, in general, the boundary curve Y = 0 is not closed in D. Let us assume 
that A,, <A) holds. We shall now describe a procedure to determine the analytic 
continuation of the solution Y outside of D. Since in the following we have also 
to operate in the supersonic region, and as we shall employ Chaplygin’s solutions, 
it is convenient to introduce, instead of /, the variable 7 = (k — 1) q?/(2 a"), where 
a is the speed of sound at a stagnation point. Let t,, = 7 (A,,), and let us introduce 
Chaplygin’s solution: 


(A, F, +.B, F,*) c!@ exp [4 (a » 6)/L], (4.8.1) 

where 
Pell (op, hs = B, 1 =) (4.8.2) 
Yeas a (1 = Tet Bly, = Ny. Vy — iy 2 oe Lo zy, (4.8.3) 


F being the hypergeometric series. Here 
‘ 1 1 
»=($+)), p= (k —1)4, x ==|($—A) +4] (4.8.42) 


6 =s|(4 6) A,| 4,=|(s7) es+y+ ze”. (4.8.4b) 


We assumed, above, that we know the solution YW (7, 0) = Y® (A(z), 6) in D. 
In order to compute the continuation Y® of Y, obtained in D by the use of the 
integral operator, into the domain B,, say, we write WY) in the form: 


Sj / o 
wo — S"[(AY FP, + BOF) cD cos (77°) + (AMF, + BOF,*) v2 sin (27°), 


va 


(4.8.5) 

and determine the coefficients AM, B™, x = 1, 2, y= 0,1, 2,..., so that on t = 15, 
—L<6<L (where + L denote the lower and upper bounds of 8), 

PONT ,.0) == l (Ty, 0); (4.8.6) 

PY) (to, 0) = PO (to, 9), (4.8.7) 


holds. It can be shown that the series (4.8.5) converges in a certain strip, say 
[vt <1 <1,,--L, <6 <1]. By replacing 1, by t,, L by L,, and repeating the same 
procedure, we may continue the function to a range of values [t, = ti = Te, 
— L, <6 <L,], and so on. If we employ the integral operator of the first kind, 
under the assumption that the function g in the expression eee Im [p (g)] is 
regular in a sufficiently large domain, we can prove the existence of ¥ in the whole 
subsonic region, as well as the fact that by this procedure the function Y can be 
determined in the subsonic portion of the domain of definition of the flow. If the 
flow can be continued to the supersonic region, then the procedure described also 
yields the continuation to that region. If, on the other hand, we replace the infinite 
series (4.8.5) by a finite one (in doing this, we may require that 
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an (To, 6) SN (To, 9)| = © (4.8.8a) 
[0 (x9, 6) — PO (zy, 8)| Se, (4.8.8b) 


where ¢ is any preassigned quantity, however small), then the finite sum is defined 
everywhere, and if ¥ = 0 is a closed curve, we obtain in this way a “mixed” solution 
of the exact equation, which together with its normal derivative has a jump (of the 
value <<) on the line t = 1,. This jump, however, can be made arbitrarily small. 


9. Tables for the determination of transonic flow patterns. 
Alternative formulas. 


The preparation of tables for the determination of transonie flow patterns needs. 
some explanation, which will be given below. For a better understanding, let us. 
begin with an incompressible fluid. A plane is introduced whose Cartesian coordinates. 
are the variables 9 and 

qd 
n = \@oq-1 dq = log q, (4.9.1) 
i 
with 0) = 1 for an incompressible fluid. The complex potential for a large class of 
symmetric obstacles in these variables can be written in the form 


co 


g Zi =D) be in (4.9.2) 


v=0 
with z= 7 +76 and 7 equal to the logarithm of the speed at infinity. In the case of 
symmetric obstacles, the b,’s are purely imaginary. The representation (4.9.2) holds 
only in a circle with origin at 7 = 7) and of radius equal to the distance from 7 = 7p, 
to the nearest singularity. However, by employing suitable summation methods, e. g., 
9 @) = lim DY) 6/0 Cee = teh o (4.9.3) 


where the symbol J’ denotes the Gamma function, a representation is obtained for g (2) 
in a much larger domain (for details, see reference 86). With the constants b, purely 


imaginary, the stream function ¥Y can be written 


Y = Im(g(z)] =>) ib, B (H, 8: 1), (4.9.4) 
v=0 
yo) (750, 4) = Re lz = yah”. (4.9.4a) 


The streamlines of the flow patterns corresponding to g(z) in the physical plane are 
given by 


a =>) ib, X°) (ivOisiaa)s (4.9.5a) 
vy=0 
y =>) ib, Y™ (7, 8; 19), (4.9.5b) 


where 
AY (1, 0; No) ad \ Oo." e7l0 {[— pp” cos 9 — 


— Y sin 0] dy + [(H cos 6 — wi" sin 6] 0}, (4.9.6a) 
Vo. (7) Og) = \ Oo 2 e”% {{— La sin 6 + 
ry (”) cs mee (7) “ — (y 

+ Wi" cos O) diy + [YP sin 6 + Y” cos 6] dd}. — (4.9.6b) 
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The coefficients b, depend on the boundary curve of the profile. In employing the 


hodograph method in the incompressible fluid case, it is convenient to prepare a set 
of tables for 


p40 (7, 9; mo), x (7, 95 mo), yo (7, 9; mo); (4.9.7) 


for a number of values of the parameter 7. For, by using a set of b,’s appropriate 


~ 


to the given obstacle, Y (7, 8; 7) could immediately be determined by making use 
of equations (4.9.4) and (4.9.7). To find the image of a streamline, 


y (72027, i= § = Const, (4.9.8) 
in the physical plane, a sufficiently dense set of points (75-0:) ="), 2eaee, is- first 


found, corresponding to the streamline ¥ (7, 0; 7) = s. From the tables, the values 
are read for 


X) = X) (7, 0,; 149), Y= YO (%,, 0,; m9), (4.9.9) 


and by using (4.9.5), the image ¥ (x, y) = s of Y (7, 0; 79) = s is found in the physical 
plane. : 

The essence of this method is (as explained in previous sections) the separation 
of the quantities X”), Y, W, that are independent of the profile from the 6,’s 
which depend on it. For the former, a set of tables can be prepared once and for all. 
The problem of determining the latter recurs with each profile. This usually is a 
relatively simple matter, in the case of incompressible flow. In the compressible case, 
the problem is much more involved, but in most cases of interest, the b,’s found for 
the incompressible case may be used in the compressible case as a first approximation 
to the exact b,’s. Since the published tables were calculated for solutions having 
forms a little different from those cited above, we shall present these alternative 
formulas. The equation corresponding to (4.5.11) now has the form 


ee ap Om Vig Oe Cl = 212 )e Oe _®, (4.9.10) 
with 7 =H, and the complex solutions of (4.9.10) are given in the form: 


teil 


P* [f (@)] =| H* (n, 05 m0) #2/2) (1 — #)1 [atl — YY, (4.9.11) 


or 
Old 2) So(— 4) 4X92) HD Ue) en + DL @ + 1)) Pm) 9g @)}. 
(4.9.11b) 


where z = 7 +-7 6, and g (z) is an analytic function of z. The function #* may be 
taken in the form 


E* (y, 03 9) = Sy (— n)~* F{1/6, 5/6, 1/2, X}, (4.9.12) 

where F(x, 6, y, X) is the hypergeometric function and the other symbols used are 
xX =#[e(— ny) +20 — nol/2¢ (— ny”, (4.9.13) 

© = (Da) (ik = 1)/2 eee eo — DY (4.9.13a) 


The quantity S, is given by (4.6.3a), and the polynomials P™ (z) are defined by 
the recursion formula and initial conditions: 
PY LOR 0.0 isi 22 2 Ft PO a0 6 (499.14) 
P (0) = 0, (4.9. 14a) 
with 
F* = (5/36) (1/4 ¢) (— )-?. (4.9.15) 
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The functions g‘” (z) may be written as 


1 


zz Lene ns t=2 ie 
gz) =)... 4g @n)dz = ((— I*"/(n—1)!] | @— tH)” g(t)dt. (4.9.16) 
00 0 t=0 


The stream function is defined as the imaginary part, i.e, Y=IJIm[p*(g)]. A 
form similar to that one in (4.9.4) may be given, namely 


p* [9 2) = 2) bp {@'— 9) }, (4.9.17) 

where oi. 
p* {(z — yo)’ —"!2} = B® [p* (2), p* (2); no] [p* (2) — mo", (49-18) 
p* (2) = p* (4 +t 6) =— €(— 1)" 170 = 3 (4.9.18a) 
B® [p* (2), p* (2); no] = B® (3, 35 to] = So(— 0)“ F [1/6, 5/6, » + 1/2, X], (4.9.18b) 
X = (3 — m)/[(— 2¢) (— n)1. (49 486) 


It is seen, then, that in the case of the stream function, the transition from the 
incompressible to the compressible case is made by multiplying each term of its power 
series expansion by a “‘compressibility factor’ B (3,3; 79) and by replacing the 
variable z—7+i%6 by the variable 3 — —c(—y)/?+726. This latter change 
implies a distortion in the 7 direction mentioned above. Corresponding to expression 
(4.9.4) for the incompressible case, the stream function for a symmetric profile can 
be written in the form 


Y (n, 6; Mo) = Dd) 1b, P™ (n, 6; no), (4.9.19) 


where 
WO) — Re [B (4, 33 9) (8 — M0) —"} (4.9.20) 


The implicit representation of the streamline Y (7, 6; 79) = s in the physical plane 
is given by 


== D/ iby XO” (9, O53 Ho) YU = tO i Caneel (4.9.21) 
with 
XO = aA (By ee cos 6 — Y sin 6) dy + (C, B® cos 6 — ne sin 6) d6], 
(4.9.22) 
Yo = | A, [(— B, YP sin 6 — Y7” cos 6) dy + (C, ¥{ sin 6 + P') cos 6) d6], 
(4.9.23) 
A, = [oe exp (| /— 1) dy], (4.9. 24a) 
B, = (1 — M?)'ts [— I* (n) Ths, (4.9.24b) 
0, = (1 — mye [f—(y)), (4.9.24¢) 
L* (yea as (4.9.24d) 


In order to evaluate the “compressibility factors’ B (yj, 6; 9), suitable power 
series representations must be found for the hypergeometric function appearing 
in (4.9.18b). It is well known that for |X| <1, F (a, B, y, X) may be represented 
by a hypergeometric series F (a, B, y, X) and for |X| >1 a combination of two 
hypergeometric series with arguments and parameters different from those appearing 
in the hypergeometric function may be used. In the case of the hypergeometric function 
appearing in (4.9.18b), the boundary separating the domains of convergence is 


[[—¢ (— 9)" = Ho + 4 Of 2e (— n)]| = 1. (4.9.25) 
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Condition (4.9.25) may be presented in the form: 


3c x? + + 2c(— nn, + 2=0, for 7 <0, (4.9. 26a) 
3 0% 75 — @ — 2c¢(— nyt ny —ne = 0, for 7 > 0. (4.9.26b) 

In these domains, the following representations are used respectively: 
F (1/6, 5/6, » + 1/2, X) = F (1/6, 5/6, » + 1/2, X), (4.9.27) 

F (1/6, 5/6, » + 1/2, X) = A, (1/X)"* F (1/6, 2/3 — v, 1/3, 1)X) + 

+ B, (1/X)"* F (5/6, 4/3 — v, 5/3, 1/X), (4.9.28) 
where the quantity X is given by (4.9.18c). The A,’s and B,’s denote the quantities: 
A, = l'(— 2/3) (v + 1/2)/[F (— 1/3 — ») P (5/6), (4.9.29a) 
B, = I (v + 1/2) IF (2/3)/(L (1/3 — ») LF (1/6)]. (4.9.29b) 


It has been shown above how it is possible to extend the methods used in the case 
of an incompressible fluid to the case of a compressible one. The analogues of 
equations (4.9.2), (4.9.4), (4.9.5), and (4.9.6) of the incompressible flow were 
obtained, namely equations (4.9.17), (4.9.19), (4.9.21), (4.9.22) and (4.9.23). 
These permit the setting up of tables for a compressible fluid of the type that was 
outlined for an incompressible fluid. 

In the literature one may find two types of tables. Type I depends upon a para- 
meter 7, and has to be prepared for a number of values of 7, say %o,, (v= 1,...,). 
In this case, one must resort to interpolation methods in order to obtain streamlines 
for intermediary values of 19: 1,4 <1 <o,»41- A table of Y (7, 6; 79) (4.9.20) 
for a specific value of 7» (jy) = — 0°02) is given in %. In the same paper, a table is 
also included for 

y” (n, 0; — 0°02) = I m {p* [(z — 0-02) ]}, (4.9.30) 


in view of possible later applications, as well as a conversion table from the variable 7 
to the local Mach number MM. 
The second type of tables which may be calculated are the tables containing the 


values of F (1/6, 5/6, » + 1/2, « +7). These values may be computed once and 


for all, as well as derivatives of F with respect to the argument for a number of 
values of « and #, where 


a = 1f2—yol[—2e(— mJ], B= O/[—2e(— mJ]. (4.9.31) 


By the use of tables of fractional powers and trigonometric functions, and the tables 
outlined above, the computation can now be made of (4.9.20), namely: 


wr) — Re {S,(— 7" [— ¢ (— 2 — no +2 OY”: P, (4.9.32) 
F = F (1/6, 5/6, » + 1/2, « +7). (4.9.32a) 


Similarly, the derivatives ah (yn, 9; 4) and a (yn, 9; 49) can be calculated. The 
stream functions in the hodograph plane are obtained by forming ¥ (, 9; 7) = 
=» ib, V™ (n, 9; Ho). The integrands of equations (4.9.22) and (4.9.23) may 
be evaluated now by using the tables for the derivatives of F. Since the integrands 
of (4.9.22) and (4.9.23) are perfect differentials, the value of the integral is 
independent of the path and it may conveniently be taken along a streamline. The 
transition to the physical plane is carried out by numerical integration of (4.9.22) 
and (4.9.23). 
Ingenieur-Archiv VII, 4. 25 
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In * one may find the application of the tables, mentioned above, to the calculation 
of the flow around an oval-shaped obstacle in the compressible flow domain 
corresponding to the flow around an ellipse in the incompressible flow domain. The 
example shows clearly the entire procedure. 


13 §, Bergman: On tables for the determination of transonic flow patterns. Hans Reissner 
Anniversary Volume: Contributions to Applied Mechanics. J. W. Edwards. 1949. 
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Die Laplace-Transformation und ihre Anwendung. Von P. Funk, H. Sagan und F. Selig. Mit 
18 Textabb., VII, 1068S. Wien: F. Deuticke. 1953. S 50.—. 


Die Laplace-Transformation erfreut sich bei den Technikern und Physikern groBer Beliebt- 
heit, liefert sie doch eine strenge Begriindung fiir den Heaviside- oder symbolischen Kalkiil. 
Diese Wertschatzung der Laplace-Transformation geht aber oft nicht mit einem wirklichen 
Verstiindnis der Theorie Hand in Hand, erfordert sie doch etwas weitergehende mathematische 
Kenntnisse, welche nicht allgemein gelaéufig sind bzw. nach dem Abgang von der Hochschule 
nicht mehr geléufig sind. Man kann auf diese Sachlage den Ausspruch eines Spaniers in bezug 
auf die deutsche Sprache anwenden: Es geht mir mit dieser Sprache wie mit meiner Frau. Ich 
kenne sie, ich liebe sie, aber ich beherrsche sie nicht. Die Literatur tiber die Laplace-Transformation 
ist sehr umfangreich. Es gibt darunter ausgezeichnete Werke, z. B. die bekannten Biicher von 
G. Doetsch, aber sie sind den Technikern oft zu umfangreich und zuviel mit Deltasparren ver- 
sehen. Andere Biicher sind wieder zu kursorisch und rezeptartig. Die Verfasser unternehmen 
nun die gewi8 nicht leichte Aufgabe, eine Hinfiihrung zu geben, welche die rechte Mitte halt. 
Das Buch geht auf Vorlesungen zuriick, die der erste Verfasser verschiedentlich tiber diesen Gegen- 
stand gehalten hat und die sich vor allem an Techniker wandten. Die mathematischen Ent- 
wicklungen werden nur so weit gegeben, da der Gedankengang ersichtlich ist. Fur alle Details 
wird auf die systematischen Lehrbticher verwiesen. Zur Erlauterung werden ausgewahlte Beispiele 
von nichttrivialem Charakter gebracht, welche dem Referenten als das Wichtigste an diesem 
Werk erscheinen. Hs sollen nun einige dieser Beispiele angefiihrt werden: Transformator, Abelsche 
Integralgleichung und ihre Anwendung auf die Bestimmung der Dichte der Erde aus der Fort- 
pflanzungsgeschwindigkeit eines Erdbebens nach G. Herglotz und E. Wiechert, Warme- 
leitungsgleichung (Bremsung eines Rades), Thomson-Kabel. Der letzte Paragraph beschaftigt 
sich mit der wichtigen Theorie der Regelvorgange. Er erscheint mir gerade an den entscheidenden 
Stellen zu kursorisch und erfordert wohl fiir eine zweite Auflage eine Umarbeitung. 

Der erste Verfasser hat dem Buch noch ein Nachwort hinzugefiigt, welches eine schéne Ubersicht 
uber die historische Entwicklung gibt und sich mit der Persénlichkeit und dem Werk von 
Heaviside liebevoll und kritisch auseinandersetzt. Hervorgehoben seien noch die (knappen 
und vielleicht nur von den Kennern ganz zu wiirdigenden) Bemerkungen, welche die Theorie 
in die allgemeine Theorie der Funktionale einordnen und ihr so erst die richtige Beleuchtung 
geben. Die Verfasser haben aus Griinden, welche sie im Vorwort darlegen, keine Korrespondenz- 
tabelle in ihrem Werk aufgenommen. Trotzdem glaubt der Referent, da8 man bei einer zweiten 
Auflage dem Buch eine Tabelle dieser Art, welche wenigstens die Korrespondenzen, die in dem 
Buch behandelt werden, umfa8t und eine Formelsammlung hinzufiigen sollte. Dies wiirde den 
Wert des Buches vergr6Bern. 

Der Verlag hat das Buch in guter Ausstattung herausgebracht. E., Hlawka, Wien. 
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Trotz der unzweifelhaften Wichtigkeit, die durch Temperaturinderungen hervorgerufene Span- 
nungen in der Praxis besitzen, gibt es bisher keine zusammenfassende Darstellung ihrer Theorie. 
Wohl finden sich in verschiedenen Lehrbiichern der Elastizitaétstheorie vereinzelte Hinweise 
und Beispiele; der iiberwiegende Teil der Untersuchungen ist aber in zahlreichen Einzel- 
veréffentlichungen iiber viele Zeitschriften verstreut. Wenn daher in diesem Buch der Versuch 
unternommen wird, eine einheitliche und zusammenfassende Theorie der Warmespannungen 
zu geben, so glauben die beiden Autoren, durch ihre Arbeit einem wirklichen Bediirfnis 
entgegengekommen zu sein. 

Im vorliegenden Buch werden ausschlieBlich stationare Temperaturfelder unter der Annahme 
der klassischen Theorie der Elastizitat, also eines linearen Zusammenhanges zwische1: Spannungen 
und Verzerrungen, untersucht, wobei die Unabhangigkeit der elastischen und thermischen 
Materialkonstanten von der Temperatur vorausgesetzt ist. Wenn auch hierdurch die Anwendung 
der Theorie auf verhaltnismaéBig kleine Temperaturaénderungen beschrinkt bleibt, so umfaBt 
sie doch einen groBen Problemkreis und liefert die Voraussetzungen fiir die Lésung vieler 
praktischer Aufgaben. Eine weitere Arbeit, die auch die nichtstationéren Temperaturfelder 
beriicksichtigt und auf elastisch-plastische Spannungszustande eingeht, ist geplant. 
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Das vorliegende Buch erméglicht dem Planer und Konstrukteur, die optimale Liésung bei der 
Projektierung einer Pumpwerksanlage zu finden. Es vermittelt ihm Kenntnisse iiber die grund- 
sitzliche Arbeitsweise von Kreiselpumpen und ihre Betriebseigenschaften, die verschiedenen 
Bau- und Betriebsarten von Pumpenanlagen und die hydraulischen und elektrischen Steuer- 
und Schaltméglichkeiten fir Elektrokreiselpumpen. Neben dem hydraulischen Arbeitsbild wird 
immer der elektrische Schaltplan und eine schematische Darstellung des Pumpwerks gebracht. 
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